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PREFACE 


OBJECTIVES 

The primary objective was to evaluate potential EOS image-position- 
ing methods using various combinations of satellite data and ground 
control points* Secondary objectives were to assess the merits of the 
RBV camera for the EOS mission; verify the assumed user requirements 
for image resolution, coverage, registration, and positioning; model 
the imaging geometry of a conical-scanning Thematic Mapper to assist in 
the evaluation of image-positioning methods; and comment on sensor and 
ground-processing features that, if developed, might enhance system 
performance* 

SCOPE 


The study made use of information obtained from literature sur- 
veys, discussions with NASA personnel monitoring the program, and dis- 
cussions with potential users of EOS data* The analyses used to esti- 
mate potential accuracies of different image-positioning methods were 
based on photogrammetric experience and state-of-the-art data. The 
discussion covers image sensor design factors, position and attitude 
considerations pertaining to external EOS geometry, and geometric error 
analysis of positioning methods* 

CONCLUSIONS AND RECOMMENDATIONS 

EOS attitude control system design must be considered as a part 
of the image positioning problem for EOS, along with image sensor de- 
sign and ground image processing system design. With suitable effi- 
ciency for ground control point selection and matching activities during 
data processing, extensive reliance should be placed on use of ground 
control points for positioning the images obtained for EOS and similar 
programs . 
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SECTION 1 

INTRODUCTION AND SUMMARY 


1.1 BACKGROUND AND GENERAL COMMENTS 

This report describes the results of a study carried out by 
Bendix Research Laboratories (BRL) , Southfield, Michigan, from April 
through August, 1972, under NASA/GSFC contract NAS5-21727. The work 
was directed toward mapping accuracy considerations for the Earth 
Observatory Satellite (EOS) program. Many conclusions are also be- 
lieved applicable to earth-resource imaging satellite programs in 
general; the primary goal assumed for such programs is the generation 
of geometrically corrected and registered images, positioned with 
respect to the earth’s surface. The three EOS sensor candidates of 
interest were the high-resolution imagers: the Thematic Mapper 

(Seven-Channel Scanning Radiometer), the Return Beam Vidicon camera 
(RBV) , and the High-Resolution Pointable Imager (HRPI) . 

The primary study tasks were: 

(1) Analyze the relative merits of image-positioning methods 
using satellite-data only, ground control points only, 
and advantageous combinations of the two. 

(2) Estimate the positional accuracy that can be obtained for 
each of the methods . 

Additional tasks included: 

(3) Assess the merits of the RBV camera for the EOS mission. 

(4) Verify the assumed user requirements for image resolution, 
coverage, registration, and positioning. 

(5) Model the imaging geometry of a conical-scanning Thematic 
mapper to assist the analyses of the primary tasks. 

(6) Comment on sensor and ground -processing features that, 
if developed, would further improve system performance. 

(7) Define attitude performance requirements in terms of 
allowable contribution to positional error. 

The need for EOS image-positioning on the earth arises from the 
uses planned for the images. Details of interest on EOS images will 
be compared with other kinds of information: ERTS images, larger- 

scale aerial images, maps, tabular, and textual information. The 
established earth system of latitude and longitude forms a common 
reference for all these comparisons. A uniform plane-coordinate re- 
ference such as the Universal Transverse Mercator grid is useful for 
some comparisons, and for measuring distance, direction, and area. 


1-1 



EOS image positions should be defined with respect to one or both of 
these coordinate systems. In addition to this absolute position, the 
17-day repeat cycle of the EOS images will be of considerable benefit 
in making temporal comparisons of phenomena from successive EOS images. 
Some analyses require comparison on a resolution-element basis . As 
a result, very accurate relative positioning of one EOS image to an 
earlier or later image of the same earth scene is sought. This form of 
relative positioning is called temporal registration here. Finally a 
single imaging observation of an earth scene is made by several sensors 
in several spectral bands. Comparison between images from different 
bands and sensors is very useful, once again on a resolution-element 
basis. Hence, accurate relative positioning between the resolution 
elements of different sensors and spectral bands is important. This 
form of positioning is called event registration . Since the EOS 
resolution elements will be as small as 20 meters, considerable care 
must be taken if temporal and event registration are to be done with 
the desired accuracies . 

It is possible to produce the EOS images with no effort toward 
providing position and registration. This task then would simply be 
passed to the user, who would have to make the necessary investment 
in equipment and personnel for his purposes. If there were to be only 
one user, this solution would be as effective as any other. Since, 
however, there are many different disciplines engaged in the use of 
earth-resource data, it is clearly desirable that the image data from 
EOS he positioned and registered to the maximum extent consistent with 
production requirements before it is released . 

Given that EOS image positioning and registration are to be per- 
formed, the next questions are, how and how well? This study is con- 
cerned with the answer to these questions. The study is based upon a 
meeting of two different areas of technology with a common goal: 
positioning all points in an image with respect to the earth's surface. 
Assuming the internal geometry of the imaging sensor is well known, 
six external parameters, three translational and three rotational, are 
needed to position the image with respect to the earth. 

One way of determining these parameters is with the use of satel- 
lite data, information derived completely independently of the image 
content. Customarily, satellite data consists of (1) the satellite 
ephemeris for the three translational parameters, (2) measured and 
telemetered spacecraft attitude angles, and (3) calibration or auto- 
collimation data to relate spacecraft attitude to image-sensor attitude. 
The space images obtained from meteorological satellites have been 
positioned exclusively from satellite data. A large body of knowledge 
has been developed in the performance of this positioning method, 
largely taken from missile and satellite guidance, control, and tracking 
technology ♦ 
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A second way of positioning an image is with the use of ground 
control information derived using only image information content. 
Ground-control data consist of (1) the known locations of a few earth- 
surface features, and (2) the measured locations of these features 
in an image. This method of positioning, with various permutations, 
has been used for years in the precise mapping of the earth f s surface 
from aerial photographs. The body of knowledge developed here has 
been taken largely from photogrammetry , surveying, and geodesy, dis- 
ciplines considerably less familiar to NASA mission planners than 
those associated with satellite-data positioning. 

The two image-positioning methods have been successfully combined 
recently in the Earth Resource Technology Satellite (ERTS) program. 
Satellite data initially are used to position all images obtained from 
the ERTS spacecraft. The attitude data are taken from horizon sen- 
sors and limit the positional accuracies obtained. Ground-control 
data are used to increase positioning accuracy for a selected number 
of images, and also to calibrate and monitor the offsets between the 
spacecraft attitude sensors and image sensors. (The ERTS combination 
should not be construed to mean that satellite-data methods are in- 
herently inaccurate; the Orbiting Astronomical Observatory (0A0) 
program incorporates a stellar tracker to achieve very high (inertial- 
reference) attitude accuracies.) 

From the ERTS experience, it is clear that both image-positioning 
methods have advantages and disadvantages. It is also apparent that 
a systems approach to the problem is preferable to one in which image 
positioning methods are incorporated after the rest of the program 
is complete. The question for EOS planning, then, is not, "which 
method is best for EOS? ", but, "how should the two positioning methods 
be combined to give the best approach for EOS?" The combinations 
selected for this study may not necessarily be the best, although some 
effort was made here to select only the most promising. More extensive 
interaction between the two technology areas is needed. In fact, the 
major conclusion of the present study has been the recognition that 
the photogrammetry/surveying/geodesy technology and the control/guidance/ 
tracking technology should work together to develop the best positioning 
scheme for EOS . Both areas have a specialized technical vocabulary 
and audience. It is not reasonable to expect a group well-versed in 
one area to master the other and at the same time conduct an adequate 
technical evaluation of a many-faceted problem. This conclusion is 
borne out not only by this report, but also in some of the references 
cited which were written from the opposite side of the technological 
fence . 


1.2 SUMMARY OF FOLLOWING SECTIONS 

Section 2 of the report presents a short description of the different 
positioning schemes that were considered in the study. In Section 3 the 
image sensors are discussed, in particular the internal geometry of the 
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sensors. This is an important subject, since the accuracy with which this 
geometry can be described is the limiting consideration in event registration 
Offset pointing is also considered, since it represents a departure from 
the ERTS geometry and carries with it certain geometric implications. Align- 
ment between attitude and image sensors is important for satellite-data 
positioning, and is discussed in Section 3.4.1. The value of the RBV camera 
as an EOS sensor is discussed in detail in Section 3.4.2. 

In Section 4, the central matters of satellite position determina- 
tion, attitude control and measurement, and ground-control point princi- 
ples and applications are discussed in some detail. 

Section 5 is concerned with the accuracy attainable for the different 
positioning methods. Error analyses are presented based on the best 
values presently available for the different components. Finally, Sec- 
tion 6 presents the significant conclusions of the study. In most cases 
these conclusions also are stated in the relevant technical discussion. 
Recommendations for additional work also are presented in this section, 
following as logical extensions of one or more of the conclusions. 

One problem in the discussion of internal sensor geometry was found 
to be a lack of knowledge about the geometry of a conical-scan Thematic 
Mapper such as that implemented in one design under consideration. 

Appendix A, written by T. A. Eppes , describes the analysis performed to 
better understand the detailed geometry of this sensor. In the course of 
the analysis, a Fortran simulation of the scanner geometry was programmed 
and exercised. The program and the results of the testing are given in 
this appendix. 

Appendix B, written by R. J. Ouellette, presents the findings of 
a user requirements survey, performed to verify the ground-resolution 
needs of EOS image-data users. The survey summarizes many conversations 
with technical people in several Government agencies . The informality 
of these conversations was useful in obtaining user opinions under differ- 
ent circumstances from those used to compile published statements of agency 
needs, plans, and programs. 
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SECTION 2 

POSITIONING TECHNIQUES 


Two kinds of exterior orientation information are needed for posi- 
tioning the EOS images: the three locational components of image-sensor 

position and the three angular components of attitude* For the RBV 
camera, this information is needed only for the moment of exposure. How- 
ever, the Thematic Mapper and the HRPI require position and attitude data 
continuously over many minutes of uninterrupted image collection. The 
goals for position and attitude accuracies are commensurate with the 
highest-resolution image sensor. The 20-meter ground resolution of the 
HRPI for EOS is equivalent to 0.0012 degree.* These numbers indicate the 
magnitudes of interest for EOS position and attitude measurement. 

Five components will be considered here in developing image- 
positioning techniques: 

(1) Best-fit satellite ephemeris, for position 

(2) Ground control points (GCP) for position and/or attitude 

(3) Stellar absolute inertial attitude sensor 

(4) Gyro short-term attitude measurement assembly 

(5) Horizon-scanning earth-oriented attitude sensor 

These components can be combined in several different ways to obtain the 
six image-positioning parameters required. Either (1) or (2) must pre- 
sent in any scheme to provide the three locational parameters. For the 
angular values, the situation is not so simple. Depending on the design 
of the attitude control system, no attitude information, a limited amount, 
or a great deal may be needed for positioning images. 

Five main combinations of the components listed above were selected 
for detailed consideration. They are summarized in Table 2-1 and described 
in the remainder of this section. Extensive error analysis of the methods 
is postponed until Section 5, after additional discussion of the differ- 
ent components. 


During this study, angles have been seen specified in arc seconds, 
decimal parts of a degree, radians, and "milli-earth-rotation units". 

In this section, degrees and decimal parts of a degree are used through- 
out for uniformity, simply because the attitude control system specifi- 
cation is so expressed in reference (1) . (Some standardization of angu- 
lar measurement units seems as desirable for reports as the standardi- 
zation of linear units . ) 
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Method 


Ephemeris 
Position Data 


Frequent 


GCP 


SIMS-E 

X 

GCP-E 

X 

SIMS- EG 

1 

X (relative) 

GIMS-E-l 

X 

GIMS-E-2 

X 

GIMS-E-3 

X 


HIMS-E 


X 


Components 


Stellar 

Sensor 

X 



Good fit of attitude angles to 
low-order curve within an image* 


Good relative ephemeris 
accuracy. 


X 


X 


Smoothly varying ephemeris 
error; even gyro drifts. 

X 


Good relative ephemeris accuracy; 
low and even gyro drifts. 






































2.1 STELLAR INERTIAL MEASUREMENT SYSTEM (SIMS) AND EPHEMERIS 

This method, referred to hereafter as SIMS-E, represents a pure 
satellite-data positioning scheme, using components (1), (3), and (4) 
from Table 2-1. Ephemeris data define the three translational parameters 
of the image sensor. The gyros provide short-term attitude information. 

If the gyros are gimballed, this information is provided directly in 
some form of local-vertical attitude angles. If the gyros are strapped 
down, the rates must be integrated to give inertial attitude angles 
which are further processed from the ephemeris to give local-vertical 
attitudes. In either case, the gyros are subject to drifts and other 
anomalies which make them unreliable over some period of time which is 
short with respect to the operational life of the EOS satellite. The 
stellar sensor provides absolute inertial-attitude fixes at the time 
intervals required to update the gyro data. 

This method further requires that the angular relationships between 
the star sensor, the gyro housing or platform, and the image sensor be 
defined to the same or higher accuracies than those needed for the stellar 
inertial— attitude fixes. This is called the attitude-transfer require- 
ment here. The simplest conceptual way of meeting this requirement is 
to mount the three types of sensors directly to one another. This is 
possible for only one image sensor, and might well severely impact the 
design difficulties for both SIMS and the imager. If direct mounting 
is not feasible, SIMS attitude values can be referred to a calibration 
block mounted to one of the image sensors. This solution shifts the 
attitude-transfer problem to the sensor designer to ensure rigidity 
and parallelism between the block and the internal geometry of the imager. 
It also is suitable for only one of the image sensors, probably the one 
with highest resolution, HRPI. 

The next simplest solution is to mount SIMS and the image sensors 
to a common rigid structure; the spacecraft would be the common struc- 
ture, but it is difficult to thus maintain a calibrated attitude relation- 
ship to the one or two arc seconds needed. The remaining alternative is 
some form of "auxiliary-auxiliary" or tertiary attitude sensor which 
monitors the attitude between the image sensor and the stellar/gyro 
package. In reference (1)* a three-axis autocollimator is suggested for 
the tertiary sensor. This may take several forms, and is discussed in 
more detail in Section 3.4*1. 

The SIMS-E method illustrates the strengths and weaknesses of the 
satellite-data technique for image-positioning: 

• Image content is not considered at all for positioning. 

• Ground processing of image data is completely automatic. 

• Reliability of the three types of attitude sensors is of 
vital concern. 


^References are listed in order of mention at the end of the report. 
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• Attitude transfer is needed to relate attitude and image 
sensors, 

• Random ephemeris errors impose the ultimate accuracy limit, 
even with zero errors in attitude, attitude-transfer, and 
internal image geometry, 

• Development costs are largely those associated with space- 
qualified reliable attitude and attitude-transfer sensors. 

• More spacecraft weight and power are needed beyond that nor- 
mally required for attitude control, 

2.2 GROUND CONTROL POINTS (GCP) AND EPHEMERIS, METHOD GCP-E 

This method represents a nearly complete ground-control positioning 
technique, using components (1) and (2) in Table 2-1. Ephemeris data are 
used for two of the three parameters, along- and across-track position. 
Ground control points are used for the other four unknown parameters, and 
also for other parameters associated with mean attitude rates and accel- 
erations. The GCP positioning technique is such that errors in the 
ephemeris values are compensated by two of the angular parameters with- 
out significant error; the technique is described in detail in Section 4.3. 
The ephemeris data are not strictly necessary for the solution, but are 
normally available for the satellite anyway, and provide a modest statis- 
tical improvement in positioning accuracy. 

This method is the only one of the five which does not require or 
use telemetered attitude data from gyros on the spacecraft. To achieve 
this freedom, some constraints are placed on the attitude behavior of 
the spacecraft, other than those originally specified in reference (1): 
attitude behavior must be known in advance to be such that deviations 
from some simple behavior model will not be in error over a 30-second 
period by more than 0.0005 to 0.001 degree in each axis. The magnitudes 
of the angles and rates are not of concern* Details are given in Sec- 
tion 4.2. 

Method GCP-E shows the typical characteristics of the GCP method 
of positioning, constrained somewhat in this case by the non-framing 
EOS sensors: 

• The primary component used for image positioning is the infor- 
mation within the image. 

• Ground processing of image data requires some manual operations, 
supported by equipment and procedures designed to obtain highest 
production# 

• The only spacecraft sensor of concern is the image sensor itself. 

• Attitude transfer is not required, 

• The accuracy limit is determined by image resolution and internal 
image geometry, and by the small irregular deviations in attitude 
from a smooth low-order curve. 
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• Development costs are those associated with an efficiently 
operating ground data-processing facility, including some 
manual steps. 

• No extra spacecraft weight and power are needed. 

2.3 STELLAR INERTIAL MEASUREMENT SYSTEM PLUS GROUND CONTROL POINTS, 

METHOD SIMS -EG 

This method is the first of two hybrid positioning schemes which 
were selected for consideration in the study. Components (3) and (4) 
are used from Table 2-1, supported by shared limited use of components 
(1) and (2). Note from the table that the GCP are used infrequently, 
for one or perhaps a few images per orbital pass, rather than in each 
image to be positioned. They serve here primarily to detect and remove 
bias in the ephemeris. This bias is taken as representative, and applied 
ahead and back along the orbital arc to all other images to be positioned. 
If the spacecraft structure is suitably stable over a single imaging pass, 
the GCP will also solve the attitude-transfer problem. 

There are several advantages to this sparing use of GCP over method 
GCP-E. First of all, the many fewer GCP to be selected and correlated 
mean much less difficulty in a high-production ground processing facility. 
This is especially true if a completely digital image processing scheme 
is to be used; the lengthy digital image correlation needs are much less 
significant than in a method where each image requires correlation of 
several GCP. Second, the GCPneedbe selected only in areas with highest 
image quality. Third, the ephemeris error effects are greatly reduced 
from what they were for method SIMS-E. Fourth, the need for some GCP 
processing also gives a back-up capability for use in the event of 
spacecraft attitude-sensor failure. Finally, attitude-transfer between 
attitude and image sensors can be achieved using GCP instead of pre- 
launch calibration. 

Unfortunately, some of the negative aspects of methods SIMS-E and 
GCP-E also are incorporated in the hybrid SIMS-EG. Development costs 
are high, since both precision spacecraft attitude sensors and a ground 
processing facility for GCP are needed. Also, as is shown in Section 5, 
some of the errors of each method are added in this hybrid method, so 
that not much of an increase in accuracy is obtained. More spacecraft 
weight and power still are needed for attitude measurements than would 
be required only for attitude control. The reliability of the stellar 
and gyro attitude devices is of concern, too, although the GCP capa- 
bility does provide some backup. 

Table 2-1 shows that an additional condition is imposed on the use 
of method SIMS-EG: good relative ephemeris accuracy. Over orbital arcs 

of a few thousand kilometers, relative accuracy should be high enough 
(see Section 4.1) to make consideration of SIMS-EG worthwhile. 
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2.4 GROUND CONTROL POINTS PLUS GYRO DATA PLUS EPHEMERIS, 

METHOD GIMS-E 

This hybrid positioning method removes the attitude-behavior con- 
dition imposed on method GCP-E by incorporating telemetered gyro data. 
There are three variations of this method to be considered, with a de- 
crease in each variation of the amount of GCP needed. 

The first, GIMS-E-1 uses ephemeris data for along-track and across- 
track position. Gyro data provide high-frequency relative attitude data. 
GCP define altitude and absolute attitude values for the three rotational 
axes. Gyro drifts are considered here to be less than 0.001 degree over 
15 seconds. If this is not the case, mean attitude rates must be carried 
as three additional unknowns in the GCP image-positioning solution. This 
means either more GCP must be used to retain a given accuracy, or lower 
positioning accuracy for the same number of GCP per image. 

The second variation is GIMS-E-2. In this case, GCP are not used 
in every image for positioning, but only at intervals along an Imaging 
pass. Between GCP images, the attitudes and positions are interpolated, 
with bias-adjusted gyro and ephemeris data used for attitude and position, 
respectively. The variation requires fewer GCP operations than GIMS-E-1 
and preserves much of the accuracy, but requires better gyro character- 
istics and low random error in the ephemeris. 

The third variation is GIMS-E-3, Here, extrapolation is used in- 
stead of interpolation. GCP are used for a single image, together with 
ephemeris and gyro data. Values are obtained for the usual six parame- 
ters plus three mean attitude rates. These values are extrapolated 
forward and backward along the orbit using ephemeris differences and the 
angle- and rate-adjusted gyro data. Even fewer GCP are required here, 
in exchange for even better-defined gyro angles and ephemeris. 

The three variations differ only in their trade of GCP frequency 
for ephemeris and gyro quality. Thus, they share these characteristics: 

• Image content is used only to the extent necessary to over- 
come gyro and ephemeris deficiencies. 

• Ground processing of image data requires some manual opera- 
tions, but fewer than method GCP-E, and provides back-up for 
gyro failure. 

• Reliability of the gyro device on the spacecraft is of concern. 

• Attitude transfer is not required, if the attitude between the 
gyro reference and the different sensors remains constant 
between GCP updates. 

• Ephemeris and gyro error effects can be traded with frequency 
of GCP use. 

• Costs are in the development of a high-quality gyro package, 
as well as manually aided ground processing. 
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• Some spacecraft weight and power are needed* beyond that 
normally required for attitude control. 

The variations of GIMS-E place somewhat greater reliance on GCP than does 
method SIMS-G, with a consequent decrease in concern over spacecraft com- 
ponents and an increase in the required manual processing. The accuracy 
is somewhat improved from SIMS-G as well, at least for variation GIMS-E-1. 
Accuracy for the other two variations can be traded with ephemeris/gyro 
quality and number of GCP needed. 

2.5 HORIZON SCANNER WITH GYRO PACKAGE AND EPHEMERIS DATA, 

METHOD HIMS-E 

This method is included to show back-up capability and a cost/ 
accuracy benchmark. Method HIMS-E is a pure satellite-data positioning 
scheme, using components (1), (4), and (5) from Table 2-1. The method 
is functionally the same as that used to position the ERTS bulk images, 
although gyro performance is considerably improved. HIMS-E operates 
in the same way as SIMS-E, except that the horizon sensor replaces the 
stellar mapper or tracker of SIMS-E for attitude updates to be applied 
to the gyro data. The horizon data for HIMS-E will result directly (after 
some modeled adjustments) in true local-vertical attitude. This is dif- 
ferent from the situation using SIMS-E, where the attitude derived is 
inertial and must be transformed to the local reference frame using 
accurate ephemeris data. The horizon sensor cannot approach the stel- 
lar sensor for attitude accuracy, although considerable improvements 
can be made in the raw data by proper oblateness-ef f ect removal. The 
characteristics of this method are qualitatively the same as for SIMS-E, 
except that (1) spacecraft-component costs, reliability, weight, and 
power are somewhat less of a problem, and (2) the horizon— sensor error 
would be expected to impose the accuracy limit, not ephemeris errors. 
Overall, this method is probably the least expensive, and it is certainly 
the least accurate. 


A 

It may be possible to use the ACS gyros for the required attitude 
information. 
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SECTION 3 

IMAGE SENSOR CONSIDERATIONS 


The following discussion is particularly concerned with the three 
high-resolution image sensors presently proposed for the EOS payload: 
the Thematic Mapper, the HRPI, and the RBV camera pair. Some of the ma- 
terial will be seen to apply also to image sensors and earth— resources 
missions in general. 

The EOS orbit and ground-swath coverage parameters are those derived 
in reference (1) : 

(a) Altitude = 980 km above equator, circular orbit 

(b) Period - 104.6 minutes 

(c) Inclination = 99+ degrees 

(d) Ground— swath width for high-resolution imagers = 185 km 
(100 nm) , with 10 percent overlap at equator 

(e) Revisit time = 17 days 

(f) Local time of descending node = 8:39 AM 
3.1 INTERNAL GEOMETRY 

3.1.1 Thematic Mapper 

In reference (1), a conical scan geometry is suggested for 
the Thematic Mapper, The conical scan results from an effort to increase 
duty cycle without severe weight penalty for a particular resolution. 

It is possible that this goal also can be met with an optical design 
which results In the more usual straight-line scan geometry, such as 
that employed by many aerial scanners. The simpler geometry may be pre- 
ferable for ease of processing or some other considerations. Still, 
there is nothing in the use of conical scan which inherently limits the 
geometric use of the data. 

The effect of one conical image-plane scanner implementa- 
tion is shown in Figure 3—1. In this figure, the cone is vertical. The 
effective imaging geometry is based upon the text and the schematic dia- 
gram contained in reference (1), except that six catoptric probes are 
assumed instead of eight, as in reference (3). The six probes divide 
the cone base circle into six equal sectors, and 48 degrees of each sec- 
tor are used for active scanning (80 percent duty cycle). Six high- 
resolution lines are swept out in each scan swath, with the detectors 
for a given spectral region arranged in a one by six array, as shown in 
the figure. To cover a strip 185 km wide on the earth with a 48-degree 
central angle of the base circle requires a cone vertex angle slightly 
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SPACECRAFT VELOCITY VECTOR 



Figure 3-1 - Vertical or Circular Conical Scan 


greater than 26 degrees. The chief characteristic of vertical-cone imag- 
ing is that all objects on the earth are scanned with the same nadir 
angle. This might appear to be advantageous for mul tispectral applica- 
tions, since the spectral reflectance characteristics of many features 
change with viewing angle. However, the angle formed by the sun, ground 
feature, and the scanner detector also is important to spectral reflec- 
tance, and this angle changes during a single scan swath. Thus the 
fixed-nadir-angle scanning may not be of significant value. 

In Figure 3-2, the scanner just described has been mounted 
in the spacecraft with a forward tilt of one-half the vertex angle. As 
a result, the arc formed on the earth by a single scan swath passes 
through the spacecraft nadir at the swath center, with increasing nadir 
angle toward the ends of the swath. This configuration supposedly results 
in a more "natural” viewing geometry, since it is more nearly that of a 
conventional scanner sweeping a straight line from side to side through 
the nadir (e.g. , the ERTS MSS scanner). A more tangible reason for tilt- 
ing the cone is to reduce the positional error caused by terrain relief 
variations. Figure 3-3 shows the projective displacement, dS, of a hill- 
top away from its orthogonal position. The displacement is equal to 
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h(S/H). For the vertical-cone EOS imager shown in Figure 3-1, the ratio 
of S/H is nearly one-fourth throughout each scan line, so that all points 
on the terrain are differentially displaced radially from the nadir by 
amounts equal to one-fourth their height differences* For the tilted- 
cone imager in Figure 3-2, no terrain relief displacements occur at the 
center of each scan line, where S/H is zero* Displacement effects in- 
crease away from scan center until, at the ends of scan, they are nearly 
as great as for the vertical-cone imager. But the near terrain displace- 
ments still are significantly less for the tilted-cone than for the 
vertical-cone imager. Such displacements are noticeable only in a posi- 
tional sense, with respect to the map locations of image features; for 
either of the scanning cone configurations, any terrain displacements 
will repeat on subsequent imaging passes. Thus temporal registration 
will not be affected. 

The instantaneous field of view (IFOV) for the Thematic 
Mapper was originally specified as 66 microradians for the six high- 
resolution channels described in reference (1). At the nominal EOS al- 
titude of 980 km, a single high-resolution scan swath will contain 6 x 
2,870 samples, and the swath thickness will be 388 m. As a result of 
the EOS Mission Review Group’s conclusions, reference (2), the original 
IFOV will be reduced. The present goal is 44 irticro radians. If this is 
achieved, a single scan swath will contain 6 x 4,300 samples, and the 
swath thickness will be 259 m. The rotation rate of the wheel containing 
the six sets of scanning probes will be about 2.7 rps for the 66-micro- 
radian IFOV and 4.1 rps if the 44-microradian goal is achieved. 

Details of the conical-scan image geometry, its interaction 
with earth-rotation effects, and the effects of attitude and altitude 
variations are described in Appendix A. 

3.1.2 High Resolution Fointable Imager (HRPI) 

This device presently is planned as a solid-state linear 
array of 1,800 detectors for each of 4 spectral bands, each detector 
with a 20 m square IFOV. The resulting geometry is that of a straight 
line 20 m wide and 36 km long, normal to the spacecraft velocity vector. 
The sensor can be rotated about the velocity vector to image any part 
of the 185-km wide Thematic Mapper coverage. If the Thematic Mapper is 
made capable of off-axis pointing, the HRPI may be required to follow; 
this requirement is not yet defined. 

The extremely simple internal geometry of the HRPI has been 
aptly compared to a pushbroom. A straight line is scanned, with a fixed 
spatial relationship between all the detectors forming the line. As the 
spacecraft moves forward, this line is "pushed" ahead, covering the area 
of interest in a continuous swath. 
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3.1.3 Return Beam Vidlcon Camera (RBV) 

The RBV camera is familiar from its application as a multi- 
spectral imager in the ERTS program. For EOS no multispectral capability 
is to be used. The camera is proposed as a single-band (panchromatic in 
one reference) sensor to obtain reference images of higher resolution 
and greater locational accuracy than those provided by the Thematic Mapper. 
The RBV images also are to "aid greatly in achieving the desired location 
accuracy" for Thematic Mapper images , reference (2). 

A pair of 2-inch RBV cameras have been recommended for EOS. 
Each has a focal length of 250 mm, twice that used for ERTS, and covers 
an area about 98 km on a side. The cameras must be tilted about 2.5 
degrees, one to the left and one to the right of the ground track. Ex- 
posure and readout alternate between the two cameras. 

A total of 4,500 television lines per image are to be scan- 
ned (the ERTS RBV tubes are scanning 4,125 lines), equivalent to a ground 
pixel (theoretical picture element) of 22 m on the earth. But image 
motion compensation will be required if effective resolution is to be at 
all commensurate with this pixel dimension . 

Since the RBV cameras are basically analog television image 
sensors, the image distortion can be expected to change during any exten- 
ded period. To enable the detection and correction of these random dis- 
tortions, each of the RBV cameras for EOS will contain a reseau, a grid 
of marks etched on the camera faceplate. The measured locations of the 
reseau marks provide geometric calibration data. The same marks will 
appear in the RBV images and can be used to remove most distortions in- 
troduced during the readout sequence. For the ERTS cameras, the reseau 
array contained 81 marks in a regularly spaced 9 by 9 array. Some modi- 
fications to this pattern are recommended below. 

3.2 GEOMETRIC DESIGN CONSIDERATIONS 

In this discussion, the specific EOS image sensors are not so much 
of interest in their own right, but as examples of the geometric problems 
and potential in the design of image sensors. A great deal of geometric 
diversity is possible in designing image sensors to be used in earth- 
resources spacecraft. This stems from the different requirements to be 
satisfied by the sensors. Typically, high duty cycle, low weight, low 
volume, and low power requirements are sought. Under some conditions, 
by permitting "unusual" imaging geometry — and any departure from a 
central perspective projection is sometimes considered unusual — some 
of the other goals can be more nearly achieved. This is the rationale, 
for example, guiding the conical-scan geometry being considered for the 
Thematic Mapper design. 

For image positioning and registration, there is nothing inherently 
"wrong” about incorporating unusual geometry in the design of an image 
sensor. However, it is important that the geometry be well known . Only 
with known internal geometry can the image data be effectively positioned . 
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Three methods can be considered for providing the image geometry 
information. The first is to correlate the image data point by point 
against a master image for which the geometry is known. The correlation 
process can be analog or digitial. The great advantage of this method 
is the lack of attention needed in image-sensor design. On the other 
hand, the method is quite slow (particularly when using digital methods), 
it requires continuous tracking to maintain the correlation process, and 
it will fail when adequate image detail is not present, references (27, 
28), or when images from spectrally different regions are being matched, 
reference (26). It is clearly unsuited to high-throughput situations. 

The second method is to provide advance calibration data on the 
image sensor geometry. This enables direct mathematical modeling of 
the sensor. The calibration data can be derived from grid images made 
in the laboratory, or from images collected over a test range containing 
a large array of ground control points. If this method is used, some 
care must be taken during image-sensor design to ensure geometric sta- 
bility. Stability does not necessarily mean geometric simplicity , al- 
though data processing may be speeded if the geometry is simple. In any 
case, data processing proceeds quite rapidly using advance calibration 
data. 


The third method of providing geometric information is to provide 
calibration data pertaining to each image. Typically, this is done by 
providing a reseau (grid) in the image. An intermediate amount of effort 
is required in both the sensor design and the image-data processing, 
compared with the first two methods. The image sensor must be geometri- 
cally "well-behaved" between adjoining reseau marks, so that some form 
of interpolation will adequately describe the image geometry. Also, the 
normal data processing operation must provide for the identification 
and measurement of the reseau marks in the image data, in order to pro- 
perly calibrate the internal geometry for each image. This is a time- 
consuming operation, although much less so than point-by-point image 
correlation over the entire format. 

In practice, the second and third methods often merge. The ERTS 
multispectral scanner (MSS) uses the second method, advance calibration. 
Even so, some analysis based on GCP residuals will be needed during the 
operational life of the MSS to detect possible changes in internal geo- 
metry. This analysis could be done easier using optional timing lines 
within each scan. The ERTS RBV camera array incorporates the third 
method of providing internal geometric information. Early ERTS data 
has shown that the random shifts in the internal geometry of the RBV 
cameras are smaller or slower than originally thought, reference (4). 

If this behavior is characteristic, the reseau data from one image 
can be applied to a number of subsequent images, thus becoming a form 
of advance calibration. 
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Several points merit special consideration by image-sensor designers 
for EOS-type missions: 

(1) An interrelationship exists between the attention paid to 
geometric properties during design and the attention that 
must be given to the geometric processing needed to pro- 
cess the image data obtained with the sensor. 

(2) Stability is the most important geometric design consid- 
eration; adequate geometric stability should be provided 
in such form that only infrequent off-line methods need 
be used to monitor and update internal image geometry. 

(3) Adherence to a particular geometric pattern for an imaging 
sensor is not always a necessity, provided the sensor gives 
the desired ground coverage and the geometry can be ade- 
quately described. This freedom in selecting geometry is 
greatest when extensive and complex processing must be 
applied to the image data for other reasons. 

(4) Sensors with moving components require special calibration 
of their interaction with the spacecraft structure; depend- 
ing on the frequencies involved, satellite attitude control 
may be degraded (Section 4.2), image-sensor geometry may 

be degraded, or image quality may be reduced. 

3.3 OFFSET POINTING 

One or more of the EOS image sensors is presently planned to have 
some capability for commanded offset pointing. Such capability has 
several advantages for data users, including avoidance of clouds or haze 
and the ability to image phenomena which are short-lived with respect 
to the normal 17-day repeat cycle. 

At the equator, to center the Thematic Mapper on the adjoining 
ground swath to the east or west of the ground swath currently beneath 
the spacecraft, a roll angle of 10.7 degrees is needed. To center on 
the second ground swath east or west requires a roll angle of 20,5 de- 
grees. The third swath requires 29 degrees of roll. These angles are 
maxima, decreasing with increasing latitude. 

For purposes of data management it is most desirable to constrain 
the offset pointing capability to integral ground swaths. This is no 
hardship for users, since the coverage will be complete, and will cor- 
relate well with the normal images for the same ground swaths obtained 
in the normal downward-pointing mode. Since the overlap between suc- 
cessive ground swaths is a function of latitude, the amount of angular 
offset to be provided for the sensors must also be variable as a function 
of latitude, if offset imaging is to cover integral ground swaths. 

Reference (1) shows that some misunderstanding of the offset point- 
ing requirements may exist. The spacecraft motion supplies one dimension 
of a scanner image. This motion is quite rapid with respect to the 
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time-scale for earth-resource phenomena and cloud coverage. For example, 
by pitching a scanner in the EOS orbit ahead by 30 degrees, an image can 
be obtained only 1.5 minutes sooner than it would have been obtained by 
leaving the scanner vertical , Hence, there is no need to provide for 
a pitch offset. Gross yaw offset is not required either, at least for 
scanners which sweep rapidly compared to the earth's rotation effect. 
(There may be some rationale for a very small yaw offset to be provided 
by the spacecraft attitude control system as a function of latitude, in 
order to eliminate very small gaps in scanner coverage; this is discussed 
in Appendix A.) As a result, only a single axis is needed for offset 
pointing . One-axis pointing is considerably easier to design and imple- 
ment than a three-axis offset. 

Control of the pointing angle is the next item for consideration. 
Control can be provided in fixed large increments, such as five degrees, 
or in much finer increments. There is no need for extremely fine angular 
control here, since the image sensors cover several degrees. On the other 
hand, as mentioned above, efficient data management (cataloging, indexing, 
processing requests, retrieving image data) of the thousands of images 
which will be obtained dictate that the pointing be confined to integral 
ground swaths, referred to normal repeat cycle. The pointing precision 
should be about the same, too. For ERTS, repeat precision is ten percent. 
If the same figure is applied to EOS, the pointing increments need be no 
finer than about 0,5 degree, A stepper motor applying angular increments 
of this size would provide ground swath repeatability of about 6 percent, 
and 128 steps would cover the 60-degree range needed to step over 3 
ground swaths to either side. To minimize the effects of disturbance 
torques on spacecraft attitude stability, it is necessary to implement 
the sensor offset pointing well in advance of any planned high-resolution 
image collection . The lead time depends on the attitude control system, 
the design of the pointing device, the moment of the pointing optics, 
and the placement in the spacecraft. 

The accuracy with which the offset angle must be measured depends 
on the method used to position the image. For GCP techniques, measure- 
ment is not strictly needed at all, although some gross indication may 
be helpful in locating GCP for correlation. Satellite-data positioning 
methods require that the offset angles be measured with the same high 
accuracy as the spacecraft attitude angles. It may be desirable to 
avoid this high accuracy measurement by using only GCP for positioning 
off-axis images, even though the normal downward-looking images are 
positioned primarily with satellite data. 

In regions of moderate terrain relief, some degradation in temporal 
registration accuracy must be expected when comparing off-axis images 
with the normal downward- looking images (Figure 3-4). The extent of 
misregistration caused by this effect depends on the height of terrain 




Figure 3-4 - Terrain Displacement from Offset Pointing 


and the differences in the two imaging angles, and is given approximately 
by 


ds 


(214) 

(980) 


N dh 


0.22 N dh 


where 

ds = misregistration 

N = number of ground swaths between the two images 
being registered 

dh = terrain elevation 

The terrain elevation dh is measured above mean terrain if GCP are used 
for positioning, or if spacecraft data are being used together with a 
world-wide mean terrain model. For spacecraft-data positioning without 
such a model dh is measured above mean sea level. 


3-9 



A positive way in which to consider terrain displacement is as a 
new capability for stereoscopic viewing. Even for a single ground-swath 
offset, an adequate stereo effect will be obtained. Terrain height deter- 
mination can be performed under these conditions, with a relative error 
of 


_ 980 _ 4.6 

h 214 N p N p 

where 

- standard error of relative height determination 
nip = standard error of position repeatability 

3.4 OTHER SENSOR CONSIDERATIONS 

3.4.1 Attitude Transfer 

Some positioning methods depend on auxiliary data to deter- 
mine the attitude of the imaging sensor. One of these methods may be 
used for EOS. If so, it is necessary to know the relative attitude be- 
tween the auxiliary attitude sensor and each of the high-resolution image 
sensors. This is the attitude transfer requirement discussed in Sec- 
tion 2. The rms error with which this attitude should be known is com- 
mensurate with the rms error expected from the attitude sensor. To keep 
the effects of error to beneath significance in a root-sum-square error 
sense, the positional contribution of error normally must be held to 
one-third of the root-sum-square positional error from all other causes. 
For EOS, the rms error goal being used here is one HRPI resolution ele- 
ment, or 20 meters. The corresponding combined effect of attitude- 
transfer error goal then is 6.7 m/980 km or 0.00039 degree, about one- 
second of arc. 


It is not necessary, or even desirable, to devote effort 
before launch to calibrating the attitude- transfer angles. They can be 
more meaningfully measured after launch using GCP in a single bias-removal 
calibration. However, some assurance of rigidity is needed so that 
these calibration values can be used with confidence. 

The most straightforward way of achieving such rigidity 
would be to mount the image and attitude sensors together. This may not 
be possible for many reasons* But if the image and attitude sensors are 
mounted separately to the spacecraft structure, it appears to be very 
difficult to design the spacecraft with adequate rigidity to maintain 
such a high alignment accuracy. This leads to a tertiary-sensor method 
such as that suggested in Section 7.7-1 of reference (1). With this 
method, an additional optical device is used to measure and record the 
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three-axis misalignment between the external attitude sensor (a star 
mapper in this case) and one or more of the high-resolution image sen- 
sors* The alignment data obtained in this way would be telemetered to 
the ground receiving station along with the attitude-sensor data. This 
solution gives the greatest freedom to spacecraft design and to the place- 
ment in the spacecraft of the different sensor packages, and is undoubt- 
edly feasible. However, it places another set of unknowns in the posi- 
tioning problem, requiring additional reliability, cost, and care during 
system design to ensure adequate performance. It may be possible to com- 
bine the attitude alignment into either the auxiliary attitude sensor 
or the high-resolution image sensors, to either automatically maintain 
some fixed three-axis relationship or else record and compensate for 
deviations. Unfortunately, this solution simply shifts the problem to 
the image-sensor or attitude-sensor designers. 

It must be noted again that a positioning method for which 
GGP are used with each image is not affected by these considerations. 

A gyro package may be needed for the scanning imagers, but such devices 
provide only attitude differences , over periods of time which are short 
enough to assume a constant angular alignment between gyro package and 
image sensor. 

A second kind of sensor alignment concerns the different 
image sensors alone. To image the same area, these sensors must be 
boresighted before launch. Sometimes an array of imagers is intended to 
operate as a single composite image sensor; the array of three raultispec- 
trally filtered REV cameras used for ERTS is an example. An inordinate 
amount of effort appears to be spent sometimes in the prelaunch bore- 
sighting of image sensors. For the relatively large areas typically 
being imaged, a misalignment of 0.1 degree is certainly permissible. 

This can be achieved with little effort. In general, the additional ef- 
fort spent aligning image sensors more accurately could be more profit- 
ably devoted to designing ways to maintain constant alignment after 
launch. 

3.4.2 Framing Camera Utility 

In reference (2) , the addition of the RBV cameras to the 
original EOS sensor complement was recommended for two reasons. The 
first is to provide higher spatial resolution than the Thematic Mapper 
(20 m for high-contrast and 50 m for low-contrast targets is quoted in 
reference (2)), with a corresponding capability for high location ac- 
curacy. The second is to somehow "relate the spectral and spatial data" 
of the Thematic Mapper, and so "aid greatly" in accurately locating 
Thematic Mapper image data (quotes are from reference (2)). Thus the 
advantages cited are (1) the effective resolution of the RBV camera 
versus other sensors, (2) the attainable positional accuracy for the 
RBV images, and (3) the ability to effectively use the RBV images as 
locational references for Thematic Mapper images. These three points 
will be considered in order. 
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First, the resolution advantage. Resolution is not a pri- 
mary concern in this study, although it is unavoidably involved in some 
considerations. Resolution is not simply given by the theoretical pixel 
(picture element) size of the RBV camera, any more than by the instantan- 
eous field of view (IFOV) of a scanner. Effective signal-to-noise ratio 
also must be considered. Initial ERTS images already have shown that the 
net image quality of the scanner with its 75 m IFOV is higher than that 
of the RBV camera with its 45 m theoretical pixel. Part of the reason 
for this may be the shading and mottling in the RBV images. A second point 
in connection with RBV resolution concerns image motion. During a nominal 
0.010-second exposure, the camera moves over 60 m. Image motion already 
is blurring the ERTS 45-meter theoretical pixel size. There is no point 
in doubling the ERTS-camera focal length to obtain a theoretical 20-meter 
ground pixel unless some form of image-motion compensation is provided. 

The second advantage cited for the RBV is that of locational 
accuracy. A framing camera is not necessarily of high geometric accuracy 
simply because it is a framing device. This is the very reason why sim- 
ple amateur cameras are not used for photogrammetric mapping. The RBV 
camera is geometrically a framing sensor, but the reseau on the faceplate 
testifies to the need for extensive and continuous monitoring of the 
geometric stability. During data processing, each of the large array of 
reseau marks must be located and accurately measured, even in the pre- 
sence of background imagery. This task cannot but add appreciably to 
the total data processing time necessary. At least some of the reseau 
marks must be checked in each image, ref erence(4) , or series of images, 
to look for shifts in the distortion pattern. In fairness, early ERTS 
RBV images show rather good short-term stability of the image distortions. 

The third stated advantage of the RBV camera concerns its 
ability to serve as a locational reference for Thematic Mapper images. 

The locational accuracy of the Thematic Mapper is discussed separately 
in Section 5. It is not in need of an RBV camera to be improved. For 
the present discussion, however, the accuracy of the Thematic Mapper 
image is not so much of interest as the mechanics of using one image to 
position another. Appendix G-3 of reference (2) recommends investigation 
of image correlation techniques to spatially relate the RBV and other 
images. Automatic image correlation techniques, both analog and digital, 
are well known, references (26, 27, 28, 29, 30), both for use in con- 
touring and profiling the earth* s surface and for simple image matching. 

The ERTS precision processing subsystem incorporates a modification of 
this electronic image correlation experience to precisely lock on sel- 
ected ground control points. It must be emphasized, however, that there 
is a tremendous difference between correlating a few small areas specially 
chosen for their pattern and contrast and correlating continuously over 
the entire image. The problem is increased when attempting correlation 
between images from different spectral regions, reference (26). Experi- 
ence gained at these laboratories and elsewhere has shorn that image 
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correlation is a useful tool when used with discretion, human judgment, 
and manual back-up. This is true for digital as well as analog or hybrid 
applications. But correlation of all image details is extremely time- 
consuming; its wholesale application has no place in a high-production 
environment, and should not be considered further. 

Summarizing the discussion above, 

• Resolution: The RBV camera would need image-motion 

compensation and some other improvements before any 
comparison of net image quality is meaningful. 

• Geometric accuracy: The reseaus require lengthy 

processing and periodic checking by methods ill-suited 
to automatic operations in order to maintain geometric 
quality; the quality is not demonstrably higher than 
that of the Thematic Mapper. 

® Locational reference: The image correlation process 

is far too time-consuming and unreliable for high- 
production application in matching Thematic Mapper 
image detail to RBV images. 

In conclusion, then, an RBV framing camera on EOS will not affect Thematic 
Mapper positioning or registration . With improvement, the RBV may have 
merits of its own in terms of resolution or as a backup sensor. Still, 
there may be other imagers with similar attributes and they merit equal 
consideration, in which cost, weight, power, and data processing complex- 
ity are balanced against usefulness of the images. Other use may be 
preferable for the weight and power entailed; an on-board data-f iltering 
or compression package to cut transmission data rates is one possibility. 

A backup Thematic Mapper or Extended HRPI are others. 

If, in spite of these arguments, the RBV camera is approved 
as an EOS sensor, some form of image-motion compensation device must be 
provided for the camera . In addition, the existing reseau pattern should 
be changed as follows: 

(1) Reduce thickness of reseau crosses by at least 50 percent. 

(2) Scribe crosses as X T s rather than + f s, with respect 
to raster. 

(3) Change spacing so reseau density is higher at edges 
and corners, but total number of reseaus remains 
the same. 

3.4.3 Data Compression 

Although it is not a subject for the present study, the 
high data rates associated with EOS cannot help but impress the observer. 
In the next few years, users can be expected to demand even greater spa- 
tial resolution, further increasing the data rate. One way of limiting 
the data rate is through onboard processing and filtering of meaningful 
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information. This approach is already being pursued. Another method 
is to sense and/or transmit data for only those spectral regions and 
resolutions currently of interest. There still appears to be some di- 
versity of opinion as to the number of spectral bands necessary to detect 
the various phenomena under investigation, and the necessary amount of 
sensitivity within the bands. The numbers of bands needed, their wave- 
lengths, and the resolution within the bands can be expected to change 
as functions of experience, data application, time of year, and geographic 
region, to name but a few variables. Thus ground-selected spectral re- 
gions and resolutions may satisfy the spectral needs of users, while 
permitting a more reasonable data rate. It has been noted from several 
user surveys, including the one carried out as a part of this study, that 
the need for multiple spectral bands to some extent complements the need 
for high spatial resolution. Thus, some advantages might be seen for an 
imager which could provide the same data rate by high resolution using 
a single panchromatic channel, or by decreased resolution but with several 
spectral channels. 
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SECTION 4 

POSITION AND ATTITUDE 


This section describes those aspects of geometric positioning 
which are not associated with the image sensors — the external EOS 
geometry. 

4.1 SATELLITE POSITION FROM EPHEMERIS 
4.1.1 Selected Orbit 

The basic requirements for the EOS orbit are nearly the 
same as for the ERTS orbit, and probably will remain so for most sub- 
sequent earth-resource satellites: a near-polar, nearly circular stable 

orbit, with regularly repeating continuous coverage of the earth T s sur- 
face. For EOS and ERTS the orbit is also sun-synchronous , with the 
orbital plane precessing about the earth in inertial space at the same 
rate as the mean 'earth rotation about the sun. Within these constraints , 
the variables are altitude, repeat cycle in days, and right ascension 
of the ascending node. The latter variable determines the local time of 
transit and whether imaging will occur on the ascending (northbound) or 
descending (southbound) node. The various combinations were discussed 
in reference (1), and also in reference (9), and resulted in the EOS 
parameters shown in Table 4-1# The corresponding ERTS values also are 
shown in the table. 


Table 4-1 - Orbit Parameters 


Mean Altitude at Equator 
Repeat Cycle 
Imaging Node 

Mean Local Time at Equator 
on Imaging Node 

Inclination of Orbit 

Distance at Equator between 
Successive Orbits 

Distance at Equator between 
Adjacent Ground Traces 

r 


EOS 

980 km 
17 days 
Descending 

8:30 am 
99+° 

685 km 

171 km 


ERTS 

920 km 
18 days 
Descending 

9:30 am 
99° 

2870 km 
159 km 
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4.1,2 Discussion of Orbit Selection 

In reference (1) , the rationale governing the selection 
of the EOS orbit was discussed in some detail. Basically, the orbit 
was governed by the previous altitude selection of 1,000 km for the 
Thematic Mapper and another lower-resolution sensor. Many prospective 
users of earth-resources satellite image data claim they want higher 
image resolution. They naturally have asked why missions cannot be 
planned for a lower orbital altitude. From the mission-planning view- 
point, low altitudes are avoided for two reasons, increased air drag 
and decreased time above the horizon for transmission to ground receiv- 
ing stations. With increased air drag, it is necessary to periodically 
carry out orbit maintenance with a small thruster. However, as refer- 
ence (1) shows, the additional fuel weight required for this maneuver 
is quite low, at least for orbit maintenance intervals of about 20 days. 
Thus, it is at least feasible to discuss altitudes as low as the 570 km 
treated in reference (1) . 

Decreased above-horizon pass time is a severe limitation 
for satellite programs organized as was ERTS, with a very few stations 
acquiring data directly from the satellite. However, the higher data 
rates to be expected with future satellites, including EOS, raise the 
possibility of using a relay satellite to provide the link between 
imaging satellite and ground stations. Such a solution if adopted 
would have many advantages , only one of which is the removal of con- 
cern about above-horizon time for the imaging satellite. 

So there are ways around the original limitations on 
lower altitudes. Other effects of lowered altitude now must be con- 
sidered, and for these ground-swath width must be considered. 

If ground swath-width in kilometers is kept constant and 
orbital altitude is decreased, the image sensors on the satellite must 
have a wider total angular field of view than at the higher altitudes. 
Since the image sensors are to maintain the same angular resolution, an 
increased data rate is needed to achieve the resulting higher spatial 
resolution over the same ground swath. Also, a larger maximum nadir 
angle will result when imaging from the lower altitude* Terrain relief 
displacements then will be larger. However, as long as the nadir angle 
remains within about 30 degrees very little terrain masking will result 
(conventional aerial cameras have half-field angles of over 40 degrees) • 
Temporal registration between images made by successive passes over the 
same earth scene will not be affected as long as the ground-trace repe- 
tition is controlled properly. 

A second matter of concern is the repeat cycle time. For 
any fixed altitude, the repeat cycle is determined solely by the selec- 
ted ground-swath spacing, reference (1): the narrower the swath width, 

the more days needed for the repeat cycle. In general, a short repeat 
cycle is desirable to give more opportunities for cloud- and haze-free 
imaging of each area of interest. Thus, a wide ground swath is in 
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order. Figure 4-1 shows this relationship for two different orbital 
altitudes, 570 and 900 km. These altitudes are associated with 15 and 
14 orbits per day, respectively, from reference (1). As the figure 
shows, a decrease in orbital altitude allows a choice between a modest 
decrease in either the repeat time or the ground-swath size. For any 
orbital altitude, a significant decrease in repeat time is possible 
only by increasing ground-swath size. Swath size can be much larger 
than the present 160-170 km and still result in a nadir angle of less 
than 30 degrees. The roll offset capability proposed for the EOS The- 
matic Mapper gives an effective shortening of the repeat cycle by in- 
creasing the available swath width on each pass. 

Summarizing the discussion, a low orbital altitude can 
provide more spatial resolution and a slight decrease in repeat cycle 
time. These advantages must be compared against (1) the need for 
periodic orbit maintenance, (2) the possible need for relay satellites 
to provide adequate transmission to ground stations, and (3) larger 
nadir angles and data rates for a given ground-swath spacing. 

Orbits other than circular can be considered for earth- 
resources missions. A highly elliptical orbit would allow imaging 
along a short arc near perigee without long-term air drag. The motion 
of perigee in the orbital plane would shift the imaging arc from one 
pass to the next. Earth rotation and orbital precession would shift 
the ground trace of the orbital plane. Such an orbit obviously is much 
more complex than the simple circular choice. 



Figure 4-1 - Relation of Repeat Cycle and Ground Swath to Attitude 
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A circular orbit does not give constant altitude above 
the earth* As a result of the earth 1 s oblateness, at the highest lati- 
tudes the satellite altitude is about 21 km greater than at the equa- 
tor. The EOS images will have about a 2 percent range in scale as a 
result. This is a small but significant amount and will require addi- 
tional work during digital image processing. The altitude range also 
causes an under- or over-scan condition with a fixed rate scanner, and 
a change to the number of scan swaths per image format. It would be 
desirable to use an orbit with an ellipticity which exactly matched the 
oblateness of the earth, thus giving constant altitude. For a sun- / 
synchronous requirement, this orbit would have motion of perigee in/the 
orbital plane, reference (7), so the earth and orbit ellipticity would 
match only occasionally. Hence, the circular orbit apparently is the 
best way to obtain nearly constant altitude above the earth. 

4.1.3 Effects of Improper Orbit 

For most imaging satellites, insertion into an incorrect 
orbit will significantly affect mission success. It is assumed here 
that departure from the planned orbit is too great for any small adjust- 
ment to be of much help. Then it is likely that sun-synchronism and a 
regular repeat cycle would be lost. The image scale, ground coverage, 
and ground resolution would depart from that designed, as would the 
sidelap between images. 

These changes from design would have severe impact on 
image data processing. For efficient design, a high-production proces- 
sing system requires that a relatively narrow range of parameters be 
assumed. This is true whether analog, digital, or hybrid processing 
schemes are being used. Digital-processing impact perhaps is not so 
obvious, but the special digital logic needed to obtain high production 
will constrain input geometry within narrow limits. The amount of high- 
speed memory provided for image-data shifting also will be only that 
needed for the planned orbit. Finally, any image printer producing a 
film output from a digital system will require special data modification 
to retain a uniform output scale or other pre-designed format, just as 
would the image printer used for any other method of processing. 

The costs are very high to safeguard against improper 
orbit in the design of an image-data processing system. Past launch 
performances indicate that the risk of improper orbit insertion is quite 
small. If this is the case, it does not seem cost-effective to require 
the image-processing system to be designed for such an unlikely 
eventuality. 

Incremental image sensors such as the Thematic Mapper and 
the HRPI collect image data as a function of time. If this time is a 
fixed pre-launch value, a bad orbit will cause along-track data gap or 
overlap. Only if the rotation rate of the Thematic Mapper or the inte- 
gration and readout time of the HRPI can be adjusted by ground command 


4-4 



can the difficulty be corrected. A framing camera, collecting an entire 
image at once, has an advantage in this respect. 

4.1.4 Ephemeris Accuracy 

Some of the image positioning methods discussed in Sec- 
tion 2 require a very accurate value for the position of the EOS satel- 
lite in space at the time each image is being sensed* Since image posi- 
tioning must be done in an earth-based coordinate system, satellite posi- 
tion must also be so defined. The usual positional reference is geodetic 
latitude, longitude, and altitude above the earth ellipsoid of reference. 
This information is provided by the satellite ephemeris, which is deter- 
mined from tracking data. The ephemeris is frequently revised to incor- 
porate the most recent tracking information. A best-fit ephemeris is 
determined for after-the-fact (past) satellite passes, and provides 
highest accuracy. A predictive ephemeris is used to extrapolate ahead 
in time with slightly lower accuracy. 

Position uncertainty is caused primarily by lack of de- 
tailed knowledge about (1) the higher tesseral harmonics of the earth f s 
gravity field, (2) air drag effects, and (3) tracking station location 
with respect to the earth’s center of mass. Air drag effects are well 
in hand, and knowledge of the gravity field’s fine structure is rapidly 
being augmented by reduction of data from many satellites. The center- 
of-mass distance also will be known more closely in the next few years. 
For the ERTS satellite, one sigma position uncertainty was estimated at 
100 meters with the best-fit ephemeris. The along-track component is 
about 80 meters. For the 1976-78 time frame, a one-sigma estimate of 
30-50 meters has been made by Cooley, reference (8) , predicated upon 
improvements in knowledge of the three factors just listed. The lower 
figure is more likely to prevail over the United States, the higher out- 
side tracking-station line of sight. One-sigma tracking station loca- 
tional uncertainties with respect to mass center are expected to be 
considerably less than 10 meters. These estimates have proven to be 
somewhat conservative. ERTS-1 along-track ephemeris errors have been 
reported as 40 meters rms. 

At least one-half the ephemeris uncertainty is estimated, 
reference (7), to be a systematic effect which will repeat for subse- 
quent passes over the same subpoint vicinity. If this estimate is cor- 
rect, repeatability uncertainty for EOS will be 20 to 30 meters, one 
sigma. Error as a function of time is said to take the form of a sine 
wave, with one cycle per orbit, reference (7). If this is the case, 
the positional ephemeris error over an arc of several thousand kilo- 
meters will change by only 5 to 10 meters. 

Other opinions on accuracy differ somewhat from those 
cited above. Stansell, reference (5) , states, ’’Today’s [1971] geodesy 
permits determining a satellite orbit with an accuracy of 10-15 meters 
rms and predicting that orbit ... to an accuracy of 40 meters rms or 
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better." Stansell was writing about the U. S. Navy navigation satel- 
lites, whose typical altitudes and inclinations are similar to that 
planned for EOS. 

Gaposchkin and Lambeck, reference (6) , describe an exten- 
sive analysis which combined dynamic satellite analysis, simultaneous 
satellite photogrammetric triangulation, deep-space network tracking 
data (for relative longitude and distance from tracking stations to the 
earth’s rotation axis), and surface gravity measurements. The combina- 
tion solution provided a check with each of the individual data sources 
and showed generally good agreement. As a result, the geocentric posi- 
tions of 15 stations are determined to 5-10 meters, one sigma, and the 
generalized geoid to 3-4 meters. Orbital residuals of 7-10 meters were 
obtained using a combination of laser and optical data. 

Finally, a completely geometric solution of a global 
station network is nearing completion. This work, reference (11, 12), 
has used photogrammetric satellite triangulation [different observa- 
tions from those described in reference (6)], combined with high-accuracy 
transcontinental baselines. One sigma position uncertainty for 45 
stations is 4.5 meters, reference (10). When these results are combined 
with other data, even better results are expected, which will be of use 
in satellite orbit analysis. 

A lower limit of 5 to 10 meters seems attainable for 
satellite positional uncertainty. The effort to achieve such accuracy 
and the time before it can be routinely obtained are not clear. Thus, 
the remainder of the analysis in this report has used Cooley’s, refer- 
ence (8), lower one-sigma repeatability error of 20 meters. 

4.2 ATTITUDE DETERMINATION REQUIREMENTS 

This section provides some additional details on the spacecraft 
attitude control system, inertial-attitude reference packages, and star 
sensors. The need for obtaining a more complete knowledge of EOS space- 
craft attitude with time is stressed. Finally, in its admitted defi- 
ciencies and erroneous interpretations, the need again is emphasized for 
more communication and closer contact between the two areas of technology 
represented by the satellite-data and control-point methods of posi- 
tioning images. 

The section is divided into two parts, the actual attitude be- 
havior of the spacecraft and the measurement of this attitude by 
auxiliary sensors • 

4.2.1 Image Continuity Requirement 

The EOS satellite must maintain a local-vertical or down- 
ward-looking attitude, at least during imaging. Yaw must be constrained 
to allow the Thematic Mapper and HRPI to scan normal to the orbital 
path. Any attitude corrections must be made slowly between the allowed 
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limits • These conditions are maintained by the attitude control sys- 
tem (ACS) . 


Attitude Error Amplitude 

The baseline ACS is described in reference (1) . Each of 
the three axes is allowed a bias of +0.5 degree maximum, plus 0.2 degree 
rms about the bias. This implies maximum errors of 1.1 degree in any 
axis, with rms values of 0.7 degree. These limits appear to be ade- 
quate, No potential users are known who require closer tolerances than 
these for data collection. Even larger pitch and roll errors may be 
permissible. If lower attitude errors would prove of benefit, the 
greatest improvement with the least effort can be made by measuring and 
removing the 0.5— degree attitude biases. The high-accuracy attitude 
measurement technique to be used for image positioning (either stellar 
sensor or GCP) provides this information. By suitable ACS design, the 
measured biases could then be removed from the ACS using ground commands. 

There are some opposing arguments concerning the allow- 
able yaw error. On the one hand, the image displacements caused by yaw 
are at most only one-tenth those of either pitch or roll for the same 
angles, so it would seem less important to control yaw as tightly as the 
other two angles. On the other hand, a digital geometric correction of 
image data for yaw is troublesome for high-speed processing, since it 
entails more than the simple data translation which is needed for pitch 
or roll adjustment. It would be desirable to eliminate this adjustment, 
or at least reduce it as much as possible. Still, extensive geometric 
processing will be needed even with zero yaw correction, to adjust for 
scale change, map projection, internal sensor geometry, and earth-rota- 
tion effects. Finally, a conical-scan Thematic Mappdr is somewhat sen- 
sitive to yaw angle in terms of continuous scan coverage. The optimum 
yaw angle differs with latitude (Appendix A) ; this adjustment could 
possibly be incorporated into the ACS. Adjustment precision of a few 
tenths of a degree would be adequate. To enforce such an adjustment 
might mean a tighter control for yaw than for pitch or roll. 

In summary, the baseline yaw variation of +0.2 degree rms 
is not obviously too large or too small for the imaging and processing 
requirements. In the absence of more critical considerations, it is 
recommended that the baseline control limits be retained. 

A horizon sensor will collect the basic information used 
to control attitude. This sensor may be more precise than that used 
for the ERTS ACS. Telemetered horizon-sensor and torquer data will be 
used to develop approximate attitude values from ground processing. The 
after-the-fact attitude determined in this way is expected to be accur- 
ate to 0.1 degree rms in each axis. Some other high-accuracy attitude 
measurement obviously will be needed to suitably accurate image 
positioning. 
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It is possible that higher attitude accuracy may be obtain- 
able from the telemetered AGS data by removing the earth-ellipticity 
horizon effects. These effects are believed to be responsible for much 
of the attitude error. One-sigma errors of 0.005 to 0.02 degree are 
thought to be feasible (20), using advanced horizon sensors and exten- 
sive ground data processing. More definitive information is needed, but 
there appears to be a potential for using ACS data as backup for a 
higher- accuracy independent attitude measurement system such as SIMS. 

For GCP positioning, such a capability would be of some limited useful- 
ness in establishing initial absolute image locations. This would be of 
benefit when analyzing images from those areas of the world without suit- 
able maps or other data from which to select GCP. 

Attitude Error Rate and Acceleration 

The baseline ACS specification defines dynamic attitude 
behavior with a single parameter of 0.05 degree per second maximum rate 
of change in any axis. This figure was chosen primarily to ensure that 
minimal overlap or underlap (less than one-fourth 1F0V) between scan 
swaths would occur with the Thematic Mapper. There may be ACS designs 
with which improved positional accuracy for images could be attained by 
increasing the maximum allowable rate. If so, it is important to know 
exactly what rates are associated with what amount of over- or underlap. 
This relationship is summarized in Figure 4-2. The EOS altitude of 980 
km has been assumed. Then the allowable rate depends only upon (1) the 
number of scan lines in a single scan swath, and (2) the fraction of a 
single IFOV allowed for overlap or underlap. As the figure shows, the 
present maximum rate is well within the limitations of six scan lines 
per swath and one-fourth IFOV lap. 

ACS accelerations are not specified for EOS, and are of 
no concern for image continuity. However, in reference (23) a maximum 
value of 0.0114 degree/s ec^ is assumed. Acceleration durations are 
very short since the reaction wheels are acted upon by impulsive torques. 
Thus, the acceleration behavior over several seconds still is not 
specified. 

The horizon-sensor/reaction-wheel ACS is basically of the 
type used for ERTS, with two main exceptions. The maximum attitude rate 
has been reduced by a factor of ten from ERTS. Three-axis gyros will be 
used for EOS, unlike the single gyro used with ERTS for roll/yaw con- 
version. Maximum accelerations are specified based on the maximum cor- 
recting torque applied to the inertia reaction wheels. When one of the 
wheels has reached maximum angular velocity, it is "unloaded 11 with a 
cold-gas impulse to nearly zero angular velocity. Based on experience 
gained from previous operational satellites, the frequency of wheel un- 
loading will be very low, perhaps once every few days. During the brief 
time required for wheel unloading, the attitude rates are much greater 
than the operational maximum given above. Any EOS image data obtained 
during this time will not be usable, and it will be very difficult to 
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Figure 4-2 - Maximum Attitude Rate for Imager Considerations 

accurately position the individual image frame collected over the 30 
seconds during which unloading occurs. As a consequence, it is quite 
desirable that unloading not occur during imaging . 

4.2.2 Image-Positioning Attitude Requirements 

The influences of the ACS on image continuity are widely 
recognized and, as seen from Section 4.2,1, have been adequately addres- 
sed by the EOS specifications (1). Unfortunately, the influences of 
the ACS on image positioning methods apparently have not been recognized; 
this subject is ignored in the specifications. ACS design directly 
delimits the positioning methods which can be used for EOS images . This 
is an inescapable consequence of any orbital mission which uses any 
imagers other than high-speed frame cameras. It is particularly import- 
ant when the image is acquired rather slowly, as is the case with the 
EOS line scanners which build up the image format over tens of seconds. 

No single set of ACS requirements can be given for image 
positioning, unless only a single positioning method is being considered. 
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Such is not the case here. As a consequence, the costs, development 
time, and success probability of the ACS design become an integral part 
of the decision to be made on image positioning method. Several 
examples will quickly illustrate this fact. 

Consider first an ACS design which makes no allowance for 
image positioning method, using only the requirements for image contin- 
uity developed in Section 4.2.1. Any of several designs will do for 
this case. Without any further conditions to be met, image positioning 
immediately requires addition of an independent gyro package for inde- 
pendent attitude angle measurement, and telemetry of the attitude values 
frequently enough to supply attitude changes of 0.001 degree in pitch and 
roll. (Yaw can be relaxed, if necessary, to several times this amount 
of change.) For the present study, such an ACS design would immediately 
eliminate method GCP-E from consideration. The only methods which 
could be expected to satisfy such conditions are SIMS-E, SIMS-EG, and 
GIMS-E-1. 

Next, consider the opposite extreme, a high-accuracy ACS. 
This possibility has been discussed in Section 7.7.1 of reference (1). 
Essentially, this control system would be operated by closing the loop 
between a SIMS package and the control reaction wheels. There are sev- 
eral factors which argue against this scheme.* Disregarding these for 
the moment, such an ACS concept would leave very little to be done by 
either an auxiliary-sensor method or a GCP method. With attitude errors 
always 0.001 degree or less per axis, independent attitude measurement 
and telemetry would be needed on the spacecraft. An initial bias removal 
of attitude-transfer errors after launch would be needed, repeated oc- 
casionally during the operational life of the satellite. One image with 
only two GCP would theoretically suffice for this task (in practice, 
several more GCP probably would be used). The largest positioning error 
then would be caused by the ephemeris data. The systematic component of 
this error can be removed by positioning a single image with GCP on each 
pass. These minimal techniques do not correspond exactly to any of the 
5 methods considered in this study, although the instrumentation includes 
some of the same components. 


The most significant probably is that users of EOS images do not need 
such absolute angular accuracy during data collection. Another reason 
is the inadequate accuracy, as given in reference (1): 175 meters is 

quoted, which is completely inadequate for temporal registration of 
images with 20— meter and 44— meter resolution elements. The onboard 
computer required to convert inertial to local-vertical attitudes is 
another difficulty with the concept in reliability, weight, and power 
requirements. The ephemeris developed by or transmitted to this com- 
puter from the ground would of necessity be a predictive ephemeris for 
real-time attitude control, with consequent lower accuracy than the 
best-fit ephemeris gives for after-the-fact determinations. 
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Between these extremes are a number of opportunities for 
ACS design which, considered along with the possible image-positioning 
methods, may well include the most cost-effective route for EOS. For 
all of these intermediates, the absolute attitude error amplitudes and 
rates are of no concern, in themselves. The behavior of interest is given 
by the time interval over which a straight line or second-order curve 
can be best fit in a least-squares sense to the attitude errors in each 
axis, with an rms error in the range 0.0005 to 0.001 degree per axis. 

The time intervals of interest range from 30 seconds to several hours. 
Obviously, the character of the error is important in this specification. 
Random errors of relatively high frequency can be removed with an inde- 
pendent gyro package used to measure short-term attitude errors. In 
this way, the achievable best-fit time intervals may be appreciably 
lengthened. However, the gyro package may be incorporated directly in- 
to the control loop of the ACS, in which case it no longer is an 
optional positioning-method component. As a consequence, it is recom- 
mended that the inclusion or non-inclusion of a gyro package, and 
whether it is a part of or independent from the ACS control loop, be 
left to the ACS design group. For the different image positioning 
methods, the time intervals and line-fittings of interest are as follows: 

(1) 30 seconds, second-order (without gyro package) 

(2) 30 seconds, first-order (with and without gyro package) 

(3) 5 minutes, first-order (with and without gyro package) 

(4) 5 minutes, second-order (with and without gyro package) 

(5) 10 minutes, first-order (with and without gyro package) 

(6) 100 minutes, first-order or second-order (with or without 
gyro package) 

(7) Several hours, first-order or second-order (with or without 
gyro package) 

There is probably little difficulty in an attitude system design which, 
with an independent gyro package, could meet the specified rms error 
requirements over 10 minutes. However, it is important in selecting 
positioning methods to know whether the same requirements can be met 
without the package, and what the costs would be. The intervals given 
are associated with angle updates once for each image, once for each 
landmass, and once for each orbital pass. 

In summary, the image-positioning attitude requirements 
for EOS differ from the image-continuity requirements in being concerned 
with smoothness of angular changes , rather than rates or amplitudes of 
attitude angles. No single set of requirements can be made for the 
ACS; requirements should properly be based upon system considerations 
which include image positioning methods as well. 
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4,2.3 Accurate Attitude Measurement by Auxiliary Sensors 


Present plans call for accurate measurement (0.001 to 
0.003 degree) of EOS attitude to be performed by a technique which is 
independent of the attitude control system, reference (1) . Continuous 
and highly accurate knowledge of attitude behavior is needed to properly 
position the line-scanning HRPI and Thematic Mapper image daga. It is 
possible that the ACS design could be such that quite smooth attitude 
behavior could be assumed. In this case, as mentioned previously, ac- 
curate measurement need be made only infrequently. In the following dis- 
cussion, however, it is assumed that continuous high- accurate attitude 
measurement is needed. Gyros on the satellite would provide this data. 
The information is not truly continuous, but the sampling and telemetry 
interval can be made adequately small compared to the spacecraft attitude 
response. For EOS, the 0.1-second interval presently being considered 
is more than adequate. 

Gyros must be updated occasionally to correct for drift. 

The frequency at which this update is required depends on the specific 
gyro design and the configuration. References (23) , (24) , and (25) 
give extensive consideration to gyro selection and configuration. This 
information would have been of greater value in the present study if (1) 
it had been available earlier in the program, and (2) the work described 
had not been optimized specifically for the SIMS technique. The gyros 
finally selected in reference (25) have performance matched to the 
update-frequency capability of the star-sensor. 

Stellar-sensor and GCP methods are both candidates for use 
in attitude update. The following comments apply to the star-sensor 
method. GCP considerations are contained in Section 4.3. 

Transformation of inertial attitudes to a local-vertical/ 
velocity-vector reference requires detailed knowledge of the earth’s 
orientation in inertial space and the location of the satellite with 
respect to the earth ellipsoid of reference. Small errors in the former 
reference frame result from polar motion and other effects. The errors 
in the latter stem from ephemeris errors, discussed in the previous sec- 
tion, and in accurately measuring the time of the stellar observations. 
Ephemeris errors will likely be the most significant. Their effect in 
positioning will be to replace the translational ephemeris errors by 
equivalent attitude errors. Additional care and complexity are required 
for the stellar-ref erence computations in defining the final image 
location in the geodetic coordinates used for geographic point position- 
ing, instead of the geocentric or reduced coordinates normally used with 
satellite operations. The GCP method requires no such care on this 
account; the observations used for positioning are supplied in geodetic 
coordinates. 


Some additional cautions in transforming from inertial to 
geodetic coordanates have been voiced by Schmid, reference (31) , in con- 
nection with a photo gramme trie triangulation using stellar control. 
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He found bias errors in his work which he ascribed to (1) lack of a 
rigorous connection between elevation and the reference ellipsoid which 
serves as the computational system for latitude-longitude determination, 
and (2) the fact that the latitude-longitude and right ascension-declina- 
tion systems are related only insofar as the primary axes of the two 
systems are parallel to each other. Schmid removed the bias by allowing 
an additional spatial rotation in his baseline orientation* 

Some estimates of stellar sensor costs, weights, and power 
requirements are given in reference (25). Power needs are very low; 
weights range from 3 to 16 kg depending on whether the sensor is a star 
mapper or star tracker. Page 6-29 of reference (25) comments on reli- 
ability: "None of the SIMS candidates has reached the development status 

that will support a comprehensive, quantitative assessment of reli- 
ability * 11 

The need for an accurate attitude-transfer device in con- 
nection with an auxiliary attitude sensor has been commented on in Sec- 
tion 3.4.1. Such a measurement Is needed to relate attitude data ob- 
tained by sensors which are independent of the image sensor. The reli- 
ability, power, weight, and telemetry requirements for such a device 
(or devices) probably could be designed to be quite satisfactory, but 
no detailed thought appears to have been given to such design. 

4.3 GROUND CONTROL POINTS 

4.3.1 Backg round 

One basic method of positioning images with respect to the 
earth's surface is through the use of ground control points (GCP) . The 
terns landmark points or geometric ground-truth points also are sometimes 
used for GCP. Such points are selected well-defined image details for 
which the earth positions are known. Positioning an image with GCP 
typically requires that several GCP be identified and measured on the 
image. The measured image locations and known ground positions of these 
points determine the geometric relationship of the entire image with 
respect to the earth's surface. For the procedure to be effective with- 
out requiring excessive GCP, internal sensor geometry must be well- 
defined. 


GCP have been used for years in photogrammetric mapping 
using aerial photographs, where the method is almost universally applied. 
(Photogrammetric methods are used for virtually all original map com- 
pilation performed today.) The GCP used for this purpose typically must 
be specially surveyed to achieve the high accuracies required. 

The mathematical GCP positioning process for a single 
photograph is called spatial resection , analogous to surveying resection 
on the earth's surface. The unique position and attitude of the aerial 
camera is determined as six orientation elements : three positional co- 

ordinates and three rotations. No auxiliary data from any other sensors 
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in the aircraft normally are used to determine these quantities, only 
the GCP. 


Two separate equations can be written for each control 
point, one for the measured x and one for the £ image coordinate of the 
point. The six unknown orientation elements appear in the same equa- 
tions, along with the three positional coordinates of the control-point 
on the earth. Thus, three GCP with their six equations are adequate to 
solve for the six unknown orientation elements. Best accuracy is ob- 
tained if the points are well distributed throughout the photo; the 
solution is indeterminate if the GCP lie on a straight line. With the 
six orientation elements known, the ground-plane positions of all other 
photo images also are known. (If height differences in the photographed 
scene are a significant fraction of camera altitude, a two-photo tech- 
nique is used, in which the intersection of conjugate image rays uniquely 
defines the ground position in all three coordinates. For satellite 
images from narrow-angle, downward- looking sensors, relief effects are 
small and potentially repeatable.) 

In adapting the use of GCP from aerial photographs to the 
ERTS satellite images, some differences were noted. First, GCP accuracy 
for image positioning need not be nearly so high as that needed for con- 
ventional mapping. As a result, it usually is possible to select GCP 
from existing maps and other available sources without special field 
survey work. A second difference is in the large number of GCP poten- 
tially available for use, with the result that extra GCP can be used in 
the spatial resection* There are several advantages in using redundant 
data; these will be discussed later. A third difference is in the use 
of electronic image correlation to match GCP against a master image. 

Other implementation changes from usual mapping practice were incorpor- 
ated. Finally, some form of auxiliary attitude-change information was 
recognized as desirable for positioning the ERTS MSS (scanner) images 
acquired. All of these differences apply to EOS as well. 

The great advantage of ground control points for position- 
ing is the elimination or great reduction in the reliance that must be 
placed on externally derived location and attitude data. In addition 
to the resulting economy and reliability increases that follow, the use 
of GCP usually results in higher positioning accuracy than can be at- 
tained using auxiliary data. As image resolution increases, the posi- 
tioning accuracy typically increases accordingly. 

Use of GCP has disadvantages. The major one is in imple- 
mentation in what otherwise may be a more or less completely automatic 
processing scheme. For high-production programs with many images to be 
processed, implementation must be given careful consideration. Equip- 
ment must be selected to assist a trained specialist in GCP selection 
and accurate measurement, in as efficient a way as possible. However, 
although equipment and procedures are important, proper personnel are 
vital. Control-point selection requires a high degree of self-motivation 
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and job satisfaction, akin to mapping and some kinds of photo interpre- 
tation. It is not difficult to find people with these characteristics, 
provided the right kind of people are sought. 

The main difficulty with GCP is in identification. In 
aerial photogrammetric practice, misidentification is the greatest prob- 
lem associated with GCP. Work with space images cannot be expected to 
be immune from the same difficulty. The use of well-trained responsible 
personnel is the most effective way of minimizing the problem, but tech- 
nology can be of considerably assistance as well; once the initial GCP 
selection has been made, automatic image correlation methods can be 
used to perform subsequent matching operations. Using extra GCP also 
helps greatly in detecting misidentif ied GCP. Initial results with 
ERTS— 1 images indicate that these techniques work very well in reducing 
the incidence and effects of misidentif ied GCP. 

4.3.2 Accur acy 

The accuracy of positioning with spatial resection depends 


upon five 

main factors: 

(i) 

Ground coordinate accuracy 

(2) 

Pointing precision 

(3) 

Image adjacency effects 

(4) 

Internal sensor geometry 

(5) 

Effective redundancy 


The equipment errors associated with measuring image-detail location and 
producing the image to be measured are not treated here, since they con- 
cern implementation, not basic limitations. The other five factors are 
discussed in the follox^ing paragraphs. 

Ground Coordinate Accuracy 

Absolute and relative accuracy are important distinctions 
here. Absolute positional accuracy of the coordinates is defined with 
respect to the earth ellipsoid of reference. Elevation normally is 
taken as distance above sea level, although this differs from the refer- 
ence ellipsoid by the separation between geoid and ellipsoid. Relative 
accuracy means that introduction of some bias into the absolute coordi- 
nates is permissible; repeatability is the important thing. 

U. S. National Map Accuracy Standards require that the 
absolute positional accuracy of well-defined map features be within cer- 
tain limits. These limits are equivalent to 0.3 mm rms at map scale. 

GCP for EOS should be selected from the standard U. S. map scales of 
1/24,000, 1/50,000 (Alaska only), and 1/62,500. Additional high-accuracy 
positional information is available on highway structures and airfield 
runways. Some of the regular survey points in the second- and third-order 
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geodetic networks of the country lie at or are sufficiently close to 
features which can be identified on satellite images. In general, con- 
trol-point accuracy need not be a problem, at least for positioning 
images in the U. S. or other nations with adequate surveying and carto- 
graphic activity. If smaller-scale maps are used for control-point 
selection, inaccuracies in the plotted map positions may significantly 
affect the final absolute position accuracy. But for many applications 
absolute position accuracy is less important than the repeatability 
needed for temporal registration. If the same GCP are used to position 
subsequent images of the same earth scene, the ground position errors 
will affect all images alike and good temporal registration will be 
achieved. Of course, it may not always be possible to use exactly the 
same GCP every time, because of cloud cover and seasonal changes. In 
this case, only that part of the ground coordinate error which is com- 
mon to all GCP in the area can be ignored. One way of overcoming this 
problem and ensuring repeatability is some form of control intensifica- 
tion . The best available GCP are used to position the first image ob- 
tained of an earth scene. Auxiliary GCP then can he selected from this 
image to use in positioning subsequent images. Absolute accuracy will 
be no better than — in fact, not quite as good as — that of the first 
image, but temporal registration will be excellent. 

Pointing Precision 

When measuring the image location of a GCP, some form of 
measuring mark or cursor must be placed directly over the image detail 
which is to be used as the control point. This mark placement is called 
pointing . 

Two kinds of pointing can be considered. The first is 
called descriptive pointing here, and occurs when a GCP is originally 
selected. Usually some written or graphic description of the GCP refer- 
ences the pointing operation. Pointing is done manually in this case, 
using optimum image illumination and magnification, together with the 
measuring mark of best shape and size. If the same GCP is used again 
on later images of the same earth scene, descriptive pointing may be 
repeated in the same way used initially. If possible, however, it is 
more efficient, reliable, and usually more accurate to use the second 
kind of pointing, called here image matching or correlative pointing . 

This is done by binocular ly comparing two images. One is the new image 
in which the GCP is to be pointed; the other is a master image on which 
the GCP has previously been precisely located using descriptive pointing. 
Special optical devices may be needed for the image matching process to 
allow one or both of the images being viewed to be rotated and scaled so 
as to permit binocular superimposition. Image matching can be done auto- 
matically as well, using digital or analog image correlation techniques 
and appropriate equipment. The accuracies and times required to carry 
out pointing with such methods vary widely, depending on many factors. 
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The ability to precisely point at the correct image de- 
tail is affected by image quality, the type of mark used, and the ability 
to make fine adjustments to the mark. The latter two matters are con- 
siderations in equipment design and need not concern us here. Image 
quality is important* Several investigators, references (15, 14), have 
tried to quantize what is intuitively obvious: the higher the image 

quality the more precisely can one point to a selected image detail. 

In general, it is possible to center a measuring mark on a symmetric 
target within one to five percent of the target size. This ability de- 
pends upon the type, size, and contrast of the measuring mark and target, 
as well as such other factors as illumination and granularity in the 
displayed image. 

From a limited amount of experience with ERTS images and 
photos from the Apollo 9 S065 experiment, it appears that when GCP are 
used which are asymmetric and poorly defined in the image, pointing 
error may be the principal factor in overall positioning error. As 
examples of this hypothesis the center of a nearly circular lake is bet- 
ter as a GCP than a peninsula extending into the lake; the center of a 
small regularly shaped island in the lake is still better. Similarly, 
the X intersection of two roads is better for pointing than a bend or 
corner in a road pattern. A Y intersection of two streams is a better 
GCP than a loop in a stream. Additional work is needed to confirm or 
refute this hypothesis, as well as on several other aspects of pointing 
to satellite images displayed in raster form. 

The precision of descriptive pointing does not appear to 
have been explored for raster images. Thus, it is difficult to develop 
the errors associated with this process. The matter is further compli- 
cated by the quality of the ERTS-1 images, which has exceeded all pre- 
launch estimates. An extremely limited amount of data obtained from the 
ERTS MSS images indicates a standard deviation in pointing precision of 
about 20 m, or about one-fourth the IFOV. The Thematic Mapper then 
should give a precision of 11 m, and the HRPI 5m. In light of the un- 
certainty in this value, a range of 10-20 m will be assumed here for the 
EOS precision. 


Image Adjacency Effects 

Some images are created using processes and materials which 
alter the position of a light-dark edge; photographic film is an example 
of a medium which exhibits these undesirable adjacency effects , refer- 
ence (16) . When it is necessary to measure GCP from such a displayed 
image, control points should be selected with special preference for 
symmetrical details, or at least details with low density difference 
from their background, in order to minimize the metric effects of edge 
migrations. 
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Internal Sensor Geometry 


Spatial resection uses some mathematical model to repre- 
sent the internal geometry of the image sensor. To a greater or lesser 
extent, this model errs in its representation, and so introduces image- 
location errors. For GCP, these errors adversely affect resection 
accuracy, resulting in positional inaccuracies for the image as a whole. 
Any image detail not used as a GCP also has an additional error associ- 
ated with the internal sensor geometry error at that particular point. 

For line scanners, the motion of the scanner platform 
provides the second dimension of the two-dimensional scanner image. To 
this extent, the description of the platform position and attitude as 
functions of time can be considered as vital parts of the scanner- image 
internal geometry. For a scanner on an orbiting satellite the platform- 
position changes can be modeled quite accurately as functions of time 
over arcs of hundreds of kilometers. For attitude behavior of the satel- 
lite, three courses are open: control the attitude precisely and con- 

tinuously with respect to the local vertical so that modeling is extremely 
simple and needs no spacecraft data; design the spacecraft attitude con- 
trol system such that attitude changes are suitably slow and regular as 
functions of time; or use telemetered changes in attitude over the 
imaging arc as auxiliary data in the spatial resection. These alterna- 
tives were discussed in Section 4.2. For any of these methods, the 
limitations on scanner internal geometry are those placed by the magni- 
tudes of high-frequency attitude changes, for which no information is 
obtainable. The second alternative is particularly attractive. If the 
spacecraft control response to a disturbing torque is such that attitudes 
can change only slowly over the imaging arc, first- or second-degree 
attitude changes can be carried as additional spatial resection unknowns. 

A few more GCP then provide adequate observational data to determine 
these unknowns. This attitude control method results in the same virtue 
of the first method above: no need for auxiliary spacecraft data. 

Effective Redundancy 

If more than the minimum number of GCP are used in a 
spatial resection, a least-squares solution is used to determine the 
unknown orientation elements. By reducing the random effects of the four 
error causes discussed above, positioning accuracy of the image as a 
whole is statistically improved in accordance with the number of GCP 
and their locations within the image. The use of redundant control 
points has several operational advantages in addition to overall ac- 
curacy improvement. First, self-checking: the least-squares solution 

gives several accuracy indicators. These result from the residuals. 

(A residual is the distance of each measured GCP location from the best- 
fitting math-model location, as determined by the least-squares solution.) 
The residual for a single control point provides an accuracy indicator 
for that point, particularly useful in detecting misidentif ied GCP. 

Also, standard deviation of all residuals can be used as an overall 
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positioning quality indicator for each individual image. This is use- 
ful as a single-number value index. Residuals also are helpful in long- 
range analyses of operator ability, equipment performance, and sensor 
calibration, to name some additional applications. 

The GCP*s used for satellite images are not specially sur- 
veyed points, but selected map points. Thus the cost per point is very 
low, making redundancy more attractive than is the case for normal 
aerial photomapping. For maximum benefit, the GCP should be selected 
in particular image areas. In general, the best image areas are those 
which (a) allow uniform GCP distribution within the image, and (b) are 
located in image areas of greatest sensitivity to orientation-element 
errors. In general, greatest sensitivity is found at the corners and 
edges of the images, where yaw and altitude changes have their maximum 
effects. If the GCP placement does not cover a large portion of the 
image, positioning accuracy will be lost. 

A mathematical analysis can be performed to show the re- 
dundancy advantage for any number and distribution of GCP. Such an 
analysis was done for the KBV (17) and MSS (18) ERTS images, based upon 
earlier work carried out during the ERTS Phase B/C study, reference (19) . 
Figure 4-3, taken from reference (19) , shows the decrease in relative 
positioning error associated with the placement of GCP in the image and 
the number of GCP used. 


TOTAL IMAGE FORMAT IS 185 KM 



Figure 4-3 - Example of Decrease in Positioning Error with Increasing 

Image Coverage and Numbers of GCP, Taken from Reference (19) 
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Additional analysis along these lines was carried out 
during the present study. Of interest was the effect of an increase in 
the number of unknown attitude parameters on the accuracy of the spatial 
resection. The method is described in reference (19). Several differ- 
ent GCP configurations were considered in the analysis, which neglected 
some minor effects. The results are summarized in Table 4-2. Maximum 
array size used for the GCP was 146 km square (Case 2 in Figure 4-3). 

The table shows that a requirement to determine mean attitude angles, 
rates, and accelerations with GCP need not degrade results significantly 
from a determination of angles and rates alone. The number in the body 
of the table represents a multiplier to be applied to the standard error 
of unit weight for a single GCP image point. The multiplied value be- 
comes the rms positioning error as a result of the spatial resection. 

The values range from about 0.35 for an angle-only solution, (three 
unknowns) to about 0.60 when angles, mean rates, and mean accelerations 
(nine unknowns) must be found. The results of the analysis accrue to 
the advantage of method GCP-E, and place in sharper focus the need for 
an ACS with predictably smooth behavior characteristics. 

4.3.3 Applications 

The GCP positioning method can be adapted for many differ- 
ent applications. Some of the applications discussed in this section 
are common in aerial photogranmietric practice. Others are of interest 
primarily to multi-sensor satellite programs. 

Perhaps the most obvious application of the GCP method 
might appear to be as an independent check on attitude (determined from 
other spacecraft sensors) and position (taken from the spacecraft ephe- 
meris, itself determined by tracking the satellite). This is not neces- 
sarily a valid application. In general, the wider the image-sensor 
field of view, the better determined are the six orientation elements. 
For downward-looking narrow-angle images, such as the ERTS RBV and MSS, 
it is not possible to isolate the effects of pitch from along-track 
translation. The effects of a change in either orientation element on 
the image are virtually identical. The same is true of roll and across- 
track translation. Such a lack of determinacy does not affect the 
accuracy with which images can be positioned with respect to the earth, 
only the accuracy in isolating the two parameter pairs which are highly 
correlated. 


The correlation between orientation elements can be used 
to some advantage for EOS images, as is being done for the ERTS program. 
The satellite ephemeris gives along-track and across-track position to 
within 100 m rms for the ERTS-1 satellite. By considering the ephemeris 
along-track and across-track values to be errorless, the spatial re- 
section computation can use two fewer unknowns. In this way the degrees 
of freedom are increased for a given GCP redundancy, resulting in a 
formal improvement in the statistical positioning quality. The small 
errors that are actually present in the ephemeris positions are 
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Table 4-2 - Relative Position Error After Thematic-Mapper Spatial 
Resection, Using Different Error Models 


Number of Unknowns GCP Configuration 



3, Angles only 

0.41 

0.42 

0.41 

0.38 

0,36 

0,38 

6, Angles and Mean Rates 

0.53 

0.59 

0,58 

0.49 

0,50 

0,61 

9, Angles, Mean Rates, and 

* 

0.60 

0.58 

* 

0.51 

0.61 


Mean Accelerations 


(1,00 = Standard Error of Unit Weight for a Single GCP Observation) 
* 

Insufficient data using this configuration. 
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compensated almost perfectly in the solution by slight increments in 
the computed pitch and roll angles. The maximum error, dS, introduced 
by this approximation occurs in the image corners of a framing sensor or 
along the edges of a scanner. For a downward-looking framing imager, 
dS is given by the usual photogrammetric differential, references (21, 
22 ) : 


2 2 

dS = f- d9 - ~ dE 
H H 


where 

S - maximum ground-distance from the nadir which appears in the 
image, ~ 130 km for EOS 

H - spacecraft altitude, - 1,000 km for EOS 

d0 = angle increment to correct for ephemeris error 

dE = position error from ephemeris, - 20 m rms for EOS 

so that dS is much less than one meter rms. Even a 200 m ephemeris error 
would be acceptable, giving a maximum error of 3 meters* 

As the image sensor is tilted farther away from the nadir, 
the correlation between the angular and translational orientation 
elements decreases in the direction of tilt. Even for an image rolled 
30 degrees off the vertical, the ERTS resection technique — carrying 
ephemeris position fixed in the solution — can be used for offset- 
pointed EOS sensors as well. 

The high correlation just discussed between angular and 
positional errors is an example of an important characteristic of GCP 
positioning: the composite nature of the results. Any error component 

which has the same image displacement effect as an angular or positional 
orientation element will be removed as a part of that orientation ele- 
ment. For example, the yaw of the image sensor may be determined with 
high accuracy, but without further information it is not possible to 
separate the component caused by spacecraft yaw from that caused by yaw 
misalignment between the image sensor and the spacecraft yaw axis. Simi- 
larly, changes in spacecraft altitude, camera focal length, film dimen- 
sions, or voltage in the image-sensor power supply all may have the same 
net effect in the image: a simple scale change. This compositing effect 

has both advantages and disadvantages. On the one hand, it automatically 
reduces the number of unknowns to only those necessary and sufficient 
to position the image on the earth 1 s surface. On the other, it requires 
some care when attempting to use GCP to isolate a particular error com- 
ponent. Nevertheless GCP are very useful in several special applications 
of concern for EOS operations. Four of these are discussed here. 
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Detailed Sensor Geometry 

Assume that more than the minimum number of GCP are used 
in a spatial resection. The image-coordinate residuals after the re- 
section represent geometric errors which could not be compensated by the 
spatial resection unknowns. Part of each residual is random, associated 
with errors in pointing, ground-coordinate measurement, and some kinds 
of internal sensor geometry. The other component is systematic, typi- 
cally resulting from uncompensated sensor behavior. It is this syste- 
matic component which shows the extent to which the mathematical model 
being used for the imaging geometry is deficient. If the random com- 
ponents can be detected and removed from the residuals, the systematic 
components remaining will show how the imaging geometry model should 
be modified. Incorporating these changes will improve the accuracy of 
all subsequent spatial resections for that image sensor, providing the 
systematic errors do not change. 

The random effects of the GCP residuals are removed by 
making an error anlysis which sums the effects of a great many residuals. 
The random components (with the normal-distribution zero-mean conditions 
which usually prevail) tend to cancel, and the systematic components 
reinforce one another. Residuals may be accumulated from only one or a 
few resections, each with many redundant GCP, or from the results of 
many resections, each with only slight redundancy. 

Examples of typical image-geometry errors detected in this 
way include lens distortion effects in cameras and sweep-angle nonline- 
arity in scanners. Sometimes the systematic errors are in the image 
processing equipment, not the original sensor, and the corrections may 
take the form of polynomial functions of the x and image coordinates. 
Analysis of residuals can be done graphically or computationally. If 
computational, some a priori model usually is assumed for the error form, 
although this may be revised in the course of the analysis. : 

Analysis of GCP residuals is a very powerful technique, 
and can be expected to see much use in satellite image analysis. This 
application is particularly useful for checking on slow changes in 
sensor geometry, such as can be expected from the MSS scanner on the 
ERTS satellites. 

Alignment of Two Image Sensors 

This application has already been used operationally with 
the ERTS-1 RBV three-camera array to precisely measure the angular 
alignment of the cameras. Such alignment accuracy was needed to produce 
event- registered multi-spectral color composite images. For this appli- 
cation, true GCP are not needed, only conjugate image details which are 
sharply defined on both images being aligned. The measured locations 
from one image serve as "ground coordinates," and a modified spatial 
resection provides the alignment of the second image sensor with respect 
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to the first* For the ERTS application, the "red 11 RBV camera was the 
reference to which the other two were aligned* 

This application of GCP positioning can give very accurate 
results, since a great many common points — 20 or more — can be easily 
located and used in the resection, and very small and sharp micro-details 
can be selected which do not appear on any map. 

Alignment Between Image Sensor and 

Auxiliary Attitude Sensor 

This application of GCP is useful for several different 
programs. A few examples are: 

• Determination of initial alignment angles after launch 

• Measurement of changing alignment angles caused by spacecraft 
flexure or attitude-sensor drift 

• Performance check on autocollimator or other device used to 
enforce the attitude-transfer requirement 

• Removal of initial bias component in gyros used as inertial 
measurement devices. 

For this kind of application the positional ephemeris error acts as an 
error in the procedure. This is once again a consequence of the high 
correlation between the pitch/ along- t rack and roll/across-track orienta- 
tion elements. The angular error, d0, as a result of ephemeris error, 
dE, is simply d0 - dE/H, where H is spacecraft altitude. With position 
fixed at ephemeris values, the error from a single spatial resection 
cannot be less than d0. For the EOS ephemeris, d0 will be about 
20 m/ 1,000 km rms, or about 0.001 degree. Interestingly, the spatial 
resection does not give the true local-vertical attitude with this 
precision, but an attitude which is locally repeatable from the bias 
errors in the ephemeris. The effective error in measuring angular 
alignment can be reduced further by performing several spatial resec- 
tions, selecting the locations in such a way that the ephemeris errors 
can be taken as truly random. This may require some care, since it 
was seen in Section 4.1.4 that ephemeris errors have an associated 
periodicity. 


Update of Gyro Drifts 

This application is sufficiently important that it forms 
the basis for one of the positioning methods being investigated in this 
study. GCP can be used for this purpose. Problems can occur if it is 
necessary to update the gyro data at some regular time interval. Clouds 
and irregular land-mass distribution prevent any guarantee of control- 
point acquisition at a desired interval. Also, the method loses much of 
its attraction if it is necessary to perform a GCP update simply to 
provide smoothing-data continuity. To be of most benefit, an update 
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should be necessary only when in the vicinity of a series of images for 
which positioning is desired. 

There are two cases where GCP may be advantageous for 
providing gyro updates. The first would result from selection of gyros 
with very low drifts. If a regular update is required only once every 
several orbital passes, or only once per day, GCP could dependably pro- 
vide the data at very low system cost. In assessing the gyro potential, 
the work on gyro selection in references (24) and (25) would have been 
more valuable to the present study had it not been so closely connected 
with the SIMS concept, as mentioned in Section 4*2. Of considerably 
more interest here is the selection of a gyro which requires the least 
frequent update, even at some greater penalty in cost or weight. 

The second case applies if it is necessary to use GCP 
gyro updates only for those passes in which image data is to be posi- 
tioned. Methods GIMS-E-2 and GIMS-E-3 in Section 2 are based on this 
situation. One image may be positioned with GCP in the center of an 
imaging sequence; the updated gyros at that point then would supply 
attitude data in both directions along the arc. If extrapolation of 
attitude data in this way is not possible, it may be feasible to use 
two sets of update GCP at the ends of an image sequence, interpolating 
attitude data between the resected end images* This interpolative 
scheme is outlined in Section 6.2 of reference (25) as feasible for one 
particular gyro selection. 

Any of the implementations discussed here has the virtue 
of greatly reducing the number of GCP required to position EOS images, 
for a given required accuracy. This is important, since the biggest 
difficulty foreseen in applying a GCP-positioning method to EOS is in 
implementation of the great number of descriptive and correlative points 
required. However, more work is needed to properly assess the potential 
for this application of GCP to positioning EOS images. In reference (25), 
the potential advantages of the GIMS-E method M landmark-inertial" method 
were recognized. But the existence of a considerably body of knowledge 
based on "landmark" techniques is not recognized* This shows once again 
the need for close contact between "landmark" and "inertial" disciplines. 
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SECTION 5 

ANALYSIS OF POSITIONING METHODS 


Section 5.1 provides an error analysis for each of the five 
positioning methods defined in Section 2. The information developed 
in previous sections was drawn upon for the numerical values attached 
to each error contribution. Section 5.2 develops the significant 
points of difference for the methods, and identifies the area in which 
further work is needed to give a clear choice. 

5.1 ERROR ANALYSIS 

In all the analyses which follow, several simplifying assump- 
tions have been made. The accuracy indicator of concern has been de- 
fined as the rms positional error in the Thematic Mapper image which is 
not repeatable from one pass to the next. Errors introduced by film 
printers or other data processing equipment have been ignored, as have 
the effects of computational simplification and round-off. Mechanical 
errors such as reference-mark mismatch also have been ignored for all 
methods. Components for which no data are available have been esti- 
mated; where uncertainty exists, an effort was made to be liberal 
for satellite-data accuracy and conservative for GCP accuracy. 

Once again, a cautionary note: the values below do not repre- 

sent the entire error budgets. 

5.1.1 Method SIMS-E 

The component error contributions for method SIMS-E are 
shown in Table 5-1. This completely satellite-data method uses ephem- 
eris data for position, a gyro package for short-term attitude behavior, 


Table 5-1 - RMS Positional Errors, Method SIMS-E 


SIMS Attitude Determination: 0.001° each axis, 

or 17 m yT 24 m 

Attitude Transfer, Gyro Package to Imager: 

0.00028° to 0.00056° each axis 7 - 24 m 

Ephemeris Repeatability: 20 m 

Internal Sensor Errors: 0.25 to 0.5 IFOV 11 - 22 m 

Root Sum Square Error 34 - 41 m 
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and a star sensor to be used occasionally to correct gyro drifts. The 
error contributions show that there is no particular pacing error; to 
that extent the system is well-balanced. The attitude-transfer esti- 
mates assumed that some form of recording autocollimator is attached to 
the different sensor packages. Internal Thematic Mapper geometric 
errors are taken as one-fourth to one-half the IFOV. This value — re- 
peated in the other analyses — represents the root sum square total 
of detector misalignment, scan-wheel wobble, timing errors, and any 
other effects which are not mathematically modeled in the solution and 
are unlikely to be determined and calibrated out during the mission. 

5.1.2 Method GCP-E 

The component error contributions for method GCP-E are 
shown in Table 5-2. This completely ground-control method uses approxi- 
mate ephemeris data for position, and ground control points for all 
other attitude errors needed to represent the spacecraft attitude be- 
havior over a 30-second image frame. The computed attitudes also remove 
the effects of ephemeris errors. The method also assumes the space- 
craft attitude control system is functioning such that attitude be- 
havior over 30 seconds can be represented by at worst a second-order 
curve in each axis, with deviations from this behavior of 0.001 degree 
rms. 


Table 5-2 - RMS Positional Errors, Method GCP-E 


Spatial Resection: (0.51 - 0.61) (16 - 35 m) 

Internal Sensor Errors: 

Attitude Departure from Second Order Model: 
0.001 degree each axis, or 17 m -Ji 

Root Sum Square Error 


8 - 21 in 
11 - 22 m 

24 tn 

28 - 39 m 


The error contribution from the spatial resection is de- 
termined from (1) the standard error of unit weight in a GCP before the 
resection, and (2) a multiplicative weight factor determined by the num- 
ber and distribution of GCP within the image. The positional standard 
error of unit weight is determined from the error contributions of 
ground coordinates, descriptive pointing, and internal image geometry. 
The effects of correlative pointing and equipment have not been included 
by proper design of the processing system, these effects can be made 
beneath significance. The errors from ground coordinates are those as- 
sociated with the selection of points from 1/24,000 and 1/62,500 maps. 
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The descriptive pointing error was developed in Section 4.3. The stand- 
ard error of unit weight, Pp, in a GCP before resection then is com- 
puted as: 

From Ground Coordinates: (0.3 mm at map 

scale, rms) 7 - 19 m 

From Descriptive Pointing: 10 - 20 m 

From Internal Sensor Errors 11-22 m 


Root Sum Square = Pp 16 - 35 m 

Note that it is not assumed here that the same GCP will always be used 
for positioning successive images. (If this could be done, the ground- 
coordinate error contribution could be deleted from considerations.) 

The weight factor with which to multiply pp is taken 
from Table 4-2 in Section 4 assuming that angles, rates, and accelera- 
tions all must be determined for the image. Depending on GCP numbers 
and distribution, the values range from 0.51 to 0.61 for the weight 
factor. 


5.1.3 Method SIMS-EG 

The component error contributions for method SIMS-EG are 
listed in Table 5-3. This method incorporates occasional GCP to remove 
ephemeris bias errors, together with the star-sensor and gyro attitude 
measurement of method SIMS-E. It is assumed that the use of GCP will 
also remove the attitude transfer error between the gyro package and 
the imager. For this method, the spatial resection errors replace the 
ephemeris errors. 


Table 5-3 - RMS Positonal Errors, Method SIMS-EG 


SIMS Attitude Determination: 0.001° each axis 

Residual Position Errors from Spatial 
Resection: 

Internal Sensor Errors: 


24 m 

8 - 21 m 
11 - 22 m 


Root Sum Square Error 28 - 39 m 


5.1.4 Method GIMS-E-1 

The component error contributions for method GIMS-E-1 are 
listed in Table 5-4. This method incorporates a spacecraft gyro package 
with method GCP-E. The gyros are needed if spacecraft attitude behavior 
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Table 5-4 - EMS Positional Errors, Method GIMS-E-1 


Spatial Resection: (0.36 - 0.42) (16 - 35 m) 

Internal Sensor Errors: 

Gyro Errors: 0.00028° to 0.00056° each axis 

(one to two arc seconds) or ( 5 - 10 m) V* 

Root Sum Square Error 


6 - 15 m 
11 - 22 m 

7 - 14 m 
14 - 30 m 


cannot be characterized to the precision required for GCP-E. For this 
method, only the attitude angles must be found in the spatial resection 
to update the gyro data. The multiplicative weight factor from Table 
4-2 can be taken as 0.36 to 0.42. The same gyro performance will be 
used here as was applied to methods SIMS— E and SIMS-EG. Attitude trans- 
fer will be done as part of the spatial resection composite angle de- 
termination, Error estimates for variations GIMS-E-2 and GIMS-E-3 cannot 
be included because of present uncertainty as to the effects of interpo- 
lated and extrapolated gyro attitudes. 

5.1.5 Method HIMS-E 

The component error contributions for method HIMS-E are 
listed In Table 5-5. This satellite-data method incorporates the 
lower-accuracy attitude sensors associated with the attitude control 
system, together with ephemeris data for position. One attitude-trans- 
fer is shown for this method, reflecting the possibility that the hori- 
zon sensor and gyros may be hard-mounted to the same structure. How- 
ever, the attitude transfer between the control-system package and the 


Table 5-5 - EMS Positional Errors, Method HIMS-E 


Attitude Determination: 0.005 to 0.02° each 




axis, (87 - 350 m) ^2 

120 - 

o 

00 

m 

Attitude Transfer of Gyros to Imager: 0.0017° 




each axis 


70 

m 

Gyro Errors: 0,001° each axis 


25 

m 

Ephemeris Repeatability : 


20 

m 

Internal Sensor Errors 

11 - 

22 

m 

Root Sum Square Error 

135 - 

480 

m 
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image sensor is taken as 0.0028 degree (ten arc seconds) instead of the 
smaller values used for the other methods. It is unlikely an autocolli- 
mator would be used with the HIMS-E method and as is seen in the table, 
there would not be significant improvement using such a device. Finally, 
the gyro performance is not assumed to be as high quality as for the 
other methods above, again with little effect in the overall error. 

Unlike the situation for the previous methods, method 
HIMS-E shows an obvious pacing error component. The values for attitude 
error assume some additional data processing to achieve these accuracies. 

5.2 COMPARISON OF POSITIONING METHODS 

The five positioning methods can be compared in many ways. Some 
of the criteria mentioned in the previous sections of this report are 
listed in Table 5-6. Even if the proper information were available, 


Table 5-6 - Possible Comparison Criteria to Apply 
to Different Positioning Methods 

1. Accuracy 

2. Cost 

3. Equipment Reliability 

4. Speed 

5. Impact on Spacecraft Weight and Power 

6. Self-Checking Operation 

7. Ability to Upgrade Accuracy During Mission Life 

8. Personnel Requirements 

9. Lead Time for Design 

10. Capability for Future Improvement 

11. Monitoring Behavior of Other Equipment Components 

12. Constraints on Spacecraft Structural Configuration 
and Component Placement 

13. Attitude Control System Requirements 

14. Computational Load of Ground Data Processing 

15. Previous Experience 


and it certainly is not, it would be futile to rank each of the methods 
studied according to this or any other list of performance trades. The 
weights to be attached to the criteria are missing; these can only be 
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supplied by NASA, and only with consideration of all aspects of the EOS 
program. Any method or methods recommended by BRL would be meaningless. 
What can be done in the way of comparison is to attempt to draw together 
the different characteristics of the different methods and summarize 
them. 

The five methods differ primarily in their reliance upon 
the different potential information components given in Table 2-1. A 
qualitative expression of this reliance y from zero to ten, is given in 
Table 5-7. Also shown is the reliance upon the spacecraft attitude 


Table 5-7 - Reliance of Positioning Methods 
on Different Components 


METHOD 


COMPONENT 

SIMS-E 

GCP-E 

SIMS-EG 

GIMS-E-1 

GIMS-E-2 

GIMS-E-3 

HIMS-E 

GCP 

0 

10 


10 

10 

10 

0 

Ephemeris (beyond 
that needed for 
tracking) 

10 

0 

5 

0 

3 

5 

10 

Gyro Package 

10 

0 


5 

. .1 



10 

Stellar Sensor 

10 

0 


0 


0 

0 

Horizon Sensor 

0 

0 

0 

0 

0 

0 

10 

ACS 

0 

10 

0 

0 

0 


0 


control system, which was used only for method GCP-E of those examined. 

As the ACS designers respond to the requirements given in Section 4.2.2, 
the methods of interest may have to be modified accordingly. The inter- 
mediate values result for those cases where only relative accuracy, low 
error rate, or near-linear error rate is demanded, or where some back- 
up is available in case of failure. 
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SECTION 6 

CONCLUSIONS AND RECOMMENDATIONS 


The discussion in Section 4.2 on ACS design has shown that control- 
system performance capabilities must be more extensively considered 
before deciding on the recommended image-positioning method for EOS. 

It is significant that, of the five methods selected for detailed 
analysis, the most attractive from the cost— accuracy aspects appear 
to be those which impose some conditions on ACS or gyro— package per- 
formance . The most important conclusion of this study thus is con- 
sidered to be that ACS design must be considered as a part of the 
image positioning problem for EOS, along with image sensor design and 
ground image processing system design. 

The second major conclusion applies only to the methods analyzed 
in detail for EOS. For these five methods, when all reasonable criteria 
are considered, the use of ground control points shows either equality 
or clear advantages over satellite-data methods in every area but one: 
speed of ground data processing. Here there is not a clear disadvantage, 
but a concern over the need for human decision-making in a high-produc- 
tion data-processing environment. The hybrid positioning methods SIMS- 
EG and GIMS-E-2 and -3 manifest this concern; the less frequent use of 
GCP with these methods means less emphasis is placed on speed in the 
manual operations and man-machine interface. The ability of the ACS 
to provide smooth attitude data enters directly into this controversy. 

The longer the interval needed between attitude updates, the more attrac- 
tive the GCP method looks in comparison with a stellar sensor. The 
efficiency of GCP selection and matching can only be evaluated for a 
specific design implementation including equipment costs, personnel 
requirements, training, and so on. Experience being gained with the 
ERTS Precision Processing Subsystem has been and will continue to be 
of value in assessing GCP problems and potential; some modifications 
that would be of benefit for EOS have already been seen. Thus, the 
second conclusion: with suitable efficiency for GCP selection and match- 
ing activities during data processing, extensive reliance should be placed 
on ground control points" for positioning the images obtained for EOS and 
similar programs . Part of this efficiency can be attained from implemen- 
tation of the operations to be performed in the ground image processing. 
ACS design can also contribute to the efficiency by requiring only infre- 
quent GCP updates of attitude values. The extent to which auxiliary gyro 
data are required follows from the attitude behavior characteristics in 
the ACS design. 

In addition to these major conclusions, a number of other findings 
from the previous sections are summarized in the following list: 
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1. Some GCP will always be needed in a satellite image program 
requiring positioning or registration. Applications include 
(a) measuring angular misalignment of image and attitude- 
measurement sensor at intervals after launch, (b) monitoring 
internal sensor geometry distortions, (c) measuring angular 
bias in inertial attitude reference devices, and (d) periodi- 
cally assessing overall positioning/registration accuracy. 
These operations should not be done by an outside agency or 
investigator, but as an operations-dedicated part of the pro- 
cessing facility. 

2. Offset pointing offers a potential for radically decreasing 
the interval between EOS image collection for selected areas. 
This capability should be implemented so that pointing is 
done by integral ground-swath stepover. These is no need 

to effect extremely precise pointing; accurate measurement 
of the pointing angle may be necessary, depending on the 
image positioning method- Stationary pointing to a single 
ground location is not needed. During and immediately after 
any image-sensor offset pointing motions, no imaging should 
occur for any high-resolution imager on the spacecraft until 
transients have decayed suitably, 

3. Reaction-wheel unloading should never be permitted during 
imaging, 

4. No special effort should be made before launch to boresight 
the image and attitude sensors; a higher-accuracy and more 
meaningful determination can be made after launch using 
GCP. The effort saved here should be spent on designing 
adequate stability in the spacecraft structure to support 
more dependable alignment of the various sensors. 

5. To assist users of EOS images as much as possible, all film 
images generated should as a minimum be corrected for all 
known internal geometric error; to the extent practicable, 
absolute positioning and registration are desirable. 

6. The sole justification for including RBV cameras on any EOS 
mission which includes either the Thematic Mapper or the 
HKPI is in one particular failure mode of the attitude 
control system. The RBV cameras will not improve Thematic 
Mapper positional accuracy. 

7. User agencies should be asked to define the extent of the 
coverage they will require, in addition to the ground resolu- 
tion and frequency. When these indicate operationally impos- 
sible amounts of information for processing, the users should 
be so advised. 


6-2 



8 . 


Determined efforts should be made to better define the neces- 
sary sensor IFOV for different ground resolution requirements; 
the ERTS MSS is expected to be useful in this regard* 

9. For image sensors requiring extensive digital or analog geo- 
metric manipulation for other reasons * stability of internal 
sensor geometry is of more concern than simplicity* 

10* In evaluating image-sensor designs for satellites, the costs 
for the geometric processing needed for the image data should 
be one of the factors considered. 

11* The need for attitude-transfer techniques, with consequent 
weight, power , telemetry, and reliability impact, should be 
given more detailed consideration for high-accuracy image- 
processing schemes involving auxiliary attitude sensors. 

12. For satellites carrying line-scanning imagers, smooth attitude 
behavior is more important that high-accuracy verticality; 
control-system designs should more fully reflect this require- 
ment . 

13* An attitude control system with high absolute accuracy is 
not needed for EOS. 
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APPENDIX A 

SIMULATION OF CONICAL SCAN GEOMETRY 
T. A. Eppes 


A.l SENSOR DESCRIPTION AND CONFIGURATION 

In an effort to examine the geometry of the conical scanner, a 
simulation of its sensing operation was developed during this study. 

The model was designed with as much flexibility as required to adequately 
explore various sensor-spacecraft configurations. For the purposes of 
describing the geometry of the sensor, it was assumed ideal in that no 
internal distortion was allowed, hence the simulation provided an external 
geometric treatment of conical-scan sensors. 

Figure A-l illustrates the basic operation of a conical scanner. 

It is characterized by some cone angle y 0 and active scan polar angle p. 
Data is collected only during the period when the scanner is "looking" 
at targets in the active scan region. This region is interrogated 
sequentially in time by the stack of detectors within the sensor. Fig- 
ure A-2 illustrates the track of the detectors in sensor coordinates as 
one swath is scanned. This figure is typical of all but one of the seven 
spectral bands. Six detector elements are shown, indicating a six line 
scan swath. For the seventh (far-infrared) band, the resolution cells 
are three times as large, hence each swath consists of only two scanning 
detectors. 

There are basically three sensor-spacecraft configurations: 

(1) forward scanning, (2) aft scanning, and (3) circular scanning. 

Forward scanning involves mounting the sensor such that the central 
axis of the cone lies in front of the spacecraft while the active scan 
region passes directly under the spacecraft. Aft scanning is identical 
to forward scanning except that the cone axis is now behind the space- 
craft and the swath curvature is opposite. Circular scanning means that 
the central axis of the cone passes through the nadir while the active 
scan area is located behind the spacecraft. Data is collected by samp- 
ling the output of the detector array at some rate that produces near 
adjacency in the resolution elements. Since the rotational motion of 
the scanner is constant, sampling is performed in terms of equal time 
or equal polar angle increments. Each data point is stored with enough 
binary bits to adequately represent the reflectivity of one instanta- 
neous field of view (IFOV) , the ground coverage of each data point being 
determined by the IFOV which is a function of the imaging optics and the 
physical size of the detector array. Data from each stack of six (two for 
infrared) detectors, although obtained at essentially the same time, are 
stored sequentially. 
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Figure A-l - Illustration o£ the Operation of the Conical Scanner 
in Sensor Coordinates 



Figure A-2 - Orientation of Detector Array During One Swath 



Because of the forward motion of the spacecraft, subsequent swaths 
interrogate new areas on the ground- These data form the basis for the 
final image which is composed of data from many differnt swaths. The 
required image format is assumed to be 185 km by 185 km which at a nominal 
equatorial altitude of 1000 km forces the scanner to have a cone angle of 
13.23 degrees if the active polar angle is +24 degrees. For the purpose 
of discussion, all specifications are based upon the use of a forward 
scanning sensor-spacecraft configuration. With a detector element angu- 
lar field of view of 44 yrad square, each image consists of approximately 
720 swaths. Correspondingly, each swath consists of about 4200 by 6 IFOVs. 
In order to produce adjacency at swath center for a latitude of 30 degrees, 
the sample period between stacks is 3.8 ysec. Hence, the time required 
to scan one swath is 0.032 second, and with a duty cycle of 80 percent, 
the time between identical points in subsequent swaths is 0.04 second. 

The following sections outline procedures whereby points in sensor 
space or earth coordinates may be mapped from one to the other. For the 
above-mentioned numerical data, examples of IFOV lay-out and swath struc- 
ture are shown. In addition, latitude-longitude grids are mapped into 
sensor coordinates. Section A. 5 at the end of this appendix provides a 
listing and operational explanation of the Fortran routines that were 
used to generate the data. 

A. 2 SENSOR-TO-EARTH DEVELOPMENT 

This section describes the basic equations associated with trans- 
ferring points in sensor space to their locations on the ground. Geo- 
metrically, points in image or sensor space will be defined in terms 
of some scanner polar angle and swath number, measured from some refer- 
ence swath number. No attention is given to the problems of geometric 
scanner distortion, detector array sampling, radiometric error due to 
IFOV distortion, or digital tape format. Hence, the problem entails 
the projection of the IFOV of the scanner characterized by some swath 
number, polar angle, and spacecraft location to corresponding points on 
the earth 1 s ellipsoid. Because the numerical values of the sensor 
parameters depend upon which sensor-spacecraft configuration is used, 
the treatment was developed for an arbitrary system configuration. 

For the purpose of this derivation, the conical scanner may be , 
thought of as emanating "rays'* that vary in space as a function of time. 
This convention in no way reflects an attempt to represent the physical 
operation of the sensor. Figure A-3 illustrates the sensor coordinate 
system where origination of the rays occurs. Based on the direction 
cosines of the ray projected from the scanner, and since the scanner 
position and orientation at future times, are known, the intersection of 
that ray with the earth may be computed. 

To achieve this end, however, the direction cosines must be trans- 
formed through several coordinates systems beginning in sensor space and 
ending in geocentric coordinates. The earth is assumed to he an ellipsoid 


A-3 



2 n 



Figure A-3 - Definition of Sensor Coordinate System 


rotating at constant rate and with shape unperturbed by surface features. 
The spacecraft position at time t is given in addition to the orbital 
parameters necessary to define its location at any later time. Orienta- 
tion elements (i.e., attitudes, attitude rates, and attitude accelera- 
tions) are given at some reference swath center, and all values needed 
prior to or after this time are calculated by extrapolating these values 
and assuming them constant over the time frame of interest. 

The development begins by defining the equation of the ray leaving 
the conical scanner in sensor coordinates shown in Figure A-3. The equa- 
tion .of this line is 


x jl + cot^ y o J 
cos n 




1 1 + cot^y I 
x _ -yL- 'oj 

A' sin n 


=, = 

B' 


zL 


cot y 
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The polar angle q may be expressed as a function of time by 


n = c x (t - t Q ) 


(A-2) 


where 

n = polar angle 

= rotational rate of scanner 
t = time 

t = time required to scan 1/2 swath 

The sensor coordinate system (x, y, z) must be rotated to account for 
spacecraft attitudes or built-in sensor offsets manifested as attitudes. 
The rotated coordinate system will be labelled (x* , y*, z*) coupled to 
sensor coordinates by a transformation matrix consisting of Eulerian 
derived matrix elements. Hence 


= A 


x 
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where 


A = transformation matrix 


and 


A 


I cos <p cos k - sin $ sin u sin k 
e e e e e 


-cos on sin k 
e e 


sin <i> cos k + cos 41 sin to sin k 
e e e' ;e e 


cos sin (c + sin sin to cos k 
e e e e e 


cos u> cos K 
e e 


s,in sin < - cos 4> sin ui cos k 

e e e e e 


-sin 6 cos w 
e e 


cos <f> cos 0) 


A- 5 



where 


$ = pitch at some time t 

k = yaw at some time t 
e 

w = roll at some time t 
e 

Two of the attitudes, roll and pitch, may be broken up into sepa- 
rately derived components. The pitch may be considered as consisting 
of a pitch offset and a space platform error. The pitch offset is used 
to generate forward, aft, or circular scanning. The space platform error 
is time-varying and may be expressed as a function of pitch at reference 
swath center, pitch rate, and pitch acceleration. Likewise, roll may be 
thought of as both offset and platform derived. The roll offset is used 
to generate scanning of ground swaths to the right or left of the present 
ground swath. The platform-derived roll is, of course, identical in 
nature to that of the pitch component. Consequently, the angles neces- 
sary to compute the transformation matrix of (A-3) are 


^ = ♦t + S o 


0) = to, + 0 

e t o 


(A-4) 


< = 
e t 
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time required to scan 1/2 swath 

time when platform attitudes are required 

$ = pitch at time t = t 
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<f) = pitch rate at time t - t 
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^ - pitch acceleration at time t = t 
E o 

eo = roll at time t = t 
o o 

m = roll rate at time t = t 
r o 

m = roll acceleration at time t = t 
u o 

k = yaw ati time t = t 
o ^ » o 

• 

k = yaw f rate at time t « t 
t : O 

I* 

k = yaw acceleration at time t = t 
c o 

<t> t = platform pitch at t 
u> t = platform roll at t 
= platform yaw at t 
B 0 = pitch offset of sensor 
= roll offset of sensor 
« e = total roll at t 
4> e - total pitch at t 
* total yaw at t 

Substituting the transformation of (A— 4) into (A— 3) and rearranging 
terms yields 


\ 


c i *’ -V c 3 z ' 

with 
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where 



element of A corresponding to the i’th row and the 
j’th column 


A’, B', C' = direction cosines of line in sensor coordinates 
defined by (A-l) 


The rotated (x f , y*, z') coordinate system must be translated to 
an earth-based system (x e> y e , z e ). The system (x 1 , y', z’) is defined 
as one whose z’ axis coincides with a perpendicular from the earth’s 
ellipsoid. Since (x e , y e , z e ) will be defined as a system tangent to the 
ellipsoid at the origination of the perpendicular, (x r , y’, z’) and(x e , 
y e , z e ) are related by 


x’ - x g y’ = y e 



(A- 6) 


where 


H = spacecraft altitude measured along perpendicular to ellipsoid 
se 

Figure A-4 illustrates the next transformation, that is, from local 
(x e> y e » z e^ to geocentric (U, V, W) coordinates. The U-axis coincides 
with the meridional plane associated with the spacecraft longitude at 
time t. The origin of the (U, V, W) system is located at the center of 
the ellipsoid. The (xe, ye» Ze) an ^ (U» V, W) systems are related by 
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a S Q instantaneous heading or azimuth of spacecraft measured from 


the south. 


= latitude of spacecraft measured at origin of local coordinate 
system (x q> y e> z e ) 


U — E cos <$ 
x se 


W = E (1 * e ) sin 4> 
x se 




a cos 4 


se , 


2 2 
1 - e sin 4 


, 1 

sel 


1/2 


yields 


e = eccentricity of the earth 
a = semi-major axis of the earth 
Substituting (A-6) and (A-7) into (A-5) and rearranging terms 


R. U + R_ - R 0 V + R, = R, W + R. 

1 2 3 4 5 6 


(A-8) 


where 
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b ii = element from the i’th row and j'th- column of the 
transformation matrix of (A-7) 

C 2 > C 2 > C 3 = coefficients associated with (A-5) 

Rewriting (A-8) in parametric form gives 


U - ^ «5 W + h - V 

V » (Rj W + E 6 - R 4 ) 


• ' (A-9) 


The equation of the ellipsoid is 


U 2 + 


2 

a 



-CA-10) 


where 


a * semi-major axis of the earth 
b = semi-minor axis of the earth 

Since the equation of the line projected from the sensor is given 
by (A-9), and the equation of the ellipsoid is given by (A-10), the inter- 
section of the two may be computed. Substituting (A-9) into (A-10) and 
solving for W gives 
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where 



In general, two values of W 0 will be found; only one is of interest. 
Since that point is the one directly beneath or in the neighborhood of 
the sensor, the value of W which is positive in sign will be selected. 
This forces calculations to bd done only in the northern hemisphere for 
the present mode. A small modification would allow operation over the 
entire ellipsoid. 

Using (A-9), U c and V D may be calculated; thus the point of inter- 
section at time t in geocentric coordinates is known. If terrain height 
variations are ignored, the latitude and longitude of the point are 
given by 
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In order to express the point in some tangent-plane coordinate 
system (x e , y e , z e ) , (A-7) may be used. The values used for $ se , L ge 
a se must be the same for all points to be expressed in that system. 

The above values are usually chosen as those corresponding to reference 
swath center. In any event,. the transformation of points in sensor 
space to their locations on the ground is now complete if all matrix 
elements are known at time t. 

To evaluate the transformation equations presented above, it is 
required that spacecraft latitude and longitude be known at all times. 
They are given at some initial reference point, and the orbital equa- 
tions may be used to extrapolate those values to some other time. Fig- 
ure A-5 is a diagram of the orbital motion of the spacecraft in two 
coordinate systems. It is seen that the orbit is retrograde with the 
sensor imaging during the descending mode. If the orbit is assumed 
circular, the orbital angle will vary linearly with time. The orbital 
angle is taken here as that angle subtended by two lines passing through 
the center of the earth (also the orbital center). One line passes 
through the terminator, or point of maximum latitude. The other passes 
through the spacecraft position at time t. The orbital angle may be 
expressed as 
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9 ' 
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where 

0' ■ orbital angle at reference swath center 
V q = tangential velocity of spacecraft 
R q = radius of spacecraft orbit 
t = time when 6 is desired 
t = time required to scan 1/2 swath 

The longitude of the spacecraft at time t is then 
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where 



G “ earth rotation rate 
o 

^seo ” ^ on 8 : *- tu ^ e spacecraft at reference point 

a = inclination of orbit measured from the poles 
seo 


Because of the transcendental nature of the equations, $ se and Hge cannot 
be computed directly and must be found iteratively* Nevertheless, con- 
vergence is rapid and excellent accuracy is attainable. Inserting some 
nominal value for H se , $ se is estimated by 
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A new value of H is computed from 
se 
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se 
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If the increment in H ge from the previous value is excessive (>1 meter), 
a return to (A-15) is in order and a new $ se is generated. This process 
continues until the increment in H se is acceptable. 

The only remaining quantity required in the transformation equa- 
tions is the spacecraft heading a Se , This is computed as the dot product 
of two unit vectors tangent to the ellipsoid at spacecraft nadir at time 
t. One vector denotes the velocity neglecting earth rotation while the 
other lies in the meridional plane or longitudinal plane of that point 
and is tangent to the ellipsoid. 

The velocity vector components in the coordinate system shown in 
Figure A-5 are 
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The meridional vector components may be expressed as 


x = -CCl • CC2 • sin ip 
m 
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The above components are not normalized to unity; hence the dot product 
of the two vectors has to be divided by the two vector lengths. The 
azimuth is then expressed as 
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Figure A— 6 is a plot of azimuth (although it includes earth rotation) as 
a function of latitude. Consequently, given some time t measured from 
the beginning of some reference swath, . the orientation elements at that 
swath center, and the spacecraft position at that swath center, the pro- 
jection of the line corresponding to the IFOV of the scanner at that 
instant on the ground may be computed. 

The only additional requirement analytically is providing a means 
of projecting all four corners of the (assumed square) IFOV for a given 
resolution element size and sample interval. That is, the input to the 
projection equations will be a lateral resolution element number and a 
corner number for the six-element detector array. Figure A-7 illustrates 
the requirement of a new cone angle and polar angle as a function of 
lateral resolution element number and IFOV corner number. It also es- 
tablishes the numbering convention for the detector array corners. The 
relationship between cone angle and polar angle changes for the case of 
"longitudinal" incrementation (in the along- track dimension) is shown. 
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Figure A-6 - Azimuth as a Function of Latitude 

A similar diagram could be drawn to show angular IFOV width in the tran- 
sverse dimension. The relations between the cone and polar angle of 
corners 1 and 3 are 
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The relationship between corners 3 and 4 or 1 and 2 is 
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With respect to (A-l), the values substituted for y 0 and n vary depending 
on which corner of the stack is of interest. Equations (A-17) through 
(A-20) represent only small changes in n and y 0 since A^ and A2 are the 
angular dimensions of the detector array (typically 44 trad by 264 prad). 
The input value of n for (A-17) or (A-19) is determined by the lateral 
resolution element number while y 0 is always the cone angle. 

A. 3 EARTH-TO-SENSOR DEVELOPMENT 

This section describes the equations and procedures associated with 
mapping points from the ground to sensor coordinates. The objectives are 
to determine (1) when a point was sensed, (2) after how many swaths from 
some initial swath did this occur, (3) what is the lateral resolution 
element number, and (4) what is the detector number within the array. 

The position of the spacecraft is given for some time t = 0. The lati- 
tude and longitude of the point to be projected are known. The unknowns 
mentioned above are determined using an iterative procedure whereby esti- 
mates of spacecraft location at sensing time are generated. Based on the 
errors in these estimates, new time increments are incorporated which 
eventually produce convergence to the time when the point was sensed. 

Initially, the location in sensor coordinates is computed for t = 0. 
The transformation is performed using the given location of the space- 
craft at that time. The point is transferred for two times, the beginning 
and ending of the swath connected with the initial spacecraft position 
shown in Figure A-8(a). Assuming a linear velocity in sensor coordinates, 
the necessary number of swath increments is calculated and the spacecraft 
is advanced accordingly. This process continues with varying degrees of 
refinement until the swath number where detection occurred is located. 

At that point, the lateral resolution element number and detector number 
are easily found by a successive iteration over some nominal range. 
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To begin, the point of known latitude and longitude is expressed 
in geocentric coordinates where the U axis passes through the meridional 
plane corresponding to the spacecraft location at that time. Hence 


U = E" cos $ cos (L - L) 
p se 


V = E" cos $ sin (L - L) 
p se 
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where 
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W = E" (1 - e ) sin $ 
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r , 2.2 ., 1/2 

[1 - e sin $] 


+ H 


a = semi-major axis of the earth 
e = eccentricity of the earth 

H = height of point above reference ellipsoid usually 
zero 


L = longitude of spacecraft at time of interest 

S6 

$ = latitude of ground point 
L = longitude of ground point 

(Up, Vp, Wp) = position of ground point in geocentric coordinates 

A transformation identical to that of Figure A-4 which expresses the point 
in local coordinates which are tangent to the ellipsoid at spacecraft 
nadir may be written as 




r 





X 



u 


u 


e 



p 


X 


y e 

= 


V 

p 

- 

0 


z 



w 


w 


e 



p 


X 


L J 



_ — J 


L- J 

J 


(A-22) 


A-22 



where 


and 
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a se ~ instantaneous azimuth of spacecraft measured from 


the south 

( x e > y e > z e ) != local coordinates components of the ground point* 


A transformation must now be performed to arrive at a sensor-based 
coordinate system. This is accomplished by translating the local system 
b y ^se (spacecraft altitude) and incorporating the Eulerian attitude ele- 
ments. Hence, 
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H g£ - spacecraft altitude along perpendicular 

(x"\ y 1 ' * , z ,M ) = sensor coordinate system components 

The Eulerian angles required in evaluating B are obtained by extrapo- 
lating the orientation elements from the initial spacecraft location 
as in (A-4). In addition $ se , H se , L se , and a se are computed in the 
same way as in the previous section. 

As previously stated, the ground point is projected for the begin- 
ning and end of a swath where detection is suspected. To facilitate 
handling in sensor coordinates, the values obtained by (A-23) are nor- 
malized in the z-direction. The two points found in sensor coordinates 
appear separated in space and time. Since, in general, the initial 
swath will not be the one in which detection occurred, some form of 
swath incrementation using the above space and time separation must 
be devised. 

The normalized sensor coordinate plane of z = 1 is divided into 
five zones as shown in Figure A-8(b). Zone one will consist of all 
areas where x >_ 0. Zone two will be that area where X 0 and X XMAX 
where XMAX is the x-coordinate of the foremost portion of the active 
scan area. It may be expressed as 

XMAX - sin y cos n + del 
' o max 


where 

del = sin y - sin (y - A- ) 

1 o o 1 

A- = angular size of detector array in along-track direction 
(A x * 6A 2 ) 

= cone angle 

n max = maximum polar angle (^24 degrees) 

If the initial transformation produces points which lie in zone 1, 
the swath increment is computed assuming a linear velocity in sensor 
coordinates. The swath increment is such that the time advance it pro- 
duces causes the points to lie in zone 2. The increment may be calcu- 
lated from 
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where 
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t QO “ time between identical points in subsequent swaths 
t = time required to scan 1/2 swath 

x 2* ! = sensor x-coordinate for point projected at the end of 
the swath where detection is suspected 

xj M - sensor x-coordinate for point projected at the beginning 
of the swath where detection is suspected 

Because the velocity is not linear and the estimate of the velocity is 
made over a small interval compared to the necessary incrementation, 
the point will not appear where it was predicted at the new times. How- 
ever, the increment needed to place the point in region 2 is not critical 
and over-incrementing by adding some fixed number of swaths to the cal- 
culated value is desirable. Since the motion of the spacecraft is in 
the +x’ M direction, the time increment will be positive, and the point 
at the times will be closer to the active scan region. 

If the point is found to be in zone 2 either initially or because 
of action taken by the above discussed procedure when zone 1 is en- 
countered, a different equation for swath incrementation is needed. 

In this base the swath increment is 
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Since the estimate of the velocity is in error, care must be taken to 
ensure that if the point tends toward a swath end, over-incrementation 
does not force the spacecraft past the true sensing time. For this 
reason, (A-24) is scaled to about 0.95. If this new time forces the 
point to again lie in zone 2, a swath incrementation of 1 is suggested. 
In this way, the point will not be passed over accidentally. 

If the point is found in zone 3, a similar equation to the above 
is used. In this case, the polar angle is estimated by projecting the 
line connecting the two points and finding its intersection with the 
foremost edge of the active scan area as shown in Figure A-8(a), The 
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polar angle n is estimated using the above procedure, and the required 
number of swath increments is 


No* Swath Increm. 
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where 

x = x-coordinate of the point in the active scan area where 
detection is suspected, from Figure A-8(a). 

To avoid over-incrementation, the above estimate is reduced by up to 3 
swaths* Each time beyond this stage when the point projects into zone 
3, the swath increment is made 1* 

As soon as either the beginning or end of a swath falls within the 
active scan area, the interswath mode is begun. The polar angle is 
easily estimated by computing it from the average point projection which 
corresponds to swath center. The lateral resolution element number is 
then known. To ensure detection, these numbers are decreased and "stepped" 
in equal time samples. Each time, the ground point is projected to sen- 
sor coordinates. At this stage, the point is projected at only one time, 
that time corresponding to the lateral resolution element number. At 
each resolution stack, the point is checked to see if it lies in the re- 
gion bounded by the 6 detectors. If not, the element number is increased 
by one (with a maximum search range of 20) and the above test is reap- 
plied. Since it is possible to lie within the active scan area and not 
be sensed during that swath, once the point is not detected using 
resolution-element searching, the swath number is increased by 1, and 
the above procedure is repeated. After three interswath mode attempts, 
the point is assumed not to have been sensed as a result of combined 
yaw/earth-rotation. During this period, the point is also checked to 
ensure that the y-coordinate is within the bound of the active scan 
region. 


To facilitate overall procedure, all points found to lie behind 
the spacecraft initially are forced into zone 2 by retreating the space- 
craft sufficiently. This is done to prevent the need for equally elab- 
orate convergence schemes for these points since convergence time could 
not be decreased sufficiently. Once the point is placed in zone 2, 
normal operation toward true detection time applying the procedure dis- 
cussed above is resumed. 

In general, the point will have been detected by this time, there- 
fore the time from initial spacecraft location, the lateral resolution 
element number, and the total number of swath increments are known. The 
only remaining unknown parameter is the detector number, 1 through 6, 
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Since all six are interrogated essentially simultaneously, the center of 
the detector lying closest to the sensor coordinate of the ground point 
at the instant of detection will be the detector number. In conclusion, 
all parameters necessary to determine the location of some point on the 
ground in the sensor data format are known, 

A. 4 CHARACTERIZATION OF CONICAL SCANNER GEOMETRY 

This section graphically illustrates the geometry associated with 
the conical scanner in terms of swath curvature, IFOV projection, and 
ground-point grid mapping. Based on the information provided previously 
and with the aid of the Fortran programs in Section A. 5, numerical ex- 
amples of conical scanner operation are shown for various spacecraft- 
ground-point configurations. ^ 

To meet the above objectives, the spacecraft or sensor parameters 
that must remain fixed are listed. It is required that the swath width 
on the ground be 185 km wide. The orbit is fixed at a radius of 9799 km 
and is circular. In addition the maximum active polar angle is 24 de- 
grees. The spacecraft will also be operated in a forward scanning 
mode since this configuration minimizes overall terrain effects and re- 
quires a lower duty cycle. Based on the above data, all other adjustable 
parameters will be determined. 

From the altitude, maximum polar angle, and required swath width, 
the cone angle is found to be 13.23 degrees. If the angular IFOV of 
each detector is 44 rad (design goal), the polar angle sampling interval 
to produce adjacency at swath center must be 0.1923 mrad at 30 degrees 
latitude. In order to match adjacent swaths at 30 degrees latitude, the 
time between identical points on subsequent swaths is 0.040225 second. 
With a scanner duty cycle of 80 percent, the time required to scan one 
swath is 0.032084 second. Since the maximum polar angle is 24 degrees, 
the sample interval of the detector array is 7.4 ysec. Hence there 
are 4340 samples per swath with each sample requiring 6 addresses to. 
the detector array. At this point all adjustable scanner parameters 
have been specified. 

The immediate point of interest is the macroscopic structure of 
one swath. Figure A-9 is a plot of the swath as it appears in a plane 
tangent to the ellipsoid at the reference swath center. The error in 
the z-direction varies radially from the spacraft nadir and is maximum 
at each end of the swath at a value of -720 m. This value is repre- 
sentative of the earth drop-off at swath end points. The swath thick- 
ness is 262 m, which is approximately that of the width of, the plotted 
line at the scale of the figure.. 

In addition to overall swath geometry, an important feature of 
the scanner, both radiometrically and geometrically, is its IFOV pro- 
jection. Figure A-10 illustrates the layout of adjacent resolution- 
element stacks for two . subsequent swaths at three regions in the swath. 
The spacecraft latitude was approximately 30 degrees. There are four 
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Figure A-9 - Swath Geometry in Local Coordinates 


possible situations which might arise with respect to the IFOVs: (1) no 

coverage, (2) single element coverage, (3) two-element coverage and, 

(4) three-element coverage. In Figure A-10, only four detector IFOV 
are annotated to illustrate the above categorized phenomena. It is 
apparent that the greatest single difference between the western and 
eastern ends is the effect of earth rotation which produces overlapping 
and underlapping, respectively. With increasing slant range, the con- 
stant angular-width IFOV projects a larger area on the earth. Since the 
sampling interval was picked to provide adjacency at swath center, there 
will be overlap at all other points with the most severe at swath ends. 
This effect is manifested by the four small rectangular areas of two- 
element coverage on the eastern and western sides. Note that the swath 
center is affected only by earth rotation which produces a shear of the 
detector arrays with no overlap. A unique effect caused by earth rota- 
tion on the western side is the appearance of small areas where three- 
element coverage occurs. 

Figure A-ll illustrates the layout of IFOVs at the terminator of 
81 degrees latitude for the same portions of the swaths. Comparison 
with Figure A-10 reveals that earth rotation perpendicular to the space- 
craft motion is not present, which is to be expected since the heading is 
due west. Nonetheless, overlap due to increased slant range is the same. 
However, since the earth rotation is parallel to but opposite spacecraft 
motion, an increase in along-track IFOV overlap is present. This effect 
is present at swath center as well, while the earlier shear is missing. 

Figure A-12 illustrates the layout of IFOVs at the equator for the 
three sections mentioned previously. It may be contrasted with Fig- 
ure A-10; most phenomena are present, but to a greater or lesser degree. 
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Figure A-10 - IFOV Geometry at 30° Latitude Under Normal Operation 
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The major feature lacking in Figure A-12 is the cell overlap along the 
horizontal edges of the IFOVs caused by earth rotation. Since the 
normal component to earth rotation is minimum at the equator, the dis- 
appearance of the overlap is expected. This figure represents a worst 
case insofar as earth-rotation shear effects are concerned. 

As a means of minimizing earth rotation effects, the sensor may 
be yawed to produce a more southerly heading. The required yaw varies 
with latitude and the effects cannot be eliminated at all points in the 
swath. Figure A-13 illustrates the effect of a 3.5-degree yaw which is 
that value required at 30 degrees latitude. The ends still exhibit 
earth rotation effects although the situation at swath center is ideal. 
Maintaining the yaw at 3.5 degrees, the IFOV layout at 81 degrees lati- 
tude was explored to see what adverse effects would occur. Figure A- 14 
illustrates what amounts to an overcorrected yaw since no transverse 
earth rotation components exists for the spacecraft at that latitude. 

At the equator, the 3.5 degree yaw compensation would reduce the over- 
and underlap, but not to the same extent as for 30 degrees latitude. 

As a final check of IFOV geometry, a 20 degree eastward roll off- 
set was incorporated at 30 degrees latitude. Figure A-15 illustrates 
the three sections of the swaths used in the previous figures. Earth 
rotation and slant-range expansion are severe in the far eastern side 
and at swath center with only moderate distortion at the near eastern 
side. Note that the near eastern side has more three-element coverage 
compared to Figures A-10 and A-14. 

Another area of investigation of conical scanner geometry on 
ground point mapping was the projection of two latitude- longitude grids 
to sensor coordinates. The analytical procedure developed in the pre- 
vious section was used in performing the transformation. Figure A-16 
shows how a grid of lat itude-longitudeliness would appear if the sensor 
data were to be printed sequentially in time across the swath and se- 
quentially in swath number in the along-track direction. The grid is 
centered at 30 degrees latitude, with equal spacing of about 28 minutes 
of arc in latitude and longitude. Figure A-17 is a representation of 
a grid mapping at 81 degrees latitude where latitude lines are much fur- 
ther spaced in distance than longitude lines, although the spacings are 
equal in arc minutes. It is apparent that the spacecraft heading has 
changed by about 75 degrees. The swath curvature still produces a bowed 
distortion in the grid as it would be printed. 

This section has discussed the basic imaging geometry of the conical 
scanner as well as providing graphical illustrations of swath structure, 
IFOV projection, and grid mapping. The overall interaction of the rota- 
ting earth with the conical scan sensor has been characterized. It Is 
hoped that this information will serve as an aid in estimating the degree 
of geometric correction required, and in understanding the problem of the 
geometric relationship between sensor and earth. 


A-32 





ROLL = 0° PITCH = 13.23° YAW = 3.5° 


10 METERS 


& n = 1.923 E-04 RAD. 
Figure A— 14 


*se = 81 ° L se = 90 ° 

- IFOV Geometry at Terminator for Yaw Compensation at 
30° Latitude 





s/c 



10 METERS 


a Roll Offset at 20° 


for 


in 

d 


S/C MOTION 



P-51-367-3 



A. 5 PROGRAM LISTINGS 


Printed teleprinter listings of the image-to-earth and earth-to- 
image computer programs are presented on the following pages. Each pro- 
gram is prefaced by a long explanatory paragraph. The program was written 
in Fortran for the SBC "Call 370" time-sharing service. 
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100 " THIS SUBROUTINE CALCULATES THE INSTANTANEOUS LOCATION 

110 "ON THE EARTH OF THE FOUR CORNERS OF THE IFOV OF A CONICAL 
120 "SCANNER. IN ORDER THAT THESE POINTS BE LOCATED* ALL S/C 
130 "ORIENTATION ELEMENTS (ATTITUDES* ATTITUDE RATES* AND ATT I - 
140 "TUDE ACCELERATIONS* AND ALTITUDE) MUST BE KNOWN AT SOME 
150 "TIME AMD THE NUMBER OF SWATHS AND RESOLUTION ELEMENTS LATER 
160 "OR EARLIER AT WHICH THE IFOV LOCATION IS DESIRED. THE PRO- 
170 "GRAM WAS DEVELOPED FOR A SPECIFIC SET OF ORBITAL PARAMETERS 
180 "ALTHOUGH THEY MAY EASILY BE ALTERED INTERNALLY. GIVEN THE 
190 "GEOMETRIC BOUNDS ON THE EARTH* THE CONICAL SCANNER PARAMETERS 
200 "WERE DERIVED USING THE ASSUMED ORBITAL DATA BY ITERATION WITH 
210 "THE PROGRAM. 

220 " AS PRESENTLY ARRANGED* THE ORBIT IS CIRCULAR WITH AN 

230 "ORBITAL RADIUS OF 7358288 M MOVING AT A CONSTANT TANGENTIAL 
240 "VELOCITY OF 7400 M/SEC. THE ORBITAL INCLINATION WAS ASSUMED 
250 "TO BE 81 DEG WITH THE S/C HEADING AT THE EQUATOR BEING EITHER 
260 "SW OR NE. THE EARTH WAS ASSUMED TO BE A ROTATING ELLIPSOID 
270 "WITH A PERIOD OF 23 HR* 56 MIN* 4 SEC AND A SEMI-MAJOR AXIS 
280 "OF 6378388 M. THE CONICAL SCANNER CURRENTLY HAS A HALF- 
290 "CONE ANGLE OF 13.23 DEG YIELDING A SWATH OF 185 KM ON THE 
300 "EARTH FOR AN ACTIVE SCAN ANGLE OF + OR -24 DEG WITH THE 
310 "ORBITAL DATA GIVEN ABOVE. 

320 ” FOR EACH OF THE VISIBLE BANDS* THERE ARE SIX DETECTORS 

330 " WHICH SIMULTANEOUSLY SENSE THE EARTH* CONSEQUENTLY* EACH 
340 "SWATH IS COMPOSED OF SIX RESOLUTION ELEMENT TRACKS. THE 
350 "ANGULAR WIDTH OF THE ELEMENTS IN THE ALONG TRACK DIRECTION 
360 "WAS ASSUMED TO BE 44 MICRORAD WHICH AT THE ASSUMED MEAN 
370- "ALTITUDE YIELDED A CELL 43.4 M LONG. THE CELL WIDTH IN 
380 "THE CROSS TRACK DIRECTION WAS DETERMINED BY THE SAMPLE RATE 
390 "OF THE DETECTOR ARRAY. THE SAMPLE RATE CAN ONLY BE DETER - 
400 "MINED IF THE TOTAL SWATH SCAN TIME IS KNOWN AND THE DUTY 
410 "CYCLE OF THE SCANNER IS GIVEN. THE TIME BETWEEN SWATHS WAS 
420 "DETERMINED BE ITERATIVELY FINDING THE TIME NECESSARY TO 3E 
430 "ABLE TO PRODUCE ADJACENCY IN THE ALONG TRACK DIRECTION AT 
440 "SWATH CENTER FOR SOME NOMINAL LATITUDE. THIS TIME WAS 
450 "FOUND TO BE .040105 SEC* AND COUPLED WITH A DUTY CYCLE OF 
460 " 80%* THE SWATH SCAN TIME WAS .016042 SEC. THE SAMPLE PER- 
470 " I OD WAS THUS 3.7 MICROSEC AT EQUAL INCREMENTS OF 1.923E-04 
480 "RAD POLAR ANGLE WHICH YIELDED 4357 SAMPLES PER SWATH BEFORE 
493 "MULTIPLYING 3Y THE 6 RESOLUTION ELEMENTS. 

500 " THE PROGRAM IS ORGANIZED SO THAT VARIOUS S/C OR SENSOR 

510 "CONFIGURATIONS MAY BE SIMULATED* E.G. CIRCULAR SCAN* PITCH 
520 "DOWN* OR PITCH UP. THE LATER TWO ARE COMMONLY REFERRED TO 
530 "AS NADIR SCANNING SINCE THE S/C NADIR POINT IS ALSO A DATA 
540 "POINT. IN ADDITION* ROLL OFFSET MANEUVERS MAY' BE PERFORMED 
550 "WHEREBY REGIONS TO THE LEFT OR RIGHT OF THE NORMAL SENSING 
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560 "AREA MAY BE INTERROGATED* HOWEVER* CARE SHOULD BE TAKEN TO 
570 "ENSURE THAT IF A PITCH DOWN OR PITCH UP CONFIGURATION IS 
580 "IN USE* THE PITCH VALUE IS MODIFIED TO FORCE THE SWATH 
590 "CENTER TO LIE ON AN AXIS TRANSVERSE TO SENSOR LOCATION 
600 "SINCE THE TRANSFORMATION MATRIX WITHIN THE PROGRAM USES 
610 "TRUE EULER ROLL* PITCH* AND YAW ANGLES. 

620 " AS PRESENTLY ARRANGED* INPUT TO THE PROGRAM IS DONE 

630 "THROUGH A DATA FILE LABELED 'FILL'. FIVE PARAMETERS ARE 
640 "PRESENTLY TRANSFERRED* SWATH NUMBER (1 OR 2>* RESOLUTION 
650 "ELEMENT NUMBER (+2178 TO -2178)* CORNER NUMBER OF THE IFOV 
660 "OF INTEREST (1 TO 4)* PITCH ANGLE FOR NADIR OR CIRCULAR 
670 "SCANNING* AND ROLL OFFSET. 

680 " ATTITUDES* ATTITUDE RATES* AND ATTITUDE ACCELERATIONS 

690 "ARE MOT INPUTTED ALTHOUGH BY MODIFYING THE 'READ' STATE- 
700 "MENT* THEY CAN BE. ON THE EARTH* FACING IN THE DIRECTION 
710 "OF INSTANTANEOUS S/C HEADING* THE CORNERS OF THE IFOV ARE 
720 "LABELLED AS 1 FOR THE LEFT-BACK* 2 FOR THE RIGHT BACK* 

730 "3 FOR THE LEFT FRONT* AND 4 FOR THE RIGHT FRONT. AS PRE- 
740 "5ENTLY ARRANGED ONLY TWO ADJACENT SWATHS MAY BE EXAMINED* 
750 "ALTHOUGH WITR SOME MODIFICATIONS MANY SWATHS MAY BE 
760 "TREATED 

770 " THE OUTPUT FROM THE PROGRAM CONSISTS OF 3 CATEGORIES* 

780 "D S/C DATA AT THE CENTER OF SWATH 1* 2) S/C DATA WHEN THE 
790 "POINT IN QUESTION WAS SENSED* AND 3) LOCATION OF THE POINT 
800 "ON THE EARTH IN TWO COORDINATE SYSTEMS. OUTPUT 1 OF THE 
810 "PROGRAM CONSISTS OF S/C LATITUDE* UNCERTAINTY IN THE ALT- 
820 "ITUDE* ALTITUDE* LONGITUDE* LONGITUDE OF THE ORBITAL PLANE 
830 "TERMINATOR* AND THE HEADING OF THE S/C. OUTPUT 2 CONSISTS 
840 "OF S/C LATITUDE* ALTITUDE* ORBITAL PLANE ANGLE MEASURED 
850 "FROM THE TERMINATOR* AND ALTITUDE UNCERTAINTY. OUTPUT 3 
860 "CONSISTS OF LATITUDE AND LONGITUDE OF THE EARTH BASED POINT 
870 " (XEE* YEE* Z EE) COORDINATES IN A CARTESIAN SYSTEM TANGENT TO 
880 "THE EARTH’S ELLIPSOID AT THE CENTER OF SWATH 1. ALL CAL- ' 
890 "CULAT IONS ARE DONE IN DOUBLE PRECISION MODE AND ALL DIS- 
90/7 "TAMCES ARE MEASURED IN METERS* ALL TIMES IN SECONDS* AND' 
910 "ANGLES IN RADIANS (EXCEPT FOR 10 WHEN DEGREES ARE USED). 

920 " 

930 " 

940 " DEFINITION OF MAIN VARIABLES 

950 " 

960 " 

970 "A= SEMI -MAJOR AXIS OF THE EARTH 
980 "D3R= CONVERTS DEGREES TO RADIANS 

990 "BETA=PI TCH VARIABLE USED TO PRODUCE NADIR SCANNING 
1000 "THETA=ROLL VARIABLE USED TO PRODUCE ROLL OFFSET 
1010 "ECCN=ECCENTRICITY OF THE EARTH 
1020 "GAMMA=HALF-C ONE ANGLE OF THE SCANNER 
1030 "ALP= COMPLEMENT OF ORBITAL INCLINATION 
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1040 
1050 
13 60 
1073 
1 080 
1090 
1 100 
1113 
1120 
1 1 30 
1 1 40 
1 1 58 
1 1 60 
1 1 70 
1 1 80 
1190 
1200 
1210 
1 220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1 300 
1310 
1323 
1330 
1 343 
1 350 
1360 
1 370 
1380 
1 390 
1400 
1410 
1420 
1430 
1 448 
1450 
1460 
1 470 
1 480 
1490 
1 500 
1510 


"VO=TANGENTIAL VELOCITY OF S/C 
”RO=RADI US OF ORBIT FROM CENTER OF EARTH 
"C 1 =ROTAT I ON RATE OF THE EARTH 

"C 3= CONSTANT USED TO DETERMINE RATE OF CHANGE OF ORBITAL ANGLE 
" (MEASURED FROM THE CENTER OF THE EARTH) 

"TO=TIME TO SCAN ONE-HALF SWATH 
"CO-ANGULAR SCAN RATE OF SENSOR 
"B=SEMI -MINOR AXIS OF THE EARTH 

"THO=ORBITAL ANGLE DEFINED AS THE ANGLE SUBTENDED BY TWO 
" . LINES, ONE DRAWN FROM THE CENTER OF THE EARTH TO THE . 

" TERMINATOR AND ONE DRAWN FROM THE CENTER OF THE EARTH 

" TO THE INSTANTANEOUS LOCATION OF THE S/C 

"XLSEO=LONGlTUDE OF S/C AT SWATH 1 CENTER 

"PHIS0=LATITUDE OF S/C AT SWATH 1 CENTER 

"HSEO= ALTITUDE OF S/C AT SWATH 1 CENTER 

”DELT=UNCERTA I NTY IN S/C ALTITUDE AT SWATH 1 CENTER 

"XLPKsLONG I TUDE OF TERMINATOR AT SWATH 1 CENTER 

"XKAPP 0= YAW AT SWATH 1 CENTER 

"PH I D=P I TCH RATE AT SWATH 1 CENTER 

"PHI DD=P ITCH ACCELERATION AT SWATH I CENTER 

"XKAPD=YAW RATE AT SWATH 1 CENTER 

"XXAPDD=YAW ACCELERATION AT SWATH 1 CENTER 

”T 00= ONE -HALF SWATH SCAN TIME OR TIME BETWEEN SWATHS DEPENDING 
" IF SWATH 1 OR 2 IS BEING EXPLORED 

"VO=ROLL AT SWATH 1 CENTER 
"WD=ROLL RATE AT SWATH 1 CENTER 


"ETA=POLAR ANGLE OF SCANNER 

»ELL=RESOLUTION ELEMENT NUMBER (-2178 TO +2178) 
"VDD=ROLL ACCELERATION AT SWATH 1 CENTER 


"T=T I ME MEASURED FROM THE START OF SWATH 1 
"PH I T= TIME VARYING CONTRIBUTION OF PITCH 
"XXAPT=T IME VARYING CONTRIBUTION OF YAW 
"WT=TI ME VARYING CONTRIBUTION OF ROLL 
"PHI=M£T INSTANTANEOUS PITCH 
"W=NET INSTANTANEOUS ROLL 
"XKAPPA=NET INSTANTANEOUS YAW 

"A ( I , J ) =TRANSFORMAT I ON COEFFICIENTS FROM S/C ORIENTATION TO 

" COORDINATE SYSTEM TANGENT TO THE ELIIPSOID 

"THT= INSTANTANEOUS ORBITAL ANGLE 

"XLS E= S/C LONGITUDE AT T 

"PHI SE=S/C LATITUDE AT T 

"HSE=S/C ALTITUDE AT T 

"B C I ,J) TRANSFORMATION COEFFICIENTS FROM LOCAL COORDINATES TO 
" GEOCENTRIC 


"XEE=X-COORDI NATE 
"YEE=Y- COORDINATE 
"ZEE-2 -COORDINATE 
"XLAT=LAT I TUDE OF 


OF POINT IN 
OF POINT IN 
OF POINT IN 
POINT BEING 


TANGENT PLANE SYSTEM 
TANGENT PLANE SYSTEM 
TANGENT PLANE SYSTEM 
SENSED 
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1 520 
1 530 
1 540 
1 550 
1 560 
1 5 70 
1 580 
1 590 
1 600 
1610 
1 620 
1 630 
! 640 
1 6^1 
1 648 
1 643 
1 650 
1 660 
1 670 
1 680 
1 690 
1 700 
1710 
1 720 
1 730 
1 740 
17 50 
1 7 60 
1 770 
1 780 
l 790 
1800 
1810 
1 820 


XLNG=L0NGITUDE OF POINT BEING SENSED 
ALPHA=HEADING OF S/C AT T 
AZO=H£ADING OF S/C AT SWATH 1 CENTER 


LATLNG AND ASE WITH LATLNG 
ALTITUDE OF THE S/C FROM 
ASE DETERMINING THE HEADING 
ALTITUDE* AND ORBITAL ANGLE 


THE TWO SUBROUTINES ARE 
’DETERMINING THE LATITUDE AND 
'SOME GIVEN ORBITAL ANGLE AND 
"THE S/C GIVEN SOME LATITUDE* 

REAL* 8 GAMMA 1 * ETA 1 * BARG* ELL 

REAL* 8 XLAT C! 8 ) *XLNG ( 1 8 ) , XEE<1 8 ) * YEE <1 6 > * ZEE C 1 8 ) 

REAL*^ A 9 DG» ft * EC CN > GA * % 

BETA* THETA* ALP* VO*R0* Cl * C3* TO * CO * B* TH0 * XLSEO, PH I SO * HSEO* % 


DELT* PH I SD* XLPR * AZO * XLPK* XLPKD* XX* X KAPP0 * PH I D*PHI DD* % 
XKAPD* XKAPDD* WO* IvD* WDD* T* ETA* TOO » CTX * CTY* CTZ*PH I T* XKAPT* % 
WT*PHI* W*XKAPPA,PHIE* CONI * CON 2* C0N3 * CON 4* CON 5* C0N6, A 1 1 *% 

A12*A13*A21*A22*A23*A31*A32*A33*CCX*CCY*CCZ*THT*THID*% 

XLSEP*XLSE* PH I SE*HSE* B11*B12*B13*B21,B22*B23*B31*% 

B32*B33*E,UX*'.>.!X*R1*R2,R3*R4*R5*R6*BIGA*BIGB*BIGC*WT *U1 *% 

V! *XLT,DIFF*AMT*P*EC0N2*EC0N*UP, VP, WP* UXX * WXX, ALPHA * HSEN 
A — 6378388 • 

CALL OPEN C 1 *’ FILL '*• INPUT ’ ) 

READ C 1 **> N CODE* ELL* NPT * BETA * THET A 

DGR=57. 29577951 31 

BETA=BETA/D6R 

THETA=THETA/DGR 

ECCN =0 .082 

GAMMA = 1 3. 23/DGR 

ALP = 9../DGR 

VO = 7400 • 

R0 = 7 358288 . 


OF 


1830 Cl =360./<861 64 . *DGR ) 

18 40 C3= V0/R0 

1 8 50 TP=0. 0160 42 

18 60 CO = 48 . / (DGR*0 • 0 3208 4 ) 

1870 B=A*DSQRTC1 .-ECCN**2) 

1880 TH0=60./DGR 
1890 XLSE0=9O./DGR 

1900 CALL LATLNG(TH0,ALP,A,B*R0*ECCN*PHIS0*HSE0*DELT) 
1910 PHISD=PHISO*DGR 

1920 30 XLPR = DATAN ( DSIN C THO ) / ( DCO S C THO )*DS IN C ALP ) > ) 

1930 CALL ASE(A,ECCN,ALP*C3*C1 *B*TH0*R0,PHISO*HSE0*AZ0) 
1940 XLPK=XLSEO-XLPR 
1950 XLPKD=XLPK*DGR 
1960 XX=XLSEO*DGR 
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1970 yRITE(6,*>PHlSD,DELT,HSE0,XX,XLPKD, AZ0 : ' f ' 

1980 XKAPPO = 0. 06128939391 l 5 ’ 

1990 PH I D=0 . 0 
2000 PH I DD=0 • 0 
2010 XKAPD=0.0 
2020 XKAPDD=0.0 
2030 TOO=TO 

2040 I F ( NCODE . EQ . 2 ) TOO = 0. 0561 47 

2050 W 0=0.0 

20 60 W D=0 « 0 ‘ 

2070 ETA=ELL* 1 . 923E-04+9 . 61 5E-05 

2080- DD=0 . 0 

2090 T=ETA/CO+TOO 

2100 ETA 1= ETA 

2110 GAMMA 1=GAMMA 

2120 IF (NPT . EQ . 4 ) GO TO 5V 

2 t3 0 ■ I F ( NPTrr 2 ) 50,51,53 

2140 51 BARG = DC OS ( GAMMA 1 ) * ( DCOS ( 4 • 4D-0 5 ) +DS IN ( GAMMA 1 >*DSINCETAl )*% 
2150 DSIN(4.4D-35>/DSQRT( 1 . -DS I N ( ETA 1 )**2*DSINC GAMMA 1 )**2)> 

2160 GAMMA= DATAN ( DSQRT C 1 . -BARG* *2 ) /BARG ) 

2 1 70 BARG = DTAN ( GAMMA 1 )* DCOS ( ETA 1 ) /DTAN < GAMMA > 

21,83 ETA= DATAN ( DSQRT (• 1 . - BARG** 2 ) / BARG ) 

2190 I F( NPT . EQ . 4 ) GO TO 53 
2200 GO TO 50 

2210 53 BAR G = DCOS ( GAMMA )*(DCOS(2*64D-04)+DSIN( GAMMA )*DCOS(ETA)*% 

2220 DSINC2. 64D-04) /DSQRT < 1 . -DCOS ( ETA ) **2* DS IN ( GAMMA ) **2 ) > 

2230 GAMMA1 =DATAN < DSQRT ( 1 » -BARG**2 ) /BARG ) 

2240 BARG = DTAM ( GAMMA ) *DS I N ( ETA J /DTAN ( GAMMA 1 > 

2250 ETA= DATAN (BARG /DSQRT ( 1 . -BARG**2 ) ) 

2260 GAMMA= GAMMA 1 

2270 50 I F C C ETA 1*ETA> .LE.0.0) ETA = -ETA 

2280 NSWAT= 1 

2290 NAB= 1 

2300 NSMAX=1 

2310 NL00P= 1 

2320 46 IF(NLOOP-l) 39,39,45 
2330 39 CONTINUE 

2340 48 ARG= DSQRT (1 .0+(DCQS( GAMMA ) /DS I N < GAMMA ) ) **2 ) 

2350 CTX=DCOS(ETA)/ARG 

2360 CTY=-1 . 3 *DS IM( ETA > /ARG 

2370 CTZ= ( DCOS ( GAMMA ) /DS I N ( GAMMA ) ) /ARG 

2380 ?HIT=PHISE0+PHID*(T-T0)+SIGN*PHIDD*(T-T0)**2/a. 

2390 XXAPT=XKAPP0+XKAPD* C T-TO ) +XKAPDD* < T-TO) **2 /a . 

2400 yT = y0+VD*CT-T0> + VDD*(T-T0)**2/2. ■■ 

2410 PH I = PH I T+3ETA 
2420 W=VT+ THETA 
2430 XKAPPA=XKAPT 
2440 PH I E=PH I 
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2450 COM 1 = DCOS ( PH IE) 

2460 C0M2 = DSIN( PH IE) 

2470 COM 3 = DCOS (XKAPPA ) 

2480 C0N4=DS IN (XKAPPA ) 

2490 COM 5=DC0S < W ) 

2500 C0N6=DSIN(W) 

2510 All = CON 1 * COM 3 -COM2* CON 6* COM4 
2520 A1 2=-C0N5*C0N4 
2530 A1 3-C0N2* COM3+ CONI *C0N6*C0N4 
2540 A2 1 = CON 1 *C ON 4 + COM2* CON 6* COM 3 
2550 A22=C0N5*C0N3 

2560 A2 3= COM2* CON 4- CON 1 * CON 6* COM3 : , 

2570 A3 1 =-COM2*COM5 
2580 A32=C0N6 
2590 A33=C0M1 +C0N5 

2600 CCX=(A1 1 /CTX-A2 1 /CTY ) / ( A2 3/ CTY-A 1 3/CTX ) - ( A 1 l/CTX-A31/% 
2610 CTZ)/(A33/CTZ-A1 3/CTX) 

2620 CCY=(A22/CTY-AI2/CTX)/ (A23/CTY-A1 3/CTX >-(A32/CTZ-A12/% 
2630 CTX)/(A33/CTZ-A1 3/CTX) 

2640 CCZ=CCX*( (A33/CTZ-A23/CTY )/ (A32/CTZ-A22/CTY ) - (A23/CTY% 
2650 -A I 3/CTX > / (A22/CTY-AI 2/CTX) ) / ( (A2 1 /CTY-A3 1 /CTZ ) /(A32/% 
2660 CTZ-A22/CTY ) - (A1 1 /CTX-A2 1 /CTY) / CA22/CTY-A1 2/CTX ) ) 

2670 THT=THO+<VO/RO)*(T-TO) 

2680 TK 1 D=THT*DGR 

2690 XLSEP = DATAN < DS IN ( THT ) / ( DCO S ( THT ) * DS I M C ALP ) ) ) 

2700 XLSE-XLPK+XLSEP+ C1*(T-T0) 

2710 IF(XLSE-3. 14159) 10,10,11 

2720 11 XLSE=XLS£-6 .2832 

2730 10 CALL LATLNG < THT, ALP, A, B, RO , ECCN, PH I SE, HSE, DE.LT ) 

2740 PHISD-PHISE+DGR 

2750 2 WBITE( 6, * )PKISD,HSE, THI D, DELT 

2760 CALL ASE ( A, EC CM, ALP, C3, Cl ,B, THT, P.0, PH I SE, HSE, ALPHA ) 

2770 B l 1 = DCOS (ALPHA ) *DS I M (PH I SE ) 

2780 31 2 = - DS 1 1J (ALPHA ) 

279 0 31 3 = - DCOS C ALPHA ) *DC OS (PH I SE ) 

2800 B2 1 = DS I N ( ALPHA )*DSIN(PHISE) 

2810 322=DC0S (ALPHA) 

2820 B2 3= -DS IN (ALPHA)* DCOS (PH I SE) 

2830 B3 1 =DCOS (PHI SE) 

2840 B 3 2 = 0 . 0 

2850 S33 = DSIN( PH I SE ) 

2860 E = A/ DSQRT ( 1 • 0 -ECCN** 2* DS I N ( PH ISE)**2) 

2870 UX = E*DCOS ( PH I SE) 

2880 '.;X = E* ( 1 . 0-ECCN**2)*DSIN(PHISE) 

2890 R1=(CCX*B1 1 - CCY*B2 1 ) / ( CCY+82 3 - CCX*E 1 3 ) - ( CCX*B1 1 -CCZ*B31 )% 
2900 /( CCZ*B33-CCX*B1 3 ) 

2910 R2 = -CCX*(B1 1 *UX+ B 1 3*WX)/( CCY*323-CCX*B 1 3) + CCX * ( % 

2920 B 1 1 +UX+B1 3*WX)/(CCZ*B33-CCX*B1 3) 
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2930 R3=CCCY*S22-CCX*B12)/<CCY*B23-CCX*B13)-<CCZ*832-CCX*B12)% 
2940 / ( CCZ*B33-CCX*B1 3) 

2950 R4=-CCY*(B2 1 *UX + 82 3*WX ) / < CCY + B2 3 - CCX*8 1 3 ) + CCZ* % 

29 60 ( B31 *UX+B33*WX + HSE) /< CCZ*B33-CCX*B1 3) 

2970 R 5 = R 1 * ( C CCY*B23 - CCX*B I 3 ) / ( CCY+82 2 - CCX*B 1 2 ) - ( CCZ*B33-CCX*% 
2980 B13)/(CCZ*B32-CCX*B12) )/( (CCX*B1 1-CCY*B21 )/( CCY*B22-% 

2990 CCX*B 12)-(CCX*31 1-CCZ*S31 ) / < CCZ*B32-CCX*B1 2 ) ) 

3000 R6=R!*( C-CCY*C321*UX+B23*WX)+CCX*(B1 1 *UX + B 1 3* WX ) % 

3010 ) /( CCY*B22-CCX*B1 2 ) -( -CCZ* CB31 +UX+B33* WX+HSE) +CCX*% 

3020 ( B 1 1 *UX+B1 3*WX) ) /< CCZ*B32-CCX*B1 2) )/C CCCX*31 1 -CCY*% 

30 30 B21 )/CCCY*B22-CCX*B12)-(CCX*Bl 1 -CCZ*B31 ) / ( CCZ *B32 -CCX*% 

3040 B12)>+R2 

3050 81 GA= 1 . 0/B**2+CR5**2/A**2)*< 1 . /Rl**2+1 ./R3**2) 

3060 31 GB= <2.*R5/A**2)*< ( R6-R2 ) /R 1 ** 2 + ( R6-R4 ) /R3**2 ) 

3070 BIGC-<R6-R2)**2/(A**2*R1**2)+<R6-R4)**2/ ( A**2*R3**2 ) - 1 .2 
3080 W1 =( -BIG8+DSQRT<BIGB**2-4.*BIGA*BI GC ) ) / < 2 . *B I GA) 

3090 U1 =< 1 ,/RI )*<R5*W1 +R6-R2) 

3100 VI =( 1 ./R3)*<R5*VJ1+R6-R4> 

3110 3 DI FF - DAT AM ( VI /U1 ) 

3120 XLT=XLSE- DI FF 
3130 DIFF=DIFF*DGR 

3140 NSUB= ( WSVAT- 1 )* 1 6+NL00P+ (NAB- 1 >*8 
3150 XLNG ( WSU3 ) =XLT* DGR 

31 60 AMT=W1 /DSQRT (A**2* C 1 . -ECCN**2 )**2+ECCN**2*Wl **2 ) 

3170 P = DATAN < AMT/ DSQRT ( 1 . -AMT**2 ) ) 

3180 XLAT(NSUB)=P*DGR 

3190 EC0fJ2=A/DS0ET ( 1 . -ECCM**2*DSIU<P ) **2 ) 

3200 ECON=A/DSQRT< 1 . - ECCM**2*DS IN (PH I SO ) **2 ) 

321 0 UP = EC0N2* DCOS (P)*DC0S( XLSEO-XLT ) 

3220 VP=EC0N2*DC0S(P)*DSIN(XLSE0-XLT> 

3230 WP=EC0N2*< 1 • -ECCN**2)*DSIN(P) 

3240 UXX=ECO!J*DCOS (PHI SO) 

3250 WXX=ECON*C 1 . -ECCN**2 > *DS IN (PHI S 0 ) 

3260 31 l=DCOS(AZO)*DSIN(PHISO) 

3270 B12=-DSIM(AZ0) 

3280 B1 3=-DC0S(AZ0)*DC0S(PHIS0) 

3290 B2 1 =DS IN ( AZ 0 ) *DS IN (PHI SO ) 

3300 B22=DC0S(AZ0) 

3310 B2 3= - DS I N ( AZ 0 ) * DC OS < PH ISO) 

3320 B3 1 =DC0S ( PH ISO) 

3330 B33=DSIN (PHI SO) 

3340 XEE(NSUB) =B1 1 * C UP-UXX ) +B 1 2* VP + 3 1 3* ( VP- WXX ) 

33 50 YEE(NSUB)=B2 I * C UP-UXX ) +B22* VP + B2 3* C UP- WXX ) 

33 60 ZEECNSUB ) =B3 1 * ( UP-UXX ) +B32+VP + S3 3* ( WP-WXX ) 

3370 T=T+0. 022222222 
3380 NLOOP=NLOOP+l 
3390 GO TO 46 

3430 45 VRITE(6,*)XLATC 1 ),XLNGC 1 )^XEE( 1 )^YEE( 1)>ZEE<1 > 
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3410. STOP 
3420 END 

3430 SUBROUTINE LATLNG ( THT, ALP, A, B, RO, ECCN, PH I SE, HSE, DELT ) 

3440 REAL *8 THT, ALP, A, B, RO, ECCN, PH I SE, HSE, DELT, KSEN 
3450 H S EN = 979900. 

3460 NCT=1 

3470 1 PH I SE= DAT AN ( < ( A**2+A*HSEN) /< B**2 + A*HSEN> ) * C DCOS ( THT >*DCOS% 
3480 (ALP)/DSQRT(DSIN<THT>**2+DC0S(THT)**2*DSIW(ALP>**2> ) ) 

3490 HSE=RO*DCOS(THT)*DCOS(ALP)/DSlN(PHISE>-A*( 1 . -ECCN**2 >/DSQRT% 
3500 ( 1 • -ECCN**2*DSIN(PHI SE) **2 ) 

3510 DELT = DABS CHSEN -USE ) 

3520 I F( DELT- 1.0) 2,2,13 

3530 13 I F ( NOT -20 ) 14,2,2 

3540 14 HSEN=HSE 

3550 NCT=NCT+1 

3560 GO TO 1 

3570 2 RETURN 

3580 END 

3590 SUBROUTINE ASE< A, ECCN, ALP, C 1 , C2, B, THT, RO, PH I S E, HSE, AZ ) 

3670 REAL *8 XO, YO, ZO, A, ECCN, ALP, Cl , C2, B, THT, RO, PHI S E, HSE, AZ, % 

3610 EPP, PHO, X, Y , Z, UARG, C 1 1 , C22, C 1 2, C2 1 , C 1 3, C23, PHI SET, PHOT, % 

3620 EPPT ,XT , YT, ZT, CCC1 , CCC2, ZETA,XETA, YETA,P 1 ,P2,P3, DOT 
3630 XO-R.O* DS I N ( THT ) 

3640 Y O=R0*DC0S (THT ) *DS I N ( ALP ) 

3650 ZO=RC*DCOS<THT)*DCOS(ALP) 

3660 EPP=A/DSQRT ( 1 . - ECCN**2* DS I N ( PH 1 S E) **2 ) 

3670 PHO= DATAN ( DS I M ( THT ) / ( DCOS C THT ) *DS I N C ALP ) ) ) 

3680 X=EPP*DCOS(PHISE)*DSIN(PHO) 

3690 Y= EPP* DC OS (PHI SE)*DCOS (PHO) 

37 0 0 Z= EPP* (1 . -ECCN** 2 )*DSIN(PHISE> 

3710 UARG= ( CA**2+A*HSE) / <3**2+A*HSE) )* ( DCOS ( THT )* DCO S (ALP ) /DSQRT% 
3720 ( DS IN ( THT )**2 + DC0S( THT ) **2* DS IN (ALP > **2 ) ) 

3730 Cl 1-1 .0 
3740 C22= I • 0 

37 50 Cl 2 = R0*DC0S< THT)* DCOS ( ALP) * DC OS (PH I SE > /DS I N ( PHI SE ) **2+ A* (% 

37 60 1 .-ECCN**2)*ECCN*ECCN**2*DSIN(PHISE)*DC0S-(PHISE>/(DSQRT( 1 

3770 ECCN**2*DSIN(PHISE)**2) )**3 

3780 C2 1 = ( 1 ./(l .+UARG**2> ) * ( A* ( A* *2 -B**2 ) / ( B** 2+HSE*A ) **2 ) 

3790 Cl 3= -RO*DS INC THT ) *DCQS (ALP ) / DS I N ( PHI S E ) 

3800 C2 3= ( ( A**2+A*HSE ) / ( 3**2 + A*HSE) )*(!./(! * + UARG** 2 ) )*(-DC0S% 

381^ (ALP ) *DSIN (THT ) / ( DSQRT ( DS I N ( THT ) **2 + DCOS ( THT ) **2* DSIN ( ALP ) % 
3020 **2))**3) 

3330 PHISET=C(C13-C23*CI 1/C21 > /( C 1 2 -C22*C 1 1 /C2 1 ))*C1 

3840 PHQT=C 1 *DS IN ( ALP ) / ( DCOS ( THT ) ** 2*DS I N (ALP ) ** 2 + DS I N ( THT ) % 

38 50. **2 ) 

3860 EPP T= A* EC CN* * 2 * DS I N ( PH I S E > * DC 0 S ( PH I S E ) *?H I S ET / DS QRT % 

3870 ( ( 1 . -ECCN**2*DSIN(PHISE>**2>**3) 

3880 XT=DC0S ( PH I SE) *DS IN ( PHO > *EPPT+EPP* DCOS (PK ISE)*DC0S (PHO)*PHOT% 
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3890 

3900 

39 1 0 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 

40 10 
4020 
40 30 
4040 
4050 


-EPP*DSIN (PHO)*DSIN<PHI SE)*PHISET 

YT=DCOS<PHISE)*DCOS<PHO>*EPPT-EPP*DCOS(PHISE)*DS1N(PHO)*PHOT% 


- EPP*DCO$ (PHO ) *DS I N(PH I SE) *PHI SET 
ZT=(1 .-ECCN**2)*(DSIN(PHI SE) *EPPT+ EPP*DC OS ( PHI SE ) *PH I SET ) 
CCC 1 ( A**2/B**2 ) *Z/ ( X* DSIN < PHO ) +Y* DCOS (PHQ ) > 


CCC2 = DSQRT C 1 . / < 1 . + CCC 1 **2 ) > 
Z ETA= - CCC2 

XETA=-CCC1 *DS1N(PHQ)*CCC2 
YETA= - CCC 1 *DCOS (PHO ) *CCC2 


P 1 =XETA*XT 
P2=YETA*YT 
P 3 -£ ET A 5^c £ T* 

DOT= (P1+P2+P3) / DSQRT ( (XT**2+YT**2+ 2T**2 ) * <XETA**2+YETA** 2 + % 
ZETA**2>> 

A2=DATAM( DSQRT ( 1 .-D0T**2)/DQT) 


RETURN 

END 


1 00 
1 10 
1 20 
1 30 
1 4''* 
1 50 
1 60 
1 70 
180 

1 90 
200 

2 1 0 
0 9 0 

230 

2.40 

950 

260 

P70 

280 

290 

30 m 

3 i a 
320 


• THIS PROGRAM CALCULATES THE VARIOUS CONICAL SCANNER I HAUL 

•PARAMETERS NECESSARY TO PROJECT POINTS ON THE EARTH TO 
•POSITIONS IN THE SEQUENTIALLY PREPARED DIGITAL TAPE. 

•TO INITIALIZE THE PROGRAM* THE LOCATION OF THE S/C NEST 
"BE KNO'-’N AT SOME TIME ’0* ALONG WITH ITS OEENTATIOW ELE- 
MENTS (ATTITUDES* ATTITUDE RATES* AND ATTITUDE ACCELERA- 
TIONS). IN ADDITION* THE LATITUDE AND LONGITUDE OF THE 
•POINT TO 3E PROJECTED MUST BE GIVEN. THE OUTPUT OF THE 
•PROGRAM IS THE TIME FROM *0* WHEN THE POINT WAS SENSED* 

•THE NUMBER OF SWATH INCREMENTS BEFORE THE POINT WAS SENSED* 
'THE LATERAL RESOLUTION ELEMENT NUMBER* AND THE DETECTOR 
"NUMBER WITHIN THE STACK OF SIX. 


* INITIALLY* THE LOCATION OF THE POINT IS FOUND IN 

•SENSOR COORDINATES FOP. THE GIVEN S/C LOCATION. BASED 
•REQUIRED ESTIMATES* SUCCESSIVE ITERATIONS ADVANCE OR 
'THE S/C TO THE POINT OF SENSING. THE SWATH INCREMENT 
•A LINEAR VELOCITY IN SENSOR COORDINATES CONSEQUENTLY 
•MUST BE TAKEN TO INSURE THAT EXCESSIVE SWATH I NCREMEN 
•DOES NOT OCCUR SINCE IN THAT CASE TILE POINT WOULD BE 


RETREAT 
S ASSUME 
CARE 
TAT I ON 


’•SKIPPED IN THE FIRST PASS. IF THE APPROACH IS MADE 
•'GRADUAL* THE POINT WILL NOT BE SKIPPED aUD THE ACCE 
"TIME MAY DP FOUND IN ONE PASS* HOWEVER* SEVERAL PHA 
"S’UVTK, INCREMENTATION MUST BE USED. 


OF 
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330 
340 
350 
360 
370 
380 
390 
400 
410 
430 
430 
440 
453 
460 
470 
480 
490 
500 
51 0 
520 
5 30 
540 
553 
560 
570 
580 
593 
630 
610 
62 3 
633 
643 
650 
663 
670 
680 
690 


THE X-Y PLANE OF THE SENSOR COORDINATE SYSTEM NAS 
"DIVIDED INTO FIVE ZONES. IT SHOULD 3E NOTED THAT THIS 
"PROGRAM ASSUMES EITHER FORWARD SCANNING OR CIRCULAR 
"SCANNING AND CANNOT HANDLE AFT SCANNING. THE REGION WHERE 
"X>3 IS LABELED AS 1, THE REGION WHERE X<0 AND X>XMAX 
"WHERE XMAX IS THE X-COORDINATE OF THE FOREMOST ( SMALLEST 
"NEGATIVE) PORTION OF THE ACTIVE REGION OF THE SWATH IS 
"LABELLED 2, THE REGION WHERE X<XMAX AND WHERE X IS NOT CON- 
TAINED WITHIN THE ACTIVE SWATH AREA IS LABELLED REGION 3> 
"THE SWATH AREA IS LABELLED AS REGION 4, AND ALL AREA 
"BEHIND THE SWATH IS LABELLED REGION 5. 

" AT FIRST THE S/C POSTTION IS KNOWN AT SOME SWATH 

"CENTER. THE LOCATION OF THE GROUND POINT IN SENSOR 
"COORDINATES IS COMPUTED FOR THE BEGINNING AND ENDING 
"OF THAT SWATH. SINCE THE SWATH SCAN TIME IS KNOWN, 

"A LINEAR VELOCITY IN SENSOR COORDINATES MAY BE FOUND. 

"IF THE POINT FALLS IN REGION 1, THE NUMBER OF SWATH 
"INCREMENTS NECESSARY TO PLACE THE POINT JUST INSIDE REGION 
"2 IS ESTIMATED USING A LINEAR VELOCITY APPROACH. 

"THE PROGRAM RETURNS AND RECALCULATES THE TWO POINTS AGAIN 
"FOR THE NEW TIME AND S/C LOCATION. IF IT FALLS IN REGION 
"2, THE NECESSARY SWATH INCREMENT TO .BRING IT INTO REGION 
"3 IS COMPUTED AS IN ABOVE. HOWEVER, BECAUSE THE SWaTK 
"INCREMENTS ARE OVER ESTIMATED (TANGENTIAL RELATIONSHIP 
" due TO INCREASED SLANT RANGE), A SCALE FACTOR OF 3.95 
"IS USED. THE PROGRAM RETURNS AND LOCATES THE TWO POINTS FOR 
"THE NEW TIME AND S/C LOCATION. THE ABOVE PROCEDURE IS 
"FOLLOWED UNTIL THE POINTS FALL INTO REGION 3. AGAIN 
"REQUIRED SWATH INCREMENTS ARE CALCULATED, THIS TIME TO 
"FORCE THE POINTS TO FALL INTO THE ACTIVE AREA OF THE 
"SCANNER. AS A PRECAUTION AGAINST OVER INCREMENTING, THE 
"ESTIMATES ARE REDUCED BY 4 SWATHS AND THE PROGRAM RETURNS. 
"IF THE POINT IS INITIALLY IN REGION 5, THE S/C IS MOVED 
"BACK UNTIL IT FALLS IN REGION -3 AT WHICH TIME THE 
"NORMAL PROCEDURE IS FOLLOWED FROM THAT POINT. 

"ITERATION CONTINUES UNTIL THE POINT ASSOCIATED WITH THE 
"END OF ONE SWATH FALLS WITHIN THE ACTIVE AREA. 


AT THIS POINT INTERSWATH MODE IS BEGUN, THAT IS 
"SAY, THE POINT WILL BE SENSED IN THAT SWATH OR THE ON 
"FOLLOWING IT. BASED ON ITS LOCATION IN SENSOR COORDI 
"THE LATERAL RELOLUT I ON ELEMENT NUMBER IS ESTIMATED. . 
’THIS POINT, THE GROUND POINT IS PROJECTED FOR ONLY OJ 
•TIME, THAT TIME BEING ESTIMATED IN THE SAME MANNER AS 
•LATERAL RESOLUTION ELEMENT. BY PROJECTING THE POINT ] 
•VARIOUS ELEMENTS ABOUT THIS TIME, THE DETECTOR STACK 4 
•3F LOCATED. IF AFTER SEVERAL ATTEMPTS (6) , THE STaC* 
'HAS. NOT BEEN LOCATED, THE SWATH NUMBER IS INCREMENTED 


AND THE ABOVE PRODEDUR: 


REPEATED. ONCE 
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810 "IS LOCATED* THE DETECTOR NUMBER Cl THROUGH 6 ) IS FOUND 
820 "3Y A MINIMUM DISTANCE TO CENTER CRITERION. IF THE POINT IS 
830 "NEVER SENSED DUE TO EARTH ROTATION GAPS* A STATEMENT IS 
8/40 "PRINTED TO THE EFFECT THAT SWATH INCREMENTATION WAS TOO 
850 "LARGE. 

860 " THE INPUT TO THE PROGRAM FROM A DATA FILE IS THE 

870 "LATITUDE AND LONGITUDE OF THE POINT WHILE THE OUTPUT IS 
880 "THE TIME FROM THE INITIAL S/C LOCATION* THE LATERAL RLL- 
890 "SOLUTION ELEMENT NUMBER* THE DETECTOR NUMBER* AND THE 
900 "SWATH INCREMENT NUMBER. 

910 " 

920 " 

93 0 " DEFINITION OF VARIABLES 

940 " 

9 50 " 

960 " DGR- CONVERTS DEGREES TC RADIANS 

970 "PLAT=LAT I TUDE OF GROUND POINT 

980 "PLUG = LONG I TUDE OF GROUND POINT 

990 "A= SEMI -MAJOR AXIS OF THE EARTH 

10:08 "ECCN = ECCENTRI CITY OF THE EARTH 

1018 “GAMMA* CONE ANGLE OF THE SENSOR 

1098 "SETA=P I TCH USED TO PRODUCE FORWARD SCANNING 

1038 "THETA=ROLL OFFSET 

1048 “ALP® COMPLEMENT OF ORBITAL INCLINATION 
1850 "VO=TAMGEMT I AL S/C VELOCITY 
I860 “RO= RADIUS OF CIRCULAR ORBIT 
1870 "Cl = EARTH ROTATION RATE 

1288 “C3= CONSTANT USED IN. CONJUNCTION WITH ORBITAL ANGLE 

1890 "T 0=1/2 SWATH SCAN TIME 

1100 "C0= SCANNER POLAR ANGLE SPEED 

1110 “3=SEMI -MINOR AXIS OF THE EARTH 

1120 “THT® ORBITAL ANGLE MEASURED BETWEEN THE TERMINATOR AND THE 

1138 ” INSTANTANEOUS LOCATION OF THE S/C 

114* "PHISEO=PITCH AT TIME EERO 

1158 "PHI D=? ITCH RATE AT TIME EERO 

1160 "PKIDD=PITCH ACCELERATION AT TIME ZERO 

1170 "XHAPPC=YaV AT TIME ZERO 

1180 ”?:XAPD=YAV RATE AT TIME ZERO 

1198 "XKAPDD=YAV ACCELERATION AT TIME ZERO 

1288 "VO=ROLL AT TIME ZERO 

1918 D=RQLL RATE AT TIME ZERO 

1228 DD=ROLL ACCELERATION AT TIME ZERO 

1238 "TOO = T I ME BETWEEN SWATHS 

1240 "XLSE=L0N3 I TUDE OF S/C AT TIME T 

1258 “TP* SWATH CENTER TIME INITIALLY 0 

1 2 60 ”XLPK=LONG ITUDE OF ORBITAL TERMINATOR 

1270 "RAD=RADIUS OF SCAN AREA IN SENSOR COORDINATES USING A 
1280 " Z-DI MENS I ON OF 1.0 
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1290 "DEL=DISPLACEMENT OF SCAN AREA IN WHICH TWO CORCLES OF 
1300 M RADA 1 1 HAD OFFSET BY DEL WITH RESPECT TO ONE ANOTHER 

1310 " FORM THE ACTIVE SCAN AREA 

1020 "T=THE INSTANTANEOUS TIME AT WHICH THE GROUND POINT IS TO 

1330 " 3E PROJECTED 

1340 "PHI SE=LATITUDE OF S/C AT T 

1 3 r '°! ”HSE=ALTITUDE OF S/C AT T 

1360 "ALPHASHEADING OF S/C AT T 

1370 "PH I T= TIME VARYING CONTRIBUTION OF PITCH 

1380 "XKAPT=TIME VARYING CONTRIBUTION OF YAW ’ 

1390 "WT=T I ME VARYING CONTRIBUTION OF ROLL 
1420 ”PHI=PITCH AT TIME T 
1410 "V=RQLL AT TIME T 
1420 ’ * X X APP A = Y A W AT TIME T 


1430 
1 44 0 
14 50 
1 460 
147 0 
1480 
1490 
150 0 
1 51 0 
1 520 
1 530 
1540 
1 550 
1 5 60 
1 570 
1580 
1 590 
1 600 
1610 
1 620 
1630 
1640 
1650 
1 660 
1 67 0 
1 680 
1690 
1 70 0 
17 10 


”NCOUNT= THE LATERAL RELOLUT I ON ELEMENT NUMBER 
"NTST=THE TOTAL NO. OF SWATH INCREMENTS 
” I TEM=THE DETECTOR NO. 

"XT=TOP OF RESOLUTION ELEMENT STACK IN X 
"X3=3GTTOM OF RESOLUTION ELEMENT STACK IN X 
"YLsLEFT BOUND OF RES. ELE. STACK IN Y 
"YR=R I GHT BOUND OF RES. ELE. STACK IN Y 

"ETAO* ETA1 * ETA2=P0LAR ANGLES ASSOCIATED WITH THE SCANNER 
REAL*8 D( 6 ) 


REAL *8 DEL ETA* PLAT * PLNG* A* DGR* ECCN* GAMMA* BETA* THETA* % 

ALP* VO* RO* Cl * C3* TO* CO* B* TUT* XX APPO* PHI SEO* PK I D* PH I DD* XaAPD* % 

XX A? DD* VO* VD* W DD* TO 0 * XLSE* TP * XLPR* XLPK* XLPKD* XX* RAD* DEL* XnAX* % 
XMIN* PH I SE* HSE* DELT* PHI SD* ALPHA* PH I T* THO* XKAPT* WT* PH I * vJ* XKAPPA* % 
El * F2* UP * VP* WP*tJX* VX*B! I * D2 1 * 33 1 * B 1 2 * B22* B32 * B 1 3* B2 3* 333* UR* % 

7R* '.;P.*XE*YE* ZE*A1 1 * A2 1 * A3 1 * A 1 2 * A22* A32 * A 1 3* A23* A33*XP* YP* EP* % 

XP2 * YP2 * Z'?7: * X? 1 , YP 1 * ZP 1 * ETA 1 * X CAL * ETA2 * VMX*XM, BB* XA* YA* EAT* % 

T* ET AO * XT* XE * ARG * GAM 1 * AP.G2 * YL* YR, PX 1 * PX2 * PX3* PX4* PX 5* PX 6* C ON* DUl-L 
PLATE* PLUG D 

INTEGER *4 ISW* MTEST* MT* I T* KMX* LOOP* !\J COUNT* ITEM* NTST 


I SW=0 

DELETA= 1 . 92 3 E- 0 4 
MTEST=0 


.’JTRY = 0 


N T = 3 
[JPT = 0 

DGR=57. 2957795131 
CALL OPEN < 1* TOLL'* * INPUT* ) 
READ ( 1 * * > PLAT * PLNG 
PLAT = PLAT/DGR 


1720 PLNG=PLNG/DGR 
1730 A=6378383. 

1740 ECCN= .082 

1750 NTST-0 

1760 GAMMA=I 3.23/DGR 
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17 79 t ’LATD=PLAT*DGR 
1780 PLNGD=PLUG*DGTl 

1790 WRITEC 6, 503) PLATD,PLNGD - 

1800 500 FORMAT (’ THE LAT . OF THE POINT IS*,F9.5//’ THE LCNG.% 
1810 OF THE POINT IS*,F9.5/) 

1820 3ETA=1 3.23/D3R 
1830 THETA=0.0 
18/40 ALP=9./DGR 

18 50 VO= 7 400 • 

I860 20=7358288. 

1870 Cl = 360 . /(861 64 .*DGR) 

1880 C3=V0/P0 

1890 TO=0. 016042 

1900 C0=48 . /< DGR*0 .03238/1) 

1913 3=A*DSGRT< 1 . -ECCN**2 ) 

1920 THT= 60 • / DGR 
1930 THO= THT 
19/10 XKAPPO = 0 » 0 

19 53 PH I SE0 = 3 * 0 
1963 PH I D=0 • 3 
1970 PH I DD=3 » 3 
19*9 XKAPD=9.0 
1990 XXAPDD=0.3 
2030 '-'0=3.0 

20 10 '. T D=0 . 3 
2023 V 00=3.3 
2939 T 00 = 3. 040 I 35 
20/10 XLSE=90 . /DGP 
2050 TP =9.3 

2960 XLPR=DATAN ( DS I N( THT ) / ( DCOS ( THT ) * DS INC ALP ) ) ) 

2073 XLPK=XLSE-XLPR 
7080 XLFKD=XL?K*DGR 
2390 XX=XLSE* DGR 
2130 I T = 3 

2110 7AD= DS I N ( 3 AMM A ) 

2123 DEL = RAD-DS INC GAMMA-2. 6/JE-0A) 

21 33 XMAX= -RAD* DC OS ( 2 4 . / DGR ) + DEL 
71/40 XMIN = - RAD 
2150 1 CONTINUE 
2160 IFCIT.E9.3) T=TP-TO . 

2170 IF(IT.EQ.l) T=TP+TO 
2180 2 THT = THQ+ C VO/P.O > *T 

7199 XL 5 E=XLPX+ C 1 *T+ DAT AN ( DS I N C THT ) / C DCOS C THT ) * DS I • J ( ALP ) ) ) 

2203 CALL LATLI3G < THT, ALP , A, 3, ftO, ECCH, PH I SZ, H5E, DELT) 

2213 PHISD=PKISL*DCR 

2223 CALL ASECA, ECCN,ALP, C3,C1, 3, THT, RO, PHI 5E, USE, ALPHA) 

7233 ?HIT=?HI SEO+PHI D*T+PHI DD*T**2/2 
2293 XXAPT=XXAPP0+XKAPD*T+XXAPDD*T**2/2 
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22 50 VT* ! -'6+'JD*T + tfDD*T**2/2 
2260 PH I = - ( PHI T+ BETA ) 

2270 w=- (-,??+ THETA) 

2280 XKAPPA=-XKAPT 

2290 El = A/DSGRT ( 1 . -ECCN**2* D5 I N (PLAT ) **2 ) 

2300 E2= A/DSQRT ( 1 . -ECCN**2*DSIN(PHI SE) **2 ) 

2310 fJP = Et*DCOSCPLAT)*DCOSCXLSE-PLMG) 

2320 VP = E 1 * DC OS (PLAT ) * DS IN ( XL SE-PLNG ) 

2330 VP=EI*( 1 . -ECCN**2)*DSlhJ(PLAf > 

2340 UX=E2* DCOS (PHI SE) 

2350 VX=E2*( 1 .-ECCN**2)*DSINCPHISE> 

2360 B 1 1 = DCOS ( ALPHA ) * DS I D ( PH I SE) ' 

2370 B2 1 = DS ID ( ALPHA > * DS IN ( PH I SE) 

2380 B3 1 = DCOS ( PH I SE) 

2390 B 1 2- - DS ID ( ALPHA ) 

242°! B22=DCOS( ALPHA ) 

2410 332=0 • 0 

2420 Bl 3= -DCOS (ALPHA)* DCOS < PHI SE) 

2430 132 3= - DS I N ( ALPHA ) * DC 0 S ( PH I SE ) 

2440 133 3= DS I D < PH I S E ) 

2450 02= UP -UX 
2460 VH = VP 
2470 7R = f *7P-VX 

2480 XF*Bl 1 *UR+3 1 2* VR+3 1 3*WR 
2490 YE=B2] *UR+B22*VR+B23*VR 
2500 ZE=33 1 *UR+B32* VR+B33*WR 

2 510 A 1 1 = DC OS < PHI ) * DC OS ( XKAPPA ) - DS I D C ?H I > * DS I D C W ) * DS I N ( XKAPPA ) 
2520 A2 1 = -DCOS ( J ) * DSI M ( XKAPPA ) 

2530 A3 1 = DS I N (PH I ) * DC OS ( XKAPPA ) + DCOS ( PH I ) * DS I D ( W ) * DS I D ( XKAPPA ) 
2540 A1 2= DCOS (PHI ) * DS IN (XKAPPA) + DS I N (PH I )* DS I N ( V )* DCOS ( XKAPPA ) 
2 550 A22 = DC OS (N)* DCOS (XKAPPA) 

2560 A32=DSIN(PHI )*DS IN (XKAPPA) -DCOS (PHI ) *DS I N ( W ) *DCOS (XKAPPA) 
2570 A 1 3= -DS IN ( PH I )*DCOS(W) 

2580 A2 3= DS ID ( W ) 

2590 A3 3 = DC 0 S ( Pi 1 1 ) * D C 0 S ( V ) 

2600 XP=A 1 1 *XE+A1 2+YE+A1 3*<ZE-HSE) 

2610 YP=A21 *XE+A22*YE+A23*(ZE-1I5E) 

2620 ZP=A31 *XE+A32*YE+A33*(ZE-HSE) 

263 0 I F ( I T . EQ « 3 ) GO TO 3 

2640 XP2=XP 

2650 YP2=YP 

2660 ZP2=ZP 

2670 I T = 0 

2680 GO TO 4 

2690 3 XP 1 =XP 

2700 YP1=YP 

2710 ZP1=ZP 

2720 I T= 1 
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0730 IFCISV.EG.l) GO TO 18 

27/40 1FC ISV.GE.2) GO TO 25 

27 571 GO TO 1 

2760 4 XP1=-XP1/ZP1 

2770 YP 1 = -Y? 1 /ZP 1 

2780 ZP 1=1.0 

2790 XP2--XP2/ZP2 

2800 YP2=-YP2/ZP2 

28 10 ZP2=1 .0 

2gor» I F (XP2 . GT. 0 . O ) GO TO 10 
2830 IFCXP2.GT.XMAX) GO TO 11 
2840 ETAl =DATAtJ<YP2/( DEL-XP2 ) ) 

2850 XCAL= - DS IM ( GAMMA ) * DC 05 ( ETA 1 ) + DEL 
2863 I F ( XP2 ♦ G T .X CAL ) GO TO 12 

28 7 0 ETAS - DAT AM ( YP 1 / ( -XP 1 ) ) 

2880 XCAL=-DS I M (GAMMA ) * DC OS ( ETA2 ) 

2890 IF(XP1 .LT.XCAL ) GO TO 14 
2900 GO TO 17 

2910 10 VMX= (XP2-XP 1)/(2»*T0) 

2920 KMX=-XP2/( VMX*TOO> 

2938 KMX = KMX +5 

29 40 TP = TP + D FL 0 AT < KMX ) * T 0 0 
2950 NTEST=1 

2953 NT ST = NTS T + KMX 
2970 GO TO 1 

2980 11 XM= ( YP2 -Y? 1 ) / ( XP2-XP 1 ) 

2990 BB=YP2 -XM+XP2 

3030 X I NT * C DEL - XM* 8 3 - DS Q R T C D S I N ( G AM MA ) * * 2 * < 1 » +XM**2 ) - {XM*DEL + BB ) 
30 1 0 **2) )/( 1 .+XM+ + 2) 

3020 VMX=fXP2-XPl )/(2.*T0> 

3370 XMX=3 .9* (XMAX-XP2 ) / < VMX*TOO) 

3343 xmx=:c:x+i 

30 53 I F ( NT . EQ • 1 ) KMX=KMX+ 1 1 0 *KMX/ ( 1 00 *KMX ) 

3060 TP=TP+DFLOAT <XMX)*TOO 

3070 NT=MT+ 1 
3080 mtst=:jtst+kmx 
3090 N?EST=1 
3100 GO TO 1 

3110 12 V MX = ( XP 2 - XP 1 >/(2.*T0> 

3 1 20 XMX= CXCAL-XP2 ) / ( VNX*TGO ) 

3130 KMX* KMX -2 

3140 I F ( : JP T . G E . 1) X MX = 1 

3150 NPT=NPT+1 

31 60 TP=T?+DFLOAT(XMX)*TOO 

3170 MTST =WTST+KMX 

3180 MTEST* 1 

3190 GO TO I 

3200 14 I F CUTEST . EQ .0 ) GO. TO 15 
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3210 WRITE* 6, 104) 

3220 104 FORMAT** THE POINT HAS BEEN SKIPPED SOMEHOW, TOO MANY % 
3230 SWATH INCREMENTS AT SOME POINT’) 

3240 STOP 

3250 1 5 VMX= *XP2 -XP 1 )/(2.*T0> 

3260 XMX= (XMIN-XP1 )/{VMX*TOO) 

3270 HMX=KMX-50 

3230 TP=TP+DFLOAT*KMX)*TOO 

3292 NTST = NTST + KMX 

3300 GO TO 1 

3310 17 WRITE* 6, 107) 

3320 107 FORMAT*’ INTERSWATH MODE HAS BEGUN') 

3330 LO OP=0 

3340 XA=(XPl+XP2)/2. 

3350 YA= ( YP 1 +YP2 ) /2 . 

3360 EAT = DATAN < -YA/XA ) 

3373 M C OUNT = EA T / DEL ET A 
3330 NCOUNT = N COUNT- 5 

3390 T=TP+ TO* 1 . 923 E-04+DFL OAT (NCOUNT) +DGR/24 . 

3400 I SW=1 
3410 GO TO 2 

3420 18 EAT= DATAN ( -YP 1 /XP 1 ) 

3430 ISW=ISW+1 
3440 N C OUNT = EAT /DEL ETA 
34^0 NCOUNT = NC0 T JNT-2 
3460 25 XPt =-XP ’ /ZPl 
3473 YP 1 =-YP 1 /ZV 1 
3480 LOOP=LOOP+l 
3490 IT=3 

3503 ETAO= DFLOAT* NCOUNT)* DELETA 
3510 XT= - DS IN ( GAMMA ) *DCOS ( ETAO ) 

3523 X3-XT+ DEL 

3530 AU3 = DCOS* GAMMA ) * ( DC OS ( 2 . 2 D- 0 5 ) - DS I N ( GAMMA ) * DS I N * ETAO ) * % 

3543 DS I M ( 2 » 2 D- 0 5 ) / DS Q E T < 1 . -DS IN ( ETAO ) **2*DS I N ( GAMMA ) **2 ) ) 

3550 GAM1 = DATAN * DSQRT ( 1 . -ARG**2 ) / ARG ) 

3 5 60 ARG 2 = DS I N ( GAMMA ) * DC 0 S * G AM 1 ) * D C 0 S ( ET AO ) / < DC 0 S ( GAMMA ) * DS I N % 
3570 ( GAM 1 ) ) 

358 0 ETA 1 = DATAN * DSQRT * 1 . -ARG 2** 2 ) /ARG 2 ) 

3590 YL = DSIN* GAM1 )*DSIM(ETA.l ) 

3600 IF*NCOUNT.LT.0) YL=-YL 

3610 ARG= DCOS * GAMMA ) * * DCOS ( 2 . 2 D-0 5 ) +DS I N ( GAMMA ) * DS I N ( ETAO > * % 

3620 DS I H < 2 . 2 D-0 5 > / DSQT1T < 1 . - DS I N * ETAO ) * *2*DS I M * GAMMA )**2) ) 

3630 G AMI = DATAN * DSQRT * 1 . -ARG**2 ) /ARG ) 

3 640 ARG2 = DS I N ( GAMMA ) * DC OS ( GAM 1 ) * DCC 5 ( ETAO )/ * DCOS * GAMMA ) * U 3 1 N ( % 
3650 GAM 1 ) ) 

3667 ETA 1 =DATAN( DSQRT ( 1 . -AEG2**2 ) /ARG2 ) 

3670 YR=DS I M * GAM 1 )*DSIU*ETAl ) 

3680 I F (NCOUNT .LT .0 ) YR=* -YR 



3690 IFCNCOUNT.SE. 0) GO TO 42 

3700 40 DUM»YL 

3710 YL=YR 

3720 Y r? = DtJM 

3730 42 CONTINUE 

3740 IFCXP1 .LT.XB) GO TO 19 

3750 T=T+T0*1 .923D-04*DGR/24. 

3760 MCOUNT*NCOUNT+1 
3770 GO TO 24 

3780 19 I F CXP 1 .GT*XT) GO TO 20 
3790 7*1 TEC 6, 1 1 1 > 

3800 111 FORMAT (' THE POINT IS HIGH, IT WILL NOT BE SENSED 
3810 STOP 

3890 20 IFCYP1 .LT.YR) GO TO 22 
3830 T=T+TO* 1 • 923D-04+DGE/24. 

3840 NC0UNT=NC0UNT+1 
3850 GO TO 24 

3860 22 IFCYP1 .GT.YL) GO TO 23 
38' 7 0 SI" I TEC 6, 1 1 3 ) 

3860 113 FORMAT < • THE POINT IS TO THE LEFT') 

3890 STOP 
3900 23 CONTINUE 
3910 PX 1=-(XT-XB)/12.+XT 
3990 PX2= - (XT-X3J/4. +XT 
3930 PX3= - CXT-XB) *5 . / 1 2 .+XT 
39/40 PX4=- (XT-XS) *7 ./I 2 . +XT 

3950 PX5=-(XT-XB)*3./4.+XT 
3960 PX6= - (XT-X3 ) * 1 1 ./12.+XT 
3970 DC 1 >=DA3S(XP1 -PX1 ) 

398 0 DC 2 ) =DA3S CXP 1 -PX2 > 

3990 DC 3 ) = DADS CXP 1 -PX3 ) 

4003 DC 4 ) =DABS C XP 1 -PX4 ) 

4010 DC 5) =DABS CXP 1 -PX5) 

4020 DC 6 ) s DABS CXP 1 -PX6> 

4030 CON=D( 1 ) 

40 40 I T EM= 1 

4350 DO 31 1 = 2,6 

4060 I FCCON.LT. DC I ) ) GO TO 31 

4378 CON-D(I) 

48 8 0 ITE1T=1 
4893 31 CONTINUE 

4198 ’.7* I TEC 6* 1 1 5 ) TP, NCOUNT, I TEX 

4113 115 FORMAT C * THE FINAL 'TIME IS’,F9.5, ' THE ELEMENT NO. 
4120 15/' THE ELEMENT WITHIN THE STACK IS', 15) 

4138 UP I TEC 6, 40 3) NT ST 

4148 403 FORMAT C* THE TOTAL NO. OF SWATH IN CHE. IS ’,16/) 

4150 STOP 

4160 24 CONTINUE 


' ) 


I S ' , % 
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41725 IF (LOOP • 3E » 6 ) GO TO 26 

41 80 GO TO 2 

4190 26 CONTINUE 

4200 ISW=0 

4210 TP=TP+TOO 

4220 IF(NTRY.EQ .2 ) STOP 

4230 MTRY=NTRY+1 

4240 GO TO 1 

4250 END 

42 60 SUBROUTINE LATLNG ( THT, ALP, A, B, RO, ECCN, PHI SE, HSE, BELT ) 

4270 PEAL *8 THT, ALP, A, B, RO, ECCN, PH I SE, HSE, DELT, HSEH 

4280 MS EM=9 7990 3 • 

42°0 UCT=1 

4303 1 PHI 5E=DATAU( ( <A**2+A*HSEitf>/ (B**2+A*HSEN> ) * < DCOS ( TUT ) *ECOS% 
43 1 0 ( ALP ) / DSGPT ( DS I N ( TKT ) **2+ DC OS ( THT ) **2 * DS I N ( ALP > **2 ) ) ) 

4320 HSE=PO*DCOS(THT)*DCOS(ALP>/DSIiM(PHlSE) -A*< 1 . -LCCJ**2 ) /D5GaT% 
4330 ( 1 . -ECCN**2*DSIM(PHISE>**2> 

4340 DELT= DABS (HSEM-H3E) 

4350 . IFCDELT-l .0) 2,2,13 

4360 13 IF(NCT-20) 14,2,2 

4370 14 HSEN=HSE 

4380 M C T = N C T + 1 

4390 GO TO 1 

4403 2 RETURN 

4410 END 

4423 SUBROUTINE AS E ( A , E C C N, AL? , C 1 , C2 , B, TH T , R 0 , PH I S E, 1 i S E , AZ ) 

4430 REAL*8 XO, YO,ZO, A, ECCN,ALP,Cl, C2, B, THT, HO, PHI SE, H SE, AZ, % 

4443 EPP, PHO, X, Y, Z, UARG, C 1 1 , C22, C 1 2 , C2 l , C 1 3, C2 3, PHI SET, PHOT, % 

4450 EPPT , XT , YT, ZT, CCC1 , CCC2 , ZETA, X ETA, YETA, P 1 , P2, P3 , DOT 
4460 XO=RO*DSlN<THT> 

4470 YO=RO* DCOS ( THT ) * DS I N ( ALP ) 

4480 ZO=RO*DCOS ( TKT ) * DCOS C ALP) 

4493 EPP=A/DSGRT( 1 . -ECCM**2*DSIM(PKISE>**2) 

45 ’^ ts HO = DATAN< B5 I M ( THT ) / ( DCOS ( THT ) *DS I M ( ALP > ) ) 

43 l 3 X= EPP* DCOS (Pit ISE)*C3I N (PHO ) 

4 523 Y= EPP * DC OS ( PH I S E ) * DC 0 S ( PH 0 > 

45 33 Z = EPP* ( 1 . -EC CM* * 2 > * D$ I M ( PH ISE) 

4 540 UARG=C CA**2+A*HSE) /(3**2+A*HSE> > * < DCOS (TiiT)*DCG3 <ALi' ) /BSN:IT% 
4553 ( DS IN ( THT )**2+DC0S( THT) **2*DS I N (ALP ) **2 ) > 

4560 Cl 1 = 1 .3 
4570 C22= 1 ♦ 0 

4 58 3 C 1 2=R0*DC0S ( THT ) * DC OS < ALP ) * LCDS C?hl SE) /DSIN ( PH ISE) **2+A* ( Z 
459-3 1 • - EC CM* *2 ) * EC CM* EC CM** 2 * DS IN (PH I SE ) * DCOS (PH l SE) / ( DS’4RT( 1 . - % 

4633 ECCN**2*D5 IN (PH ISE)** 2) )**3 

4613 C21=(l»/(1. + UARG* *2 ) ) * ( A* ( A* *2 -3**2 ) / ( 3**2 + a SE* A) * * 2 ) 

4623 C 1 3= -RO* DS I N ( THT ) *DCOS ( AL? ) / DS I i ) (PH I SE ) 

4630 C2 3= < (A**2 + A*HSE)/(B**2+A*HSE) ) * ( 1 ./Cl . + UARG**2 ) ) * (-DCOS % 

4643 (ALP > * DS IN < THT ) / ( DSQRT ( DSIN ( THT )**2+DC0S(THT)**2*DSIM (ALP ) % 
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4650 **2>>**3> 

4660 PHISET=( (Cl 3-C23*Cl 1 /C2I )/<C12-C22*Cl 1/C21 ) )*C1 

467 0 PHOT=C 1* DSI M ( ALP ) / (. DC OS ( THT ) **2* DS I N (ALP ) **2 + DS I N ( THT ) % 

4680 **2) 

4690 EP?T=A*ECCN**2*DSIN< PHI SE)*DCOS( PHI 5E) *PH1 SET/DSiiPTS 
4780 ( ( 1 .-ECCN**2*DSIN<PHI SE)**2)**3> 

/■71 - YT = PCOS ( PK I 5E ) *DS I N ( PIIO ) * EPPT+ EPP* DC 0 5 ( PH I 5E ) * DCOS ( PH 3 ) *PiiOT% 
4720 -epp*dsi:;(pho>*dsi:jcphise)*phiset 

47 38 YT=DC03 ( PHI 5 E) *DCOS ( PliO ) * EPP T - EPP * DC 05 ( PH I SE) * DS I N ( PHD ) *PiiGT % 
4740 - EPP* DC 0 S < ?H0 ) * DS I N ( ?H I S E ) *PH I SET 

47 50 Z 7= ( l . -ECC>{**2)*<BSIN<PHISE)* EPPT + EPP* DC OS (PHI SE)*P*»I SET ) 

4760 CCC 1 = - (A* *2 /B**2 ) *Z / (X*DS I M ( PHO ) +Y* DCOS ( PHO ) ) 

4770 CCC2=DSQHT ( 1 . / ( 1 • + CCC1 **2 ) ) 

4780 Z ETA=- CCC2 

4798 XETA=-CCC1 *DS I M (PHO > *CCC2 
4888 YETA- - CCC 1 * DCOS ( PHO ) * CCC2 
4810 P 1 -XE7A+XT 
4820 ?2=YETA*YT 
4830 P3=ZETA*ZT 

4840 DOT= (PI +P2+P3 ) /DSQPT ( ( KT**2+YT**2+ZT**2 ) * (XETA**2+YETA**2+% 
4858 ZETA**2 ) ) 

48 68 AZ = DAT AM ( DSQRT (1 .-BO 7**2)/ DOT) 

4870 PE TURN 

4380 END 
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APPENDIX B 

USER REQUIREMENT SURVEY 
R. J. Ouellette 


B.l INTRODUCTION 

The purpose of the user requirement survey was to determine the 
ground resolution— cell sizes required for the different anticipated 
observation needs of major users. A second purpose was to provide a 
unified base of data from which sensor resolution, registration, and 
positioning requirements can be determined in the future. 

The survey first determined what observations are of interest to 
users. The user group consisted of ten agencies in four government 
departments and represented the majority of government users. The sur- 
vey consisted of telephone interviews with fourteen technical personnel 
within these four departments. These interviews are listed as entries 
in the list of sources at the end of this section. The observations 
given by this user group were classified by basic shape and size. From 
this shape and size information an estimate of the required ground 
resolution cell size was derived. The views expressed in this section 
summarize the information received from the telephone conversations. 

In the summarization some amount of subjectivity has inevitably been 
introduced. It is hoped that as little misinterpretation as possible 
has occurred. 

In Section B.2 user observations are grouped into three major 
categories: land use, hydrology/geology, and cartography. Each of 

these categories is discussed in detail in terms of specific observa- 
tions and their respective ground sizes and shapes. 

In Section B.3 ground observations are categorized by shape and 
dimension, and a summary table is given with estimates of the percent- 
age of the total observations each categorized group represents. Also 
the size and shape groups are related to required ground resolution 
cell sizes for these observations. 

As a preface to the rest of this study two preliminary considera- 
tions merit discussion here. The first is related to the fact that 
there has been a basic lack of satellite imagery available to users to 
evaluate. This lack of experience with the raw data sometimes results 
in the users not knowing what is required of satellite imagery for a 
desired observation. This lack of experience is also the main reason 
for this study considering observations in terms of ground size and 
shape since these are values readily known by the user. 
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The second consideration relates to the types of users of satel- 
lite images. One type of user is interested in high resolution imagery 
over a large area but does not need multiple coverage. Another user is 
interested in repetitive coverage over large areas but does not need 
high resolution. Still a third type of user needs high resolution and 
repetitive coverage but only over small areas o The second group is the 
major user group. 

The second type of user generally represents the largest portion 
of the user group and is best suited to having his observation needs 
met by satellite imagery. The two other types of users can have some 
of their observations needs met by satellite imagery, and thus can pro- 
fitably participate in the satellite program by sharing image data with 
the second group. 

The most significant outcome of the survey is that the great 
majority of the user observational needs would be met with a 30— to— 70 
meter ground resolution cell size. To have a significant further im- 
pact on the user observation needs, a ground resolution cell size less 
than 10 meters would be required, (Ground resolution cell size as used 
here refers to the linear dimension on the ground defined, in a pure 
geometric sense, by the resolution element of the sensor.) 

B .2 USER OBSERVATION NEEDS 

B.2.1 Background 

United States Government agencies interested in aerial 
image data were contacted to evaluate their observation needs for pos- 
sible satellite image data. These government agencies include the 
Department of Interior, Department of Agriculture, Department of Com- 
merce, and Army Corps of Engineers. 

The Departments of Agriculture and Interior will probably 
be the largest users within this group. Of these two. The Department 
of Interior has the largest amount of equipment with which to analyze 
raw satellite data. The Department of Agriculture is working with the 
-NASA/HSC Science and Applications Directorate, Houston, Texas, to 
develop the data analysis techniques for the observations that the 
Department is interested in. 

Major government users regard satellite image data as 
part of a multi-layered data acquisition system. Satellite data re- 
presents the coarsest upper layer of data while aircraft imagery and 
ground data acquisition represent the medium and fine resolution layers 
of the system. 

Three qualifying statements should be made before pre- 
senting the data of this section. The first is that this survey re- 
stricts itself to land-based as opposed to ocean-based, observational 
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needs, because land-based observations have the most stringent resolu- 
tion registration and positioning requirements. 

The second is that for this study observations that re- 
quire size, shape, texture , and shadow analyses of objects 10 meters or 
less in diameter have not been considered. Satellite imagery of higher 
resolution (10 meters diameter or less) is impractical for large-area 
coverage because of the high data volume. If high resolution image 
data were acquired over the whole earth surface, image data volumes 
would be extremely high. For example, the continental United States, 
mapped at 10 meter square ground areas, would yield 6.6 x lQll bits of 
image data (6 bit digitization) . At anticipated satellite mapping 
rates this surface would be remapped every 18 days. There is an in- 
creased computer cost with an increase in image volume. Data rates in- 
crease as the inverse square of the decrease in resolution element size. 
Very high resolution pointable sensors would solve some of the problems 
of large amounts of image data but would make it more difficult to 
cover small specific areas of interest situated throughout a large 
ground surface area. 

The third qualifying statement is that only the geometric 
measurements involved in spectral signature analysis were considered 
for this survey. Spectral signature analysis is the use of radiometric 
signals from a surface separated into several spectral regions to define 
surface characteristics, such as crop species identification. The 
radiometric aspects of signature analysis were not considered a part of 
this task. Geometric aspects were used to the extent they were needed 
for the detection of spot, linear and area features by temporal and 
spatial signature differences. By temporal signature difference is 
meant change detection of the radiometric surface signal from one satel- 
lite pass to another, over the same area. Spatial signature difference 
refers to difference between a particular area and the surrounding 
area. u 

Observational needs are considered in three basic areas: 
land use, hydrology/geology, and cartography. These basic areas and 
their sub-groups are given in Table B-l. Within each area the observa- 
tions of interest are discussed below, and size and shape information 
are presented. Observations will be categorized into one or more of 
the three basic ground shapes of spot, linear, and area features. 


B.2,2 Land Use 


General 

Land use involves the specific areas of agriculture 
forestry , .and private/public/commercial land development* A general ob- 
servational need common to each of these specific areas is boundary 
delineation* 


B-3 



Table B-l - Observational Need Categories 


1* Land Use 

A* Agriculture 

B. Forestry 

C. Private/Public/Commercial 

2. Hyd ro logy / Geology 

A. Control 

B . Survey 

3. Cartography 

Boundary delineation of the following interfaces are of 
interest : grass/brush/ timberland/ desert/water/crop/bare soil/ culturally 

developed. Some interface boundaries are sharp, due to inherent dif- 
ferences between the areas (e.g., water/grass interface) or the pres- 
ence of linear objects at the boundaries (e.g», rivers, roads, or tree 
rows) * Other interface boundaries are less well defined because of 
smaller inherent differences between the areas (e.g., boundaries between 
different crop types). 

Boundary delineation involves linear and area features. 

For these observations linear features have ground dimensions from ten 
meters (tree rows, small rivers) to fifty meters (large rivers, major 
highways). In cases where surface differences (spatial and spectral) 
define the boundaries the boundary width will appear to be at least 
several meters wide in the image data. 

Agriculture 

Observational needs in the area of Agriculture include: 
inventory, crop species identification, stress analysis, and soil sur- 
vey (both type and moisture content). The Department of Agriculture 
is heavily involved in crop inventory and species identification through 
the Agricultural Stabilization and Conservation Service. Species iden- 
tification measurements include surface crop identification, surface 
area measurement, and crop species identification. 

Crop fields range in size from four tenths of a hectare 
(one acre) while ownership plots range from four hectares to hundreds 
or even thousands of hectares. Average field size varies from area to 
area within the country and is dependent upon the type of crops planted. 
High-value crops such as fruits and vegetables tend to be planted in 
smaller fields, while grain crops are planted in larger fields. Fields 
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in coastal and heavily populated regions tend to be smaller and dis- 
continuous, while fields in the Midwest tend to be larger and more 
grouped together. 

Stress analysis detects infestation in plant life. Spec- 
tra signature changes are the basic data measurements. Surface areas 
o interest vary from single plant areas to entire crop fields. This 
corresponds to spot features of a fraction of a meter in diameter to 
area features as large as few tens of hectares. Larger area features 
represent the major portion of stress analysis observations. 

Soil surveys are carried out in large part by the Soil 
Conservation Service of the Department of Agriculture. Soil survey 
measurements (for the type of analysis that can be provided by remote 
imagery) involve identification of soils in terms of surface area, loca- 
tion, type, and moisture content. Presently, black and white aerial 
photos are used as base maps and for delineation of different soil sur- 
face areas. Direct ground samples are the basic sources for present 
soil classification. Area features for soil survey observations have 
sizes that range from a fraction of a hectare to tens of hectares. The 
distribution of observations through this range is fairly balanced. 

Forestry 

Observational needs for forestry include tree inventory 
species identification, and stress analysis. Forest inventory and 
species identification are done within the Department of Agriculture 
The Forest Service is involved with Federal forest lands whereas the 
Soil Conservation Service is involved with farmland forest management. 

Measurement for forest inventory and species identifica- 
tion include surface area delineation and area measurement. Species 
identification from space imagery can be done by signature analysis. 
Sizes of the areas of interest for tree inventories will range from 
four tenths of a hectare (one acre) to hundreds of hectares. Observa- 
tions involving Federal forest lands have area sizes generally larger 
than observations involving farmland forest management. Also Federal 
forest land observations represent the majority of the total 
observations. 

. Stress analysis will require smaller analysis areas than 
species identification. One of the specific questions to be answered 
by the first satellite imagery will be: How small an infestation area 

can be detected from space? When the infestation areas are small, low- 
contrast spot detection from space will be required. Infestation areas 
range from three meters to three hundred meters in diameter. The 
majority of the observations occur in the middle of this size range. 
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Private /Pub lie/ Commercial 


A third land area concerns private/public/commercial 
development , All of the government agencies contacted have some inter- 
est in this area* Observational needs are transportation and communi- 
cation route studies which require inventory and development informa- 
tion, Other needs include private and commercial land use and develop- 
ment studies which require change detection analysis in large part. 

Data measurements here include spot and linear feature detection and 
surface area measurement. 

Spot features range in size from twenty to a few hundred 
meters in diameter with the major number of observations occurring in 
the lower part of the range. Linear features for transportation and 
communication route studies range in width from five to seventy meters 
with the majority of the observations falling in the mid to upper part 
of this range. Area sizes of interest from four tenths of a hectare 
(one acre) to a few tens of hectares with the majority of the observa- 
tions grouped in the lower half of this size range, 

B, 2, 3 Hydro logy /Geology 

Flood Control 

Government agencies interested in flood control include the 
Army Corp of Engineers and the Department of Agriculture, Flood con- 
trol involves dam construction, snow studies, flood plain delineation, 
and ice jam detection. Observation needs in this area are water basin 
information for dam construction, snow cover change detection to pre- 
dict water run off, and ice jam detection for flood control. 

Spot, linear, and area feature measurements are needed 
for general control studies in the areas of hydrology /geology. Spot 
features range from three to a few hundred meters in diameter. Linear 
features range from five to a few hundred meters in width. The lower 
end of these size ranges is related to small water basin and small dam 
construction on farmland while the upper end of the size range relates 
to large water basins and country and state dam construction. Area 
features range in size from a tenth of a hectare (one quarter acre) , 
for example, for small water basins, to a few hundred hectares for 
snow cover studies. Spot and linear feature observations are dispersed 
uniformly throughout respective ranges. The majority of the area fea- 
ture observations occur in the lower part of the size range. 

Inventory Survey 

Land surface survey consists of inventory and thematic 
mapping. Inventory survey is directed at answering questions such as 
follows. How much of a particular measured quantity is present (e,g,, 
water surface survey)? What is the observed change in a measured 


B-6 



quantity (e.g., soil erosion survey)? Thematic mapping is more involved 
with the surface location of some measured quantity (e.g., geological 
feature mapping) . 

Inventory surveys measure ground water, erosion, sedi- 
mentary patterns, water pollution, mineral/petroleum, heavy metal, con- 
struction material, and natural disaster phenomena. Thematic mapping 

includes geological features, river basins, water-sheds, and coastal 
areas. 

Data needs for hydrology/geology surveys include spot, 
linear, and area features. Spot feature sizes range from about ten 
meters to one hundred meters. The great majority of measurements fall 
in the upper part of this size range. Geologic features represent the 
major observations. Linear feature widths range from five to one 
hundred meters. The majority of the measurements fall in the upper 
part of this size range. Geologic features and linear boundaries com- 
prise the major observations. Area features, including spectral signa- 
ture ground surface areas, fall in the size range from four tenths of 
a hectare (one acre) to over one hundred hectares. The largest number 
of measurements occur in the upper half of this range. The majority 
of the observations involve ground water, medium to broad geological 
structures, and material location. 

B.2.4 Cartography 

The Department of Interior is the government agency most 
interested in the use of earth resources satellite imagery to make or 
revise standard topographic maps. These maps will be used as data 
bases for most large-area remote sensing observations. 

At the present the United States is completely mapped 
at 1/250,000 scale. There is considerable interest among the different 
users for satellite imagery at 1/250,000 scale. This is interesting, 
in light of some of the high resolutions and small dimensions stated 
to be of interest. Ideally, the satellite imagery should be made 
available at the largest scale consistent with resolution and the scales 
of data base materials to be used for comparison. Information avail- 
able from well-detailed imagery at 1/1,000,000 scale is believed -to be 
appropriate for the 1/250,000 scale maps in areas that are poorly 
mapped. 

Table B-2 lists the ground surface widths of a thin line 
(0.1 mm) drawn on maps of different map scales. These ground dimen- 
sions are not necessarily the minimum feature widths to be shown on the 
maps, as can be seen from the table. Some care is needed in attempting 
to relate image resolution, map symbol size, and importance of the fea- 
ture for a particular map. If the image is substituted for a line map, 
the higher information content of the image must be made available to 
the map user. Linear features of smaller widths can be seen in the 
Image because of the eye's inherent ability to pick out line features. 
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Table B-2 - Map Scale Versus Projected Ground Size 
For A 0.1-MM Line Brawn On A Map of That Scale 


Map Scale 

Ground Size 
(m) 

1/1,000,000 

100 

1/250,000 

25 

1/100,000 

10 

1/62,500 

6 

1/50,000 

5 

1/24,000 

2.4 


Area measurements can be made where minimum breadth is only a few times 
the minimum spot resolution* 

B.3 SUMMARY 

Measurement Needs of Users 

Satellite imagery for land surface measurements will be used for 
the specific measurement categories of spot, linear, and area features. 

Table B-3 lists the three classes of features and the ground 
dimension ranges associated with fine, medium, and coarse observations 
of these features. For a particular feature type an estimate is given 
of the percentage of total observation each of the groups represents* 

An estimate is given of the required ground resolution element cell 
size range for each group. The numbers in this table are not absolute 
values, but trends indicated by the stated observational needs. Now 
consider the size class groupings in Table B-3 in terms of the major 
observations that make up each size group* For spot features the 
major fine size observation is agricultural stress measurement* Medium 
size spot features are made up of a broad set of observations covering 
each of the basic observational need areas covered in Section B.2* The 
same is true of the observational make-up of coarse size spot features. 
The two sub-groups within the medium size group represent agricultural 
stress analysis, private/ commercial land use, and hydrology /geology 
control observations for the smaller sub-group and forest stress analy- 
sis and hydrology/geology survey observations for the larger size sub- 
group. The great majority of spot feature measurements are represented 
by the medium and coarse size classes. 

Fine size linear features are made up of observations from pub- 
lic and commercial land use and the hydrology/geology flood control 
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Table B-3 - Ground Surface Dimensions, Percentage of Total 
Number of Observations, and Required Ground Resolution 
Element Range for Spot, Linear, and Area Features 


Size Classes 


Feature 

Type 

Fine 

Medium 

Coarse 

Spot 

Dimension (m) 

< 3 

3-10 

10 - 30 

> 30 


GRE Range* 

A 

A - B 

A - B - C 

B - C 


Percent of Total 

5 

15 

30 

50 

Linear 

Width 

< 10 

10 

- 30 

> 30 


GRE Range 

A 

B 

- c 

> C 


Percent of Total 

10 


70 

20 

Area 

Dimension on a Side (m) 

< 65 

65 - 200 

200 - 400 

> 400 


GRE Range 

A 

B - C 

C 

> C 


Percent of Total 

5 

25 

35 

35 


^Ground Resolution Element Range 

(dimensions given are for one side of a square ground cell) 

A B C 

< 20 m 20 m to 40 m 40 m to 70 m 


and inventory survey measurements. The great majority of linear feature 
measurements are represented by the medium size class. Boundary deline- 
ation and communication and transportation route studies represent the 
major observations. The majority of the observations for coarse size 
linear features are hydrology/geology flood control and inventory sur- 
vey measurements. 

Observations of fine-size area features include hydrology/geology 
flood control measurements, and agriculture and forestry stress analysis. 
The great majority of area feature measurements are represented by the 
medium size class. Observations in this size class are a broad repre- 
sentation of each basic observation area. The breakdown of the medium 
size group into two sub-groups represents the separation between small 
crop fields, farmland forest management, and public/commercial obser- 
vations for the smaller size sub-group, and large crop field. Federal 
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forest management, and hydrology/geology flood control observations for 
the larger size sub-group. Coarse size area feature observations cover 
the entire set of basic observation areas but represent a smaller part 
of the total number of area feature measurements than the medium size 
class. 

The final discussion in this subsection will concern itself with 
the relationship between the three types of data measurement features 
(spot, linear, and area) and ground resolution cell sizes. 

Spot feature measurements made on satellite imagery can include 
spot features whose surface dimensions are smaller than a ground- resolu- 
tion-element cell size (e,g,, ice jam detection in inland waterways for 
flood control) , Since the spot feature is smaller in dimension than 
the ground resolution cell size, the inherent ground contrast of the 
spot feature is degraded by the fact that the spot energy is averaged 
over the ground area of the resolution cell. Spot features which are 
several times smaller than the ground resolution cell will be detect- 
able only when the spot feature has a very high ground contrast. 

Linear feature detection involves location of features whose 
short dimensions can be several times smaller than a ground resolution 
element. The width of the linear feature can be measured when its 
dimension is equal to or greater than a ground resolution element. De- 
tectability of a linear feature whose short dimension is smaller than a 
ground resolution element will be determined by its resultant image 
plane contrast ratio. 

For linear boundaries between well defined areas the boundary 
location will be no worse than one ground resolution cell size. Linear 
boundaries between areas with little surface or area signal differences, 
can have an error in locating the linear boundary, of several ground 
cell sizes. 

Recent studies indicate that linear feature detection for roadway 
analysis is completely reliably only when the width of the roadway is 
at least 1/3 of the ground resolution cell (Source 17, p. 855). How- 
ever, Apollo space photography has detected roadways that were 1/10 
ground resolution cell size (Source 6, p. 33). Thus very high contrast 
linear features whose widths are as much as one tenth the ground resolu- 
tion cell size are detectable from satellite imagery. 

Measurement of ground surface areas requires that the area be re- 
solved by the sensor. Areas are delineated by linear boundary detection 
and ground surface differences such as surface structure and spectral 
signature. Signature analysis of large ground surface areas is rela- 
tively independent of the ground resolution cell size. Spectral signa- 
ture identification of large area agricultural crops is the same for 66 
to 100 meter ground resolution images as it is for 3 to 6 meter ultra- 
high resolution images. Spectral signatures of large area crop are 
distinct from one another and therefore lend themselves readily to multi- 
spectral signature analysis. (Source 11, p. 1193). 
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In general for spot, linear, and area features, if the surface 
area of a feature is not resolved, dimensional measurement cannot be 
made. It is still possible to detect the presence of the feature if 
the resultant image contrast is enough to make the feature discernible 
from the background. It should be noted that area measurements near 
the resolution limits are prone to a larger percentage of error since 
they require near-perfect transmitting and collection media. Also a 
general comment to spot, linear, and area features: recent studies 

indicate that gross geological area features, fine linear features, and 
high contrast spot features can be identified from satellite imagery 
which has a ground resolution size from 50 to 100 meters (Source 2, 
p. 499). 

B.4 LIST OF SOURCES CONSULTED FOR USER REQUIREMENT SURVEY 

1. Addess, S., Topographic Division, U. S. Geological Survey, 

U. S. Department of Interior, Goddard Space Flight Center, 
Greenbelt, Maryland, (Telephone conversation 5/15/72). 

2. Amsbury, D. L., "Geological Comparison of Spacecraft and 
Aircraft Photographs of the Potrillo Mountains, New Mexico, 
and Franklin Mountains, Texas," Proceedings of the Sixth 
International Symposium on Remote Sensing of Environment , 
Volume 1, pp. 493-501, (October 13-16, 1969). 

3. Eorgason, W., Topographic Division U. S. Geological Survey, 
U. S. Department of Interior, McLean, Virginia, (Telephone 
conversation 5/15/72) . 

4. Clifton, J. W., Agricultural Stabilization and Conservation 
Service, U. S. Department of Agriculture, Washington, D. C. , 
(Telephone conversation 5/17/72). 

5* Colvocoresses , A. P., "Image Resolutions for ERTS, SKYLAB 
and GEMINI /APOLLO," Photo gramme trie Engineering , pp. 33-35, 
(January, 1972) . 

6. Fischetti, T. L. and Bondeen, W. R. , "Earth Observatory 
Satellite Mission Review Group (EOSMRG) NASA Washington, 

D. C. , Final Report (1971). 

7. Fisher, W. A., "Projected Uses of Observations of the Earth 

from Space The EROS Program of the Department of the 

Interior," Proceedings of AIAA Earth Resources Observation 
and Information Systems Meeting , (March 2-4, 1970). 

8. Heller, R. C., Forest Service, U. S. Department of Agricul- 
ture, Berkely, California, (Telephone conversation 5/18/72). 

9. Goddard Space Flight Center, "Earth Observatory Satellite 
(EOS) Definition Phase Report," Greenbelt, Maryland, 

(August, 1971). 
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10. Jarmon, J., Civil Works Directorate, U. S. Army Corp of 
Engineers, Washington, D. C. (Telephone conversation 5/19/72). 

11. Lauer, D. T. and Thamon, R. R. , "Information Content of 
Simulated Space Photographs as a Function of Various Levels 
of Image Resolution," Proceedings of the 7th International 
Symposium on Remote Sensing of Environment , Volume II, 

pp. 1191-1204, (May 17-21, 1971). 

12. McLaurin, J., Topographic Division, U. S. Geological Survey, 

U. S. Department of Interior, McLean, Virginia, (Telephone 
conversation 5/11/72). 

.13. Miller, R. H., Research Department, U. S. Department of 
Agriculture, Washington, D. C., (Telephone conversation 
5/9/72). 

14. Mohr, J. E., Van der Meer, C., and Krishnaununni , K., 
"Evaluation of Hyperaltitude Photography for Geological 
Mapping," Proceedings of the 7th International Symposium 
on Remote Sensing of Environment , Volume III, pp. 2307-18, 
(May 17-21, 1971). 

15. -Orvedal, A., Soil Conversation Service, U. S. Department of 
Agriculture, Washington, D. C. , (Telephone conversation 
5/19/72). 

16. Scheps, B. B., Engineer Topographic Laboratory, U. S. Army 
Corp of Engineers, Fort Belvair, Virginia, (Telephone 
conservation 5/18/72). 

17. Simonett, D. S., Henderson, F. M. , and Egbert, D. D. , "On 

the Use of Space Photography for Identifying Transportation 
Routes: A Summary of Problems," Proceedings of Sixth 

International Symposium on Remote Sensing of Environment , 
Volume II, pp. 855-877, (October 13-16, 1969). 

18. Smith, C. R. and Wood, P., "Multi-Discipline Applications 
for Earth Resources Satellite Data," Presented at AIAA 
Earth Resources Observations and Information Systems 
Meeting, Annapolis, Maryland, (March 2-4, 1970). 

19. Tueller, P. T. and Garwin, L., "Environmental Analysis of 
the Lake Tahoe Basin from Small Scale Multispectral Imagery," 
Proceedings of the 7th International Symposium on Remote 
Sensing of the Environment , Volume I, pp. 453-63, (May 17-21, 
1971) . 

20. Von Steen, D. H., Statistical Reporting Service, U. S. 
Department of Agriculture, Washington, D. C., (Telephone 
conversation 5/16/72). 
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21. Weinstein, 0., Miller, B. P., and Barletta, J., "Simulation 
of ERTS RBV Imagery," Proceedings of 7th International 
Symposium on Remote Sensing of the Environment , Volume II, 
pp. 1177-89, (May 17-21, 1971). 

22. Wiegand, C., Agricultural Research Service, U. S. Department 
of Agricultural, Wesaco, Texas, (Telephone conversation 
5/10/72) . 

23. Wiesnet, C., National Environment Satellite Service, U. S. 
Department of Commerce, Suitland, Maryland, (Telephone 
conversation 5/11/72) . 

24. Wood, P., "User Requirements for Earth Resources Satellite 
Data," Proceedings of Electronic and Aerospace Systems 
Convention (EASCON) , Washington, D. C., pp. 14-20, (October 
26-28, 1970). 
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